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20.  Abstract  (continued) 

-  may  interact  with  the  electromagnetic  fields  to  yield  a  coupled  problem. 

This  thesis  presents  an  integrated  study  of  magnetoelasticity  problems  for 
thin  nonferrous  conducting  plates.  Three  phases  of  the  study  have  been 
emphasized:  theoretical  modelling  of  interaction  problems,  finite  element 
eddy  current  calculations  for  rigid  conducting  plates,  and  finite  element 
numerical  studies  of  coupled  magnetoelastic  problems.  Experimental  results 
by  others  and  analytical  solutions  for  limiting  cases  have  been  used  to  verify 
the  numerical  results  obtained  at  each  stage  of  the  investigation.  The 
Fortran  programs  developed  for  each  part  of  the  study  are  also  described. 

The  electrodynamics  and  continuum  mechanics  bases  of  the  problem  are 
presented  in  the  modelling  part  of  the  study.  v.The  Eulerian  nature  of  the 
problem  is  emphasized.  All  the  field  relation^  are  linearized  in  an  updated 
Lagrangian  sense.  A  stream  function  method  is  introduced  for  the  eddy 
current  calculations  on  thin  moving  plates.  The  couplings  between  the  eddy 
current  problem  and  the  vibration  of  the  plates  are  discussed  in  detail. 

Both  linear  and  geometrically  nonlinear  problems  are  formulated  for  the 
coupled  magnetomechanical  systems. 

Eddy  current  calculations  on  rigid  conducting  plates  form  the  basis  for 
later  consideration  of  the  coupled  magnetoelastic  problems.  Steady  state 
harmonic  eddy  current  analyses  are  performed  for  both  finite  and  infinitely 
long  conducting  plates.  Transient  time-integration  analysis  is  also  performed 
for  infinitely-long  plates.  In  all  cases,  the  nonlocal  effect  of  the  self¬ 
field  on  the  eddy  current  is  included.  Infrared  and  search  coil  experimental 
results  by  others  are  employed  to  verify  the  finite  element  numerical 
solutions. 

Linear  vibrations  of  infinitely  long  conducting  plates  excited  by  pulsed 
magnetic  fields  are  studied  via  finite  element  time  integration  analysis. 

Low  inductance  strain  gages  are  used  for  experimental  verifications  of  the 
numerical  results  obtained.  The  effects  of  pulse  duration  and  various 
kinds  of  interactions  on  the  coupled  systems  are  investigated.  A  finite 
element  formulation  for  the  induced  nonlinear  vibrations  of  infinitely 
long  conducting  plates  is  presented.  Numerical  difficulties  and  some 
limited  results  for  nonlinear  static  problems  are  discussed. 
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Chapter  1 


INTRODUCTION 


Materials  Interact  with  electromagnetic  fields  In  a  variety  of  ways. 
Some  electromagnetic  Interaction  phenomena,  such  as  photoelasticity  and 
piezoelectricity,  have  found  wide  application  In  experimental  mechanics 
and  In  Industry.  But  It  Is  only  recently  that  a  corresponding  theory  of 
sufficient  generality  has  been  developed  for  electromagnetic  Interaction 
phenomena  of  elastic  solids.  This  field  Is  generally  called  electro- 
magnetoelastlclty. 

Special  branches  of  the  electromagnetoelasticity,  such  as  piezo¬ 
electricity,  have  been  well  developed  and  apolled  widely  In  the  electronic 
Industry.  Cther  branches  of  the  general  theory  remain  unnoticed  and 
seldomly  Interest  structural  mechanicians.  However,  recent  developments 
In  the  high  energy  Industry  are  gradually  changing  this  situation. 

High  energy  devices  such  as  fusion  reactors,  superconductive  energy 
storage  devices,  MHD  generators,  and  magnetically  levitated  vehicles  are 
all  subjected  to  strong  electromagnetic  fields.  The  design  of  these 
devices  calls  for  an  understanding  of  the  nature  of  the  electromagnetic 
forces  and  the  Interactions  between  the  fields  and  the  structure.  Thus, 
although  still  relatively  unnoticed,  electromagnetoelasticity  Is  gradual¬ 
ly  emerging  as  a  new  field-structure  Interaction  theory  for  structural 
mechanicians. 

The  electromagnetoelastic  materials  of  Interest  to  structural 
mechanicians  are  conductlno  metals  and  superconducting  materials.  Both 
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types  are  used  as  major  structural  materials  In  the  high  enerqy  Industry 
and  sometimes  they  appear  In  composite  form.  The  electromagnetic 
phenomenon  In  such  materials  Is  basically  a  maqnetlc  one.  The  Inter¬ 
action  theory  between  the  field  and  the  elastic  body  Is  therefore  called 
maqnetoelasticlty. 

Scattered  studies  of  magnetoelastic  structural  problems  exist  In  the 
literature.  Knowledge  of  this  theory  Is  still  restricted  to  a  relatively 
small  group  of  specialists.  In  no  small  measure,  this  fact  Is  due  to  the 
complexity  of  the  subject,  requiring  knowledge  of  mechanics,  electromagnetic 
theory,  and  thermodynamics.  Even  for  relatively  simple  materials  like  non- 
ferrous  conductors,  the  Interactions  among  these  three  fields  are  by  no 
means  trivial. 

This  chapter  gives  a  brief  Introduction  to  the  ma<n  features  of 
magnetoelastic  structural  problems.  Section  1.1  describes  the  multi¬ 
disciplinary  nature  of  the  theory  of  maqnetoelasticlty  for  nonferrous 
conductors.  Section  1.2  Introduces  the  plate  as  a  special  structural  form 
In  the  general  theory  and  describes  the  modelling  approach  of  magneto¬ 
elastic  plate  problems  used  In  this  thesis.  Section  1.3  presents  the 
organization  and  an  overview  of  the  later  Chapters. 


1.1  The  Multidisciplinary  Nature  of  Maqnetoelasticlty 

Magnetoelasticity  for  nonferrous  conductors  encompasses  three  differ¬ 
ent  fields:  elasticity,  electrodynamics,  and  thermodynamics.  The  cou¬ 
pling  of  these  three  fields  Is  shown  schematically  In  Figure  1.1  and 
occurs  through  the  field  equations,  constitutive  eouatlons,  and  boundary 
conditions.  The  quantification  of  the  mutual  Interactions  and  their 
effects  on  the  behavior  of  the  elastic  structure  Is  the  purpose  of  this 
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Details  of  the  theory  of  magnetoelasticity  will  be  presented  In 
Chapter  Two.  However,  from  the  simple  figure  here,  one  notes  that  the 
electric  current  Interacts  with  the  magnetic  field  to  create  a  magnetic 
body  force  on  the  conducting  structure.  The  deformation  and  motion  of 
the  structure  then  perturb  and  change  the  electromagnetic  field.  The 
electric  current  generates  Joule  heating  which  diffuses  through  the 
conductor  body  and  creates  thermal  stresses.  The  thermoelastic  effects 
further  change  the  motion  of  the  conductor  structure  and  affect  the 
electromagnetic  fields.  During  the  whole  process,  the  distributions  of 
the  magnetic  body  force  and  the  electric  heating  depend  on  the  variation 
of  the  magnetic  field  and  all  other  field  variables.  The  resulting  body 
force  distribution  can  be  markedly  different  from  the  gravitational  body 
force  more  familiar  to  structural  mechanicians  and  may  cause  deformation 
gradients  Inside  the  structure  that  differ  from  the  conventional  kin¬ 
ematic  simplifying  assumptions  of  a  particular  structural  form.  Further¬ 
more,  coupled  problems  are  three-dimensional  as  the  magnetic  fields  can 
permeate  through  both  the  structure  and  the  free  space  surrounding  It. 

The  coupled  problem  is  also  a  dynamic  one  since  the  interactions  come 
through  the  motion  of  the  conductor  and  the  time  variation  of  the  elec¬ 
tric  and  thermal  fields  Inside  the  conductor  body.  The  solution  of  any 
but  the  simplest  structural  configurations  Is  difficult,  and  one  generally 
has  to  resort  to  numerical  means. 

l.£  Modelling  of  Magnetoelastic  Plate  Problems 

All  the  above  mentioned  features  of  magnetoelasticity  may  be  found 


for  the  problem  of  a  conducting  thin  plate  In  a  time-varying  magnetic 
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field.  But  here  the  thinness  of  the  plate  may  cause  the  electromagnetic 
and  mechanical  field  variables  to  have  certain  restricted  variations  across 
the  thickness.  Simplifying  assumptions  regarding  these  variations  may  then 
be  made  to  model  the  problem  In  a  form  more  amenable  to  analysis.  The 
study  of  the  Interactions  between  time-varying  magnetic  fields  and  the 
motion  and  deformation  of  thin  conductor  plates  is  the  topic  of  this  report. 
The  thermoelastic  effects  are  not  considered  In  this  work.  The  modelllno 
approach  Is  described  In  the  following  paragraphs. 

The  magnetoelastic  problem  of  a  thin  nonferrous  conductor  Is  basically 
an  electromagnetic  Induction  problem  coupled  with  a  problem  of  dynamic 
elasticity.  For  structural  problems,  one  considers  low  freguency  theories 
In  both  subsystems.  When  the  frequency  of  the  electromagnetic  field  Is 
sufficiently  low,  the  penetration  depth  of  the  magnetic  field  Into  the 
conductor  Is  large  compared  to  the  thickness  of  the  plate,  and  the  Induced 
current  (eddy  current)  Is  approximately  uniform  across  the  thickness  and 
flows  parallel  to  the  mid-surface.  One  can  therefore  assume  that 

(1)  the  normal  to  the  mid-surface  remains  unstretched  and  normal  to  the 
deflected  mid-surface  of  the  plate, 

(2)  the  eddy  current  density  Is  constant  across  the  thickness  of  the 
plate,  and,  from  this  assumption, 

(3)  the  normal  component  of  the  Induced  maonetlc  field  Is  constant 
across  the  thickness  of  the  plate. 

The  Induced  eddy  current  flow  Is  therefore  perpendicular  to  the  normal  and 
rotates  with  the  normal  when  the  plate  Is  vibrating. 

With  these  three  basic  assumptions,  and  a  few  others  Introduced  In 
Section  2.3.2,  the  plate  problem  Is  modelled  as  an  eddy  current  problem 
on  a  moving  conducting  sheet  coupled  to  a  vibration  problem  of  the  thin 


elastic  plate.  The  equations  governing  these  two  subproblems  will  all  be 
referred  to  the  mid-surface  of  the  plate  and  coupled  through  velocity  and 
the  calculated  magnetic  force. 

To  achieve  such  a  two-dimensional  model,  the  current  density  Is 
chosen  as  the  primary  variable  In  the  Induction  problem.  A  stream  function 
method  Is  Introduced  to  reduce  the  vector  eouatlon  Into  a  scalar  equation. 

The  Blot-Savart  law  Is  used  to  calculate  the  normal  component  of  the  In¬ 
duced  magnetic  field,  which  Is  basically  the  stream  function  for  the  In¬ 
duced  current  density.  The  singular  nature  of  the  Blot-Savart  law  Is 
avoided  at  the  mid-surface  of  the  plate  by  transforming  the  volume  Integra¬ 
tion  Into  a  surface  Integration.  The  expression  obtained  Is  then  substi¬ 
tuted  Into  Faraday's  law  to  obtain  the  desired  governing  equation. 

The  magnetic  force  calculated  will  have  components  normal  and 
tangential  to  the  plate.  The  equations  of  motion  of  the  plate  will 
contain  the  effects  of  both  components. 

1.3  Thesis  Overview 

The  theoretical  side  of  this  thesis  Is  treated  In  Chapter  Two.  To 
provide  a  proper  background  to  the  discussion  of  the  stream  function  method, 
the  theory  of  magnetoelasticity  for  nonferrous  conducting  materials  Is 
summarized,  together  with  the  linearization  procedure  of  the  various  electro¬ 
magnetic  and  mechanical  field  relations.  The  stream  function  method  Is  then 
discussed  In  detail.  Various  assumptions  used  are  examined  and  the  electro¬ 
magnetic  field  problem  reduced  to  a  single  Integro-dlfferentlal  eouatlon 
for  the  stream  function.  A  literature  survey  of  magnetoelastic  plate 
problems  Is  also  given. 
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Chapter  Three  Is  devoted  to  the  eddy  current  calculation  on  rigid, 
stationary  conducting  plates  since  this  subproblem  Is  not  trivial  and 
Is  In  Itself  a  significant  aspect  of  electromaonetlc  field  computation. 

The  governing  equations  for  eddy  currents  derived  In  Chapter  Two  are 
applied  to  finite  (two-dimensional)  and  Infinitely  long  (one-dimensional) 
plates  subjected  to  various  exciting  magnetic  fields.  The  numerical 
results  are  compared  to  some  experimental  data  obtained  by  using  Infra¬ 
red  sensing  techniques.  The  finite  element  method  Is  used  to  solve  the 
Integro-dlfferentlal  equation  for  the  stream  function. 

Chapter  Four  Is  concerned  with  the  linear  vibration  of  Infinitely 
long  magnetoelastic  plates.  Small  deformation  linearizations  are 
employed,  and  various  types  of  coupling  Investigated.  Forced  vibrations 
Induced  by  transient  magnetic  fields  generated  by  pulsed  currents  In 
nearby  colls  are  studied.  A  finite  element  stagqered  transient  analysis 
procedure  Is  used  for  the  coupled  set  of  equations.  Parameter  studies  are 
presented  to  explore  the  effects  of  the  different  time  characteristics  of 
the  magnetic  and  mechanical  subsystems.  Some  comparisons  of  the  calcula¬ 
tions  with  experimental  results  are  given. 

Chapter  Five  discusses  the  nonlinear  vibration  of  Infinitely  lonq 
magnetoelastic  plates.  The  formulation  and  algorithm  for  solution  are 
presented  based  on  the  linearized  equations  for  the  updated.  Incremental 
analysis  procedure.  A  two-dimensional  nonlinear  beam  element  Is  proposed 
for  the  geometrically  nonlinear  problems.  Some  limited  results  for  the 
uncoupled,  static  nonlinear  elastic  problems  are  presented,  and  difficul¬ 
ties  In  the  convergence  of  solutions  for  the  nonlinear  dynamic  cases  are 
discussed.  Some  suggestions  for  the  Improvement  of  the  numerical 
technique  are  also  given. 


Conclusions  on  the  results  of  the  present  study  on  magnetoelastic 
plates  are  drawn  In  the  last  Chapter,  together  with  some  sugaestlons  for 
further  research  In  this  area. 
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Chapter  2 


GENERAL  THEORY 


Several  features  of  the  theory  of  electromaanetic  interactions  in 
elastic  solids  make  It  unique  In  continuum  mechanics.  In  the  first  place, 
there  Is  no  generally  aoreed  upon  Interpretation  of  some  of  the  electric 
and  magnetic  vectors  In  continuous  bodies.  These  dlfferlno  Interpretations, 
together  with  the  fact  that  the  laws  of  classical  mechanics  are  Galilean- 
Invariant  while  the  Maxwell  eouatlons  are  Lorentz-invarlant,  olve  rise  to 
a  variety  of  formulations  of  electrodynamics  in  a  movlnn  medium.  In  the 
second  place,  the  separation  of  total  contact  force  Into  mechanical  and 
electromagnetic  parts  Is  not  unloue.  The  definitions  of  the  stress  tensor 
and  the  electromagnetic  body  force  thus  become  mutually  dependent. 

Finally,  for  materials  such  as  piezoelectric  solids,  the  Interaction 
enerqy  Is  linear  In  strain.  If  the  stress-strain  relations  are  to  be 
obtained  to  the  first  order  in  the  strains  by  differentiation  of  an  energy 
function,  then  that  enerqy  function  must  be  correct  to  the  second  order 
(Ref.  1).  Finite-strain  theory  with  all  the  complications  of  the  associated 
two-point  tensors  must  then  be  used  in  some  problems,  even  when  the  defor¬ 
mation  Is  small.  Such  features,  as  the  three  mentioned,  make  the  qeneral 
theory  much  more  complicated  than  a  simple  superposition  of  the  Maxwell 
equations  and  the  theory  of  linear  elasticity. 

Toupin  (Refs.  2,3)  was  the  first  one  to  use  finite-strain  theory  to 
clarify  some  problems  In  the  theory  of  Interaction  of  an  electrostatic 
field  with  a  perfectly  elastic  dielectric.  Tlersten  (Refs.  4,5)  then 
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developed  a  similar  theory  for  the  deformation  of  maanetlcally  saturated 
Insulators.  A  rioorous  treatment  of  the  maonetoel astlc  Interactions  has 
been  qlven  by  Brown  In  his  monoaraph  (Pef.  1).  The  treatment  Is  mainly 
for  saturated  non-dlsslpatlve,  and  non-conductive  materials  under  static 
magnetic  fields.  The  model  of  electrodynamic  theory  used  bv  Toupln  and 
Tlersten  Is  also  criticized  In  Brown's  work.  The  Interest  In  all  these 
studies  has  been  toward  an  understandlna  of  the  phenomena  of  electro- 
strlctlon  and  maonetostrlctlon  as  physical  problems. 

Electrodynamics  in  a  movlno  media  has  Iona  been  a  controversial 
subject.  Pao  (Ref.  6)  has  qiven  a  detailed  comparison  of  four  different 
theories  In  the  literature:  the  Minkowski  formulation,  the  Lorentz 
formulation,  the  statistical  formulation,  and  the  Chu  formulation.  Field- 
matter  Interaction  theories  for  stationary  and  movino  media  based  on  these 
different  models  are  also  presented.  Discussed  also  are  the  constitutive 
equations  and  boundary  conditions.  It  Is  shown  that  under  appropriate 
transformations  of  the  different  definitions  of  the  electric  and  maanetlc 
vectors,  the  four  models  can  be  made  equivalent  to  each  other.  However, 
the  expressions  for  the  electromagnetic  body  force,  body  couple,  and  energy 
supply  for  polarizable  magnetizable  moving  media  based  on  these  models 
differ  from  each  other.  For  nonferrous  metals,  however,  the  differences 
In  these  expressions  vanish  and  the  results  are  also  much  simpler. 

Field-matter  Interaction  theory  Is  a  special  branch  of  continuum 
mechanics.  The  treatment  of  this  subject  In  the  literature  basically 
follows  two  different  procedures.  The  stress  method  uses  the  stress  concept 
and  the  expression  for  the  electromagnetic  forces  on  a  material  volume. 

The  energy  method  uses  a  stored  enerqy  function  and  a  variational  principle. 
Parallels  exist  between  the  two  procedures,  and  in  theory  they  should  yield 
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equivalent  models.  However,  since  it  is  Impossible  to  measure  the  field 
variables  Inside  the  material,  some  assumptions  always  have  to  be  made 
in  either  method.  Several  self-consistent  theories  exist  In  the  litera¬ 
ture.  Some  classical  monographs  and  articles  on  this  subject  include  Brown 
(Ref.  1),  Penfield  and  Hans  (Ref.  7),  and  Truesdell  and  Toupin  (Ref.  8). 
Some  more  recent  ones  Include  Parkus  (Refs.  Q,  10,  11)  Hutter  and 
Van  de  Ven  (Ref.  1?),  Alblas  (Ref.  13),  and  Eringen  (Ref.  14).  In  each 
of  these,  extensive  references  to  other  papers  may  be  found. 

Beside  the  Interest  In  pure  theoretical  studies,  there  has  been 
research  work  on  this  subject  toward  Its  practical  applications  In  diverse 
fields  such  as  geophysics,  optics,  acoustics,  solid-state  devices,  and  so 
on.  Knopoff  (Ref.  15),  Wilson  (Ref.  If),  Kaliskl  (Ref.  17),  and  Dunkln 
and  Eringen  (Ref.  18)  studied  the  effects  of  maqnetlc  fields  on  elastic 
wave  propogatlon.  The  maanetoelasticity  theory  they  use  is  a  simple 
combination  of  Maxwell  equations  and  the  equations  of  linear  elasticity. 

The  Influence  of  the  electromagnetic  field  on  the  strain  occurs  solely  by 
means  of  Lorentz  forces  appearing  in  the  equations  of  motion.  The  motion 
of  the  material  affects  the  electromagnetic  field  only  through  the  velocity 
of  the  particles  in  Ohm's  law.  Electromagnetic  wave  propoaatlon  usually 
is  not  considered.  The  use  of  quasi-static  theory  further  simplifies  the 
field  problem.  Many  papers  of  this  type  have  been  reviewed  by  Paria 
(Ref.  19). 

On  engineering  aspects,  Montgomery  (Ref.  20)  and  Brechna  (Ref.  21) 
discussed  stress  analysis  in  maonet  design.  Woodson  and  Melcher  (Ref.  22) 
treated  electromechanical  problems  in  electric  machines  and  other  devices. 
Melcher' s  new  book  (Ref.  23)  directs  toward  electromagnetic  field-continuum 
interactions.  The  treatment  in  this  book  is  general,  but  the  problems 


4 


discussed  are  mainly  fluid  flow  ones.  Many  other  books  on  electric 
machines  and  energy  conversion  also  Include  electromechanical  problems. 
Quasi-static  theory  is  qenerally  used  for  analysis.  The  circuit  analoqy 
provides  a  convenient  tool  In  such  a  case.  Many  of  the  studies  are  for 
the  uncoupled  cases  only,  where  the  knowledge  of  the  magnetic  field 
determines  the  body  force  and  hence  the  motion  and  deformation  of  the 
solid.  Becker  and  Plllsbury  (Ref.  24)  Is  one  of  the  few  cases  where  the 
mutual  coupling  Is  taken  Into  account. 

Moon's  recent  article  (Ref.  25)  provides  a  broad  survey  of  possible 
engineering  applications  of  magneto-solid  mechanics.  Many  problems 
discussed  find  applications  In  high  energy  technolooy.  Desion  of  fusion 
reactors,  superconductive  energy  storage  devices,  MHD  generators,  etc., 
all  Involve  considerations  of  strong  magnetic  forces  imposed  on  the 
structural  components.  Moon  and  his  co-works  (Refs.  26,  27)  have  studied 
the  stability  and  vibrations  of  superconducting  maonet  systems  extensively. 
Mlya,  et  al  (Ref.  28)  studied  conductive  cylindrical  shells  In  pulsed 
axlsymmetrlcal  maqnetlc  fields.  A  new  monograph  edited  by  Moon  (Ref.  29) 
has  addressed  the  mechanics  of  superconducting  structures.  Papers  pub¬ 
lished  In  the  nuclear  enqlneerlnq  literature  on  the  structural  considera¬ 
tions  of  the  devices  are  also  numerous.  However,  apart  from  the  few 
studies  already  cited,  they  are  all  limited  to  the  uncoupled  cases.  The 
mutual  Interaction  between  the  field  and  the  response  of  the  structure 
has  received  little  attention. 

The  purpose  of  this  report  Is  to  study  the  mutual  Interactions  between 
the  time-varying  magnetic  field  and  the  Induced  motion  of  thin,  nonferrous 
conducting  plates  on  a  rigorous  continuum  mechanic  basis.  The  dynamics  of 
the  coupled  magneto-mechanical  system  Is  emphasized.  The  eddy  current 


problem  and  the  deformation-dependent  nature  of  the  magnetic  force  are 
given  special  attention.  Both  linear  and  nonlinear  deformations  of  the 
conduction  plate  are  considered. 

In  this  chapter  the  theory  of  maonetoelastlc  plates  and  the  formula¬ 
tion  approach  used  In  this  thesis  are  presented.  Section  2.1  summarizes 
the  theory  of  maqnetoelastlclty  for  nonferrous  conductors  based  on  the 
Lorentz  force  concept.  The  Eulerlan  description  Is  used.  Section  2.2 
discusses  the  linearization  procedure  of  the  field  equations.  The 
description  of  the  coupled  system  is  then  referred  to  a  known  conflqura- 
tlon  of  the  conductor  body.  The  approximation  of  the  procedure  Is  also 
discussed.  Section  2.3  addresses  thin  nonferrous  conduction  plates.  A 
stream  function  method  for  eddy  current  calculation  Is  Introduced.  The 
problem  Is  then  modelled  as  an  eddy  current  problem  on  a  movlna  conducting 
sheet  coupled  to  a  vlbratfon  problem  of  a  thin  plate  subjected  to  deforma¬ 
tion-dependent  load.  The  various  coupllnn  effects  and  the  mannetlc  body 
force  in  the  plate  are  discussed.  A  state-of-the-art  survey  of  magneto- 
elastlc  plate  problems  Is  also  alven. 

2.1  Maqnetoelastlclty  for  Nonferrous  Conducting  Materials 

The  material  treated  In  this  study  Is  assumed  to  be  homogeneous, 
isotropic,  and  possessing  finite  conductivity.  Polarization  and  magnetiza¬ 
tion  are  not  considered.  Furthermore,  no  thermal  effects  are  Included. 

The  theory  presented  below  Is  basically  a  theory  of  the  electrodynamics  of 
a  moving  media  coupled  to  the  theory  of  nonlinear  elasticity.  The  coupllna 
effects  come  from  the  Lorentz  force,  the  modification  of  the  electromagnetic 
field  relations  by  the  velocity  of  the  material,  and  the  boundary  conditions. 
The  presentation  follows  that  of  Hutter  and  Pao  (Ref. 30).  Details  of  some 
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derivations  are  omitted  when  the  reader  may  be  referred  to  this  and  other 
cited  papers. 

Consider  first  the  kinematics  of  a  material  body  In  three-dimensional 
space.  The  position  of  a  particle  In  the  reference  (material)  confiquratlon 
will  be  denoted  by  JL,  while  Its  position  In  the  present  (spatial)  confiqura¬ 
tlon  by  x.  Both  £  and  £  are  referred  to  a  cornnon  Cartesian  coordinate 
system  fixed  In  space,  as  shown  In  Flqure  2.1.  The  component  of  £  and  of 
all  other  vector  and  tensor  quantities  In  the  reference  conflouratlon  are 
indicated  by  Greek  subscripts  (a,  8  ■  1,2,3),  and  those  In  the  present 
conflouratlon  by  Latin  subscripts  (1,  j  *  1,2,3).  The  motion  and  deforma¬ 
tion  of  the  body  are  then  described  by 

x  «  x  (X,t)  or  x.  =  x.  (X  ,t),  1,  a  B  1,2,3  (2.1) 

%  11  'N*'1 

The  particle  velocity  and  acceleration  are 

v  =  x  (X»t)  -  -rr-  x  (X,t)  (2,2a) 


v  =  x  (X,t) 
^  dr  ^  ^ 


3 

3l 


V  +  V  V 
%  % 


V 


(2.2b) 


respectively.  Introduce  the  displacement  vector  u  with  components 


ui  =  xi  '  61a  Xa’  or  u~  =  6«-  x 


la  A1 


(2.3) 


The  particle  velocity  can  then  be  calculated  by 


v  =  u  +  v  v  *  u  (2.4) 

The  deformation  may  be  measured  by  either  the  Green's  strain  tensor  in 
the  coordinate  system  Xq 


!  b 

or  the  Almansi's  strain  tensor  in  the  coordinate  system 

ei j  =  **  ^u1,j  +  u.1,i  ‘  uk,1  ukj)  (?.5b) 

The  balance  laws  will  be  applied  to  quantities  measured  in  the  coordinates 

x^  for  the  present  conflouration  only. 

2.1.1  Electromagnetic  Field  Eouatlons,  Constitutive  Relations,  and 
Boundary  Conditions 

The  electromaqnetic  field  equations  for  movino  media  may  be  derived 


from  the  following  postulated  qlobal  laws  (Ref.  6). 

Gauss  -  Faraday  :  =  0  (2.6a) 

Faraday  :  K  h  *  d$  *  ‘  h  U  S  ’  (2.6b) 

Gauss  -  Coulomb  :  *  /y  Q  dV  (2.6c) 

Ampere  -  Maxwell  :  *c  «e  *  3  M  U  8  *  dS  +  's  *  dS  (?*6d) 

In  addition,  there  is  the  law  of  conservation  of  charoe 

is  4,  •  ‘  JF  •  dv  ’  0 


In  the  above,  £  Is  the  maonetic  Induction,  £  is  the  electric  displacement, 
the  £e  and  are  the  effective  electric  field  and  maonetic  field,  respec¬ 
tively,  and  is  the  effective  current,  q  is  the  free  charoe  density. 

All  integrations  are  taken  over  material  volume  V,  material  surface  S,  or 
material  circuit  C,  which  vary  with  time.  The  time  differentiations  should 
be  carried  out  by  applylno  the  transport  theorem 

/v  ♦  dV  -  /y  [||  +  v  •  (*♦)  1  dV  (2.7a) 

at  U  *  •  *1  ■  u  *  • d* 


(2.7b) 


+  V  x 


(2.7 c) 


1. 


!  7 


* 

In  which  £  =  fy£  +  x  v  •  £ 


<6  *  *)• 


For  nonferrous  metals  with  qood  (hut  finite)  conductivity,  the  free 
charqe  density  q  can  often  be  dropped.  Also,  since  the  frequencies  related 
to  vibrations  and  mechanical  waves  are  much  smaller  than  the  frequencies  of 
electromagnetic  waves  with  the  same  wave  lenoth,  one  can  regard  the  electro¬ 
magnetic  fields  as  quasistatic  when  investiqatlno  the  dynamic  response  of 
the  conductor  structures.  Mathematically  it  means  that  the  electric  dis- 
placement  5  =  0  and  J)  =  0.  The  field  equations  so  derived  will  be  of 
the  diffusion  type  rather  than  the  wave  type.  All  the  electromagnetic 
variables  inside  the  material  are  then  treated  as  convected  quantities. 

The  electromagnetic  field  equations  in  local  form  may  be  obtained  by 
substitution  of  Eqs.  (2.6)  into  Eos.  (2.7)  and  application  of  the  dlverqence 
and  Stokes  theorems  to  the  moving  volume  V  and  surface  S.  With  the  above 
simplifications,  Eqs.  (2.6)  become 


v  .  B  =  o 

(2.8a) 

it 

v  *  £e  3  -  «  =  -  It  l  -  v  *  x 

(2.8b) 

5  =  0 

(2.8c) 

*  *  lie  =  ie 

(2.8d) 

o 

II 

Q> 

> 

(2.8e) 

The  D,  and  in  Eqs.  (2.6)  are  fields  in  the  moving  body 

as  measured  by  an  observer  following  each  material  particle.  After  the 
application  of  the  transport  theorem,  the  £  and  £  in  Eos.  (2.8)  are  those 
In  the  moving  body  measured  by  an  observer  in  a  reference  frame  fixed  in 
space.  The  effective  fields  and  can  also  be  expressed  in  terms  of 
the  corresponding  fixed-frame  quantities  (J  and  J£,  respectively.  The 
expressions  are  different  in  various  formulations  for  general  materials. 


k 
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However,  for  nonferrous  metals  the  same  results  may  be  obtained  and  are 

t,  ■  l  *  *  "  S  «e  ■  «  -  *  «  *  l 

The  effective  current  density  measured  by  an  observer  moving  with  the 
charged  particle  Is 

4  *  n  - «  x  <?-9b> 

In  the  usual  case  of  quasistatic  field  problems,  no  free  charoe  can  exist 
inside  the  conductor  and  one  has  =  jJ.  1°  problems  that  Involve  conduc¬ 
tors  with  slits  or  cracks,  free  charge  may  be  accumulated  on  the  two  sides 
of  the  narrow  openino,  and  electric  arcing  may  occur  when  the  voltaoe 
across  the  opening  becomes  too  high.  In  such  problems,  the  effect  of  the 
free  charge  cannot  be  neolected.  In  the  present  study,  however,  free 
charge  accumulation  and  electric  arclno  are  not  considered,  and  thus  the 
relation  ^  ^  is  assumed. 

Note  also  that  since  £  =  c  j[,  the  quasistatic  approximation  effective¬ 
ly  Implies  that  the  permittivity  e  *  0  and  hence  £e  =  Jhl. 

The  electromaanetlc  constitutive  relations  In  this  study  are  limited 
to  the  classical  linear  relations  for  homogeneous,  isotropic  materials. 

For  slowly  moving  media  these  relations  are  invariant  with  respect  to  the 
transformation  of  reference  frames.  Namely, 

«  •  e  £.  1}  ■  »U.  <  '  C.'Oa) 

in  the  fixed  frame,  and 

«■*«,.  S’"*,,  le  - 

In  the  moving  frame.  In  the  above  y  Is  the  magnetic  permeability  and  a 
the  electric  conductivity  of  the  material.  In  view  of  Eqs.  (2.9),  one  has 
G  ■  *  <£  *  *  x  (J.lla) 

?■»(#-*««> 
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i  -  «  (£  +  *  *  V  <*."<:) 

Again,  because  of  the  quasi-static  approximation,  only  Fq.  (?.11c)  and 
Eq.  (2.11b)  with  =  0  will  be  needed.  For  nonferrous  metals  u  can  be 
replaced  by  u0  of  the  free  space  without  introduclna  any  sioniflcant  error. 

The  electromaanetic  boundary  conditions  may  be  derived  bv  applylno 
the  global  laws,  Eqs.  (2.6),  to  volumes  and  surfaces  which  contain  sinqular 
surfaces  and  lines  across  which  the  field  quantities  may  not  he  continuous 
or  may  suffer  a  jump.  The  transport  theorem,  Fos.  (2.7),  in  such  cases 
should  be  modified  to  account  for  the  discontinuities  of  the  field 
quantities.  The  modified  transport  theorem  is  then  applied  to  olobal  laws 
for  the  small  volume  or  "pillbox"  and  for  the  circuit  shown  in  Flo.  2.2. 

By  takinq  the  limits  of  vanishlno  volume  for  the  pillbox  and  vanishino 
enclosed  area  for  the  circuit,  one  may  derive  the  jump  conditions.  For 
nonferrous  metals  with  finite  conductivity  in  a  auasl-statlc  field,  they 
are 

5  •  £  f  3  *  0  (2.12a) 

i r  c  { ♦  x  *  %  ^ ' 0  (2-w 

pH  y  u  =  o  (2.12c) 

z  •  II  i  J  *  a  •  MS +  x  *  • 0 

In  which  jj  is  the  unit  vector  normal  to  the  boundary  surface  of  the  body 
in  the  present  conflouration,  £  £  1  =  ‘  denotes  the  jump  of  the 

quantity  A  from  the  positive  side  (positive  direction  of  ^)  to  the  neoative 
side  of  the  boundary  surface. 

The  modified  transport  theorem  and  the  details  of  the  derivations  may 
be  found  In  Hutter  (Ref.  31).  Dunkin  and  Frlnqen  (Ref.  IB)  used  a  different 
approach  of  derivation.  Fouations  (2.12c,d)  have  the  followlno  form  in 


Hutter  ( Pef .  31 ) : 

fl  x  H  II  -  X  x  C  J  x  4S  ’  qs  *  (?.13a) 

H  •  I  i  d  -  (*  •  *))  I $  J  + «  *  *  *  (jc  x  4S)  +  ^  ^  Qs> 

■-Its 

In  which  a  and  J  are  the  free  surface  charce  density  and  free  surface 
s 

current  density,  respectively.  For  aood  conductors  the  free  surface 
charoe  and  current  densities,  althouoh  achievable  in  hloh  freouency  cases, 
can  only  be  an  Idealization  that  occurs  in  conjunction  with  the  assumption 
of  a  perfect  conductor  (o  =  ”  )  (Pef.  3?). 

Equations  (2.8),  (2.9),  (2.11),  and  (2.12)  complete  the  description  of 
the  electromaanetlc  field  inside  movlna  nonferrous  metals  in  ouasl-statlc 
cases. 

2.1.2  Mechanical  Field  Fauations,  Maonetlc  Body  Force,  Constitutive 
Relations,  and  Boundary  Conditions 

The  mechanical  field  of  a  nonpolar  material  is  ooverned  by  the 
followlno  balance  laws  of  continuum  mechanics. 

Mass  :  /y  p  dV  *  0  (2.14a) 

Linear  Momentum  :  ^  /y-p  v  dV  =  ^  t^  dS  +  /y  p  f  dV  (2.14b) 

Anaular  Momentum:  ^  /y  £  *  o  v  dV 

*  £  *  t(n)  dS  ♦  /v  l  x  p  f  dV  (2.14c) 

:  3T  f\l  0  K  ‘  K  *  0  dV 


Fneroy 
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mm 


In  which  p  Is  the  mass  density,  t^  the  stress  vector  at  the  surface  with 
unit  vector  £,  f  the  body  force  per  unit  mass,  1)  the  strain  enerqy,  and  t 
the  energy  supply.  The  Inteorations  are  over  the  material  volume  V  and 
material  surface  S,  both  movlnq  with  the  material  contained  in  them. 

Thermodynamic  processes  are  not  considered  in  this  study.  The  heat 
flux  and  heat  energy  supply  terms  have  been  dropped  in  the  balance  law  of 


enerqy.  Introducina  the  Cauchy  stress  tensor  t  with 

iM  Tt  <”•  t,'"’  ■  "j  (MS) 

and  applvlnq  the  divergence  and  transport  theorems  to  Eos.  (2.14),  one 
obtains  the  mechanical  field  equations  in  local  form, 

^  p  +  p  v  •  v  =  0  (2.16a) 

p  3t  %  *  V  *  l  +  p  l  (2.16b) 

l-  l*  0  (2.16c) 

p  =  t  :  v  v  +  *  (2.16d) 


in  which  t*  is  the  transpose  of  the  stress  tensor  t,  and  t  :  v  v  =  t4.  v,  . 

<v  \  ^  \  1 J  1 , 1 

is  the  trace  of  the  tensor  product  of  t-  and  v.  .. 

*  j  j  * ' 

The  body  force  p£  and  energy  supply  +  arise  from  the  electromagnetic 
fields.  For  nonferrous  metals  they  have  the  followina  expressions: 

(2-15) 

Apnlyinq  the  electromaonetic  field  equations  (2.8)  and  some  vector  opera¬ 
tions  to  Eq.  (2.17),  one  obtains  another  form  for  the  body  force, 

p  t  s  *  *  l  or  P  f(  5  Ti;jJ  (2.17a) 

in  which  J  is  the  Maxwell  stress  tensor  for  the  maonetic  system  (Ref.  6), 


(2.19) 


T1j  ’  TJ1  *^Bi  Bj  '  Bk  Bk 

Substitution  of  Eq.  (2.17)  Into  Fo.  (2.18)  olves 

♦  =  *!e  ’  £e  (2.18a) 

Eouatlons  (2.1 6d )  and  (2.18)  represent  a  conversion  of  electromaonetlc 
energy  Into  mechanical  eneroy.  Substitution  of  Eo.  (2.17a)  into  Eq.  (2.14b) 
and  application  of  the  dlveroence  theorem  yields  the  balance  law  of  linear 
momentum: 

jfr  'v »  *  dv  ■  h  i  ‘  «  ds  *  h  I  •  *  «  t*.*") 

which  states  that  the  flow  of  electromaonetlc  momentum  Into  the  material 
body  is  completely  converted  into  a  kinetic  momentum  which  affects  the 
motion  of  the  body.  Detailed  discussions  of  the  electromaonetic  momentum 
flow  and  force  density  and  of  the  transfer  of  these  into  mechanical  momentum 
and  body  force  are  oiven  in  Pao  (Ref.  8),  Penfeld  and  Hans  (Ref.  7), 

Jackson  (Ref.  32),  and  Landau  and  Lifshltz  (Ref.  33).  Here,  it  is  only 
noted  that  the  force  actino  on  the  material  body  is  eoual  and  opooslte  to 
the  "force"  (transfer  of  the  electromaonetic  momentum)  actino  on  the  electro¬ 
maonetlc  field. 

No  useful  boundary  conditions  follow  from  the  global  laws  Eos.  (2.14a,c) 
By  applying  the  same  technique  as  used  in  Subsection  2.1.1,  one  may  derive 
two  more  boundary  conditions  from  Eos.  (2. 2D)  and  (2.14d). 

z  •  I*  +  ID  -  0  (2.21a) 

z  •  tt  i '  x +  zr  ■  V  x  •  Jr  l  J  =  0  (2.21b) 

uo  po 

in  which  the  Poynting  vector  in  Eq.  (2.21b) 
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represents  a  density  of  eneroy  flow  by  the  el ectromaonetic  field  across 
the  boundary  surfaces. 

The  mechanical  constitutive  equations  are  assumed  to  be  unaffected  by 
the  presence  of  the  electromagnetic  fields.  For  isotropic,  homooeneous, 
linear  elastic  materials,  they  are  the  usual  Hooke's  Law 

T1j  =  X  51j  ekk  +  eij  (?.?3) 

in  which  X  and  G  are  the  elastic  constants  independent  of  the  spatial 
coordinates  and  e^.  is  the  Almansi's  strain  tensor,  (Eo.(2.5b). 

2.1.3  Summary  of  Equations 

For  convenience  of  future  reference  the  pertinent  eauatlons  in  this 
section  are  cirouped  tooether  and  listed  below. 


Electromagnetic  Subsystem 
Field  Equations 

v  •  B  =  0  (2.24a) 

*  »  #  =  i 

tooether  with 

v  •  £  =  0  (2.24d) 

Constitutive  Equations 

<*>■*  (2.?5a,b) 

Jump  Conditions 

JC  •  IT  |U  =  0  (2.26a) 

JG  x  1 1  *  V  x  £  }  =  0  (2.26b) 

5  *  II  JtJ  D  =  0  (2.26c) 

c  •  U  n  =  «  •  (J[  +  *  *  ?)  U  *  0  (2.2«d) 


Mechanical  Subsystem 
Field  Eauations 


^  p  +  p  V  .  V  =  0 

(2.27a) 

p  at  *  E  7  *  i  +  p  Z 

(2.27b) 

p  HTU  =  T1j  vj.1  +  * 

(2.27c) 

where  = 

Constitutive  Equations 

T1j  3  x  51j  ekk  +  2G  eij 

(2.28a) 

T.  .  *  —  B.  B.  -  I—  6..  B.  B, 
ij  u0  i  J  2vq  k  k 

(2.28b) 

where  e(j.  -  %  <u(J  *  uJ>(  -  hk>)  ukJ) 

Jump  Conditions 

jc*I;i  +  II!  =  0 

(2.29a) 

(2.29b) 

2.2  Linearization  Procedures 

The  field  relations  and  boundary  conditions  presented  in  the  last  sub¬ 
section  are  formulated  in  Eulerian  coordinates  and  refer  to  the  as  yet  un¬ 
known  deformed  conf iouration.  To  render  the  equations  amenable  to  analysis, 
one  needs  a  known,  piven  reference  configuration  about  which  one  linearizes 
all  the  field  relations  and  boundary  conditions.  For  small  deformation 
cases,  the  undeformed  initial  configuration  of  the  conductor  is  a  natural 
choice.  For  geometrically  nonlinear  problems,  one  can  choose  for  reference 
the  last  known  configuration  of  the  conductor  during  the  solution  process. 


Such  a  procedure  Is  standard  for  nonlinear  mechanical  problems  and  is  called 
the  "updated  taqranglan  formulation"  in  the  literature  (Ref.  34). 

The  initial  undeformed  configuration  or  any  known  intermediate  con¬ 
tinuation  durino  the  solution  process  can  also  be  used  for  reference  for 
the  mechanical  subsystem.  To  have  a  consistent  linearization  for  the 
coupled  problem,  however,  the  same  procedure  should  be  applied  to  the 
electromagnetic  subsystem.  Since  a  purely  Lagranqlan  description  Is  not 
available  for  the  Maxwell  equations,  the  updated  Laaranolan  formulation 
is  a  more  natural  choice. 

Hutter  and  Pao  (Ref.  30)  have  developed  a  linearization  procedure  for 
movlnq  magnetizable  elastic  solids.  Alblas  (Ref.  13),  Van  de  Ven  (Ref.  35), 
and  some  Russian  authors  (Ref.  36)  also  suggested  similar  procedures.  In 
all  these  studies,  it  is  assumed  that  the  deformation  is  small  even  with 
very  large  electromagnetic  fields,  and  that  the  deformation  has  only  a 
minor  influence  on  the  fields.  The  dominant  effect  on  the  field  Is  the 
presence  of  the  undeformed  body.  The  linearization  of  the  Maxwell  equations 
is  then  carried  out  as  a  perturbation  on  the  undeformed  body.  All  the 
authors  considered  polarizable  magnetizable  materials.  The  constitutive 
equations  and  other  field  relations  must  therefore  also  be  linearized. 

For  the  geometrically  nonlinear  problems  in  this  thesis,  the  deforma¬ 
tion  Is  allowed  to  be  finite,  but  the  displacement  Increment  within  each 
incremental  solution  Is  assumed  to  be  small  enough  to  allow  the  field 
equations  to  be  linearized  in  a  similar  way.  The  approximations  made  In 
the  linearization  procedure  will  be  discussed. 

For  the  purpose  of  linearization,  three  configurations  of  the  body  are 
distinguished  as  shown  in  Figure  2.3, 


orespnt 

configuration 

t+At 


X 


2 


X1 


Figure  2.3  Initial,  reference,  and  present  configurations 
of  a  deformable  body  In  electromagnetic  field. 


(a)  the  Initial  undeformed  confiouratlon  with  coordinates  X, 

% 

(b)  the  Intermediate  reference  confiouratlon  with  coordinates  £, 

(c)  the  present  confiouratlon  with  coordinates 

All  X,  and  x  are  measured  in  the  same  Cartesian  coordinate  system.  The 
intermediate  state  can  either  be  the  initial  undeformed  state  in  the  small 
deformation  cases,  or  the  last  known  movino,  deformed  state  In  the  non¬ 
linear  cases.  All  the  field  vectors  and  tensors  In  this  state  are  assumed 
to  be  known  and  are  denoted  by  a  left  subscript  t. 

The  displacement  of  a  particle  from  its  intermediate  confi duration 


to  its  present  confiouratlon  Is  denoted  by 

U1  =  X1  -  Vi  V  =  6<~  x<  '  u~ 


la  A1 


(2.30) 


u^  Is  assumed  to  be  small  such  that  Ju.  J  <<1  and  |u.||  «  vQ,  where  vQ 
is  some  characteristic  wave  speed,  e.p.,  vQ  =  /  t/” p 

The  derivatives  in  the  present  confiauration  referred  to  the  inter¬ 
mediate  confi duration  then  become 


o _  3 _ a  _  /.  3 

3x .  "  3£  3Xj  '  V  la  "  3x . '  3C 
1  a  1  la 


(2.31) 


In  view  of  the  small  displacement  assumptions,  Eo.  (2.31)  may  be  replaced 


!_ .  .  3_ 

3X.  “v  51a  3£ 

1  a 


(2.32) 


However,  if  coupling  terms  that  are  linear  in  the  deformation  qradlent 

u  .  are  Important,  Eo.  (2.32)  must  be  replaced  by  the  more  accurate  relation 
a,  p 


3Xj  “  ^1a  ’  Ua,6  6l6^  3£ 
i  a 


(2.33) 


29 


2.2.1  Electromaqnetlc  Field  Eouations  and  Boundary  Conditions 

The  electromaqnetlc  field  variables  are  decomposed  Into  two  parts. 

The  first  part  corresponds  to  the  Intermediate  state  and  satisfies  Fos. 

(2.24)  -  (2.26).  The  second  part  represents  the  perturbation  which  accounts 
for  the  chanqes  of  the  externally  applied  field  and  the  effect  of  deforma¬ 
tion.  These  perturbations  are  denoted  by  lower  case  letters  and  are  assumed 
to  be  small  so  that  all  the  equations  may  be  linearized  with  respect  to  them. 
With  ^  +  fc*  £  *  +  £»  ancl  ^  *  t^  +  4»  one  obta1ns»  using  Fqs. 


(2.33)  and  (2.24), 


bo,a-  “„,S  <t'f  *  '  ° 


(S  +  ej 


u,  .  (*E  +  O  ,  + 


aPy  ' t  y  T  V  ,8  “  eaPy  U6,P  'tcy  T  V,«  T  W 


GaPy  by,  g  Eapy  U5,8  ^t^y  +  by\fi  V0  ^ 


^a,a  "  Ua,6  Mp  +  ^  8^  ,a  *  ° 


(2.34a) 

(2.34b) 

(2.34c) 

(2.34d) 


For  the  small  deformation  cases,  one  obtains,  using  Fq.  (2.32) 
instead  of  (2.33) , 


b  =  0 
a, a 


Eapy 


=  0 


GaPy  by,8  U0  ~  ° 


a,  a 


=  0 


(2.35a) 

(2.35b) 

(2.35c) 

(2.35d) 


which  are  the  same  as  Eqs.  (2.24). 

For  the  small  deformation  cases,  it  Is  assumed  that  the  boundary 
conditions  may  be  satisfied  on  the  surface 
S  (ca,  t)  *0 


(2.36) 


30 


with  unit  normal  vector  $ 


In  the  Intermediate  state.  The  boundary  conditions  are  then 


(2.38a) 

8  *  C  tl  +  a  +  $  *  tS  5  1 0 

(2.38b) 

JK  *  C  V»a  J  =  0 

(2.38c) 

£  •  IT  4  5  ■ 0 

(2.38d) 

When  the  second  order  theory  must  be  used,-  for  example ,  in  elastic  stability 

problems,  JJI  must  be  replaced  by  the  unit  vector  normal  to  the  deformed 

surface  S  (x. ,  t).  n.  is  connected  to  N  by  the  following  relation 
11  a 


n,  =  N  («.  +  N0  N  e0  -  u  .  )  (2.39) 

1  a  la  By  By  la  a,1 

This  derivation  Is  aiven  by  Parkus  (Pef.  11)  and  Hutter  and  Pao  (Ref.  30). 

The  reduction  of  Eqs.  (2.34)  to  (2.35)  in  small  deformation  cases  Is 

obvious  since  u  „  ^0  In  such  cases.  The  nonlinear  terms  In  the  perturbed 
a,B  ^ 

quantities  In  Eqs.  (2.34)  should  be  dropped  to  have  a  linearized  set  of 
equations.  The  linearized  equations  also  allow  for  the  inclusion  of  trans¬ 
port  current.  In  such  cases  and  ^  are  penerated  by  some  external 

electromotive  sources  connected  to  the  conductor  body  (Pef.  37). 


2.2.2  Mechanical  Field  Equations,  Magnetic  Body  Force,  and  Boundary 
Conditions 

The  mechanical  field  equations  (2.27)  may  be  linearized  by  the  same 
procedure  as  used  for  the  electromagnetic  equations.  For  the  small  de¬ 
formation  case,  one  has,  using  Eo.  (2.32), 
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*  -  n  ¥■ 


d  a 
3t  *  =  at  * 


ea6  =  ^  (u«,P  + 


(?.40) 

(2.41) 


The  eauations  of  motion  In  this  case  become  the  usual  Navier's  equation, 
2 

(2-42) 
at1' 


P  7?  “a  =  P  Ua,PB  +  (X  +  R)  Vfa  +  P  fa 


which  is  the  basis  of  the  maonetoelastlclty  theory  used  by  Parla  (Ref.  14) 
and  Knopoff  (Ref.  15).  The  boundary  conditions  In  this  case  refer  to  the 
surface  S  and  unit  normal  vector  Ji  in  the  orioinal  conflouration. 

For  geometrical ly  nonlinear  problems,  the  principle  of  virtual  work 
provides  a  more  convenient  formulation  method.  In  the  updated  Laorangian 
formulation,  the  principle  of  virtual  work  gives 

t+At-  ,  t+At  t...  t+At„ 

S„6  e  „  dV  =  P 


aP 


'«P 


(2.43) 


t+At 


where  R  Is  the  external  virtual  work  expression, 

t+Atp  _  t+At  t+At..  .  t+At  , t+At,  ;;  ,  t+At ...  ,9 

R  "  rt+AtA  *1  *U1  dA  +  ;t+Atv  p(  fi  ■  V  5ui  dV  (2-44) 

In  Eo.  (2.43),  t+AtS  are  the  Cartesian  components  of  the  2nd  Piola- 
Klrchhoff  stress  tensor  corresponding  to  the  conflouration  x^  at  time 
t+At  but  measured  i*-  the  confiauration  £  which  occurs  at  time  t, 

a 


t+4t<.  „  Jp _  -  t+A t_  , 

SaP  t+3TTea,i  ijCB,j 

P 

.t+At 


(2.45) 


and  s  e  .  are  the  virtual  variations  in  the  Cartesian  components  of  the 

Op 

Green- Lagrange  strain  tensor  in  the  configuration  x..  (at  time  t+At) 
referred  to  the  configuration  z,  (at  time  t). 


sUtt'„6  *  ■>  S  <Ua,6  *  UB,a  * 


(2.46) 


The  Integral  of  Piola-KIrchhoff  stresses  times  variations  in  the  Green- 


Laqranqe  strains  Is  defined  over  the  conflquration  r.  (at  time  t)  of  the 

a 

body. 

In  Eq.  (2.44)  is  a  virtual  variation  in  the  current  displacement 
Increment  components  u..  When  measured  In  the  configuration  f;  , 

i  a 


<5u  =  6U ,  6, 

a  11a 


(2.47) 


u.j  is  the  acceleration  in  the  configuration  x^.  t+Att.|,  t+Atf^,  and  t+Atp 
are  the  surface  traction,  the  body  force  per  unit  mass,  and  the  specific 
mass  in  the  configuration  at  time  t+At.  The  integrations  are  defined 
over  the  unknown  configuration  x^  of  the  body. 

When  referred  to  the  configuration  of  the  body  at  time  t,  the 
applied  forces  are  evaluated  usino  (Ref.  38) 


t+itt 

a  a,1  1 

(t+Atfa  -  “a>  tdV  ”  ^  ^+Atfi  -  V  «1a 
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t+At 


dV 


(2.48a) 

(2.48b) 


The  surface  traction  t+Att^  and  body  force  t+Atf^  are  generated  by  the 
electromagnetic  fields  and  depend  on  the  deformation  of  the  conductor. 

To  linearize  Eq.  (2.43),  the  Piola-KIrchhoff  stress  tensor  Is  first 
decomposed  Into  two  parts, 

t+At, 


~*a8  tTag  ^a8 


(2.40) 


where  tTQg  are  the  Cartesian  components  of  the  Cauchy  stress  tensor  that 
satisfy  Eq.  (2.27)  -  (2.29)  at  time  t,  and  S  .  are  the  Cartesian  components 

Ct" 

of  the  2nd  Piola-KIrchhoff  stress  Increment  tensor  referred  to  the  config¬ 
uration  i  at  time  t.  The  strain  Increment  t+Ate  ,  is  separated  Into 
a  a8 

linear  and  nonlinear  parts 


(2.50a) 
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t+At 


e«P  "  fn6 


Ca6  “  enB  +  naB 


(2. 50h) 


where 


'  11  *  V.' 


n  „  *  S  u  u  „ 
a6  Y  ,a  y.B 


(2.51) 

(2.52) 


The  constitutive  relation  between  stress  and  strain  Increments  used  now  Is 

S  -  C  .  c  .  (2.53) 

a?  aPyS  Y* 

Eq.  (2.43)  can  then  be  rewritten  as 

/,  C  Q  .  e  .  6e  „  tdV  +  /.  .t  0  Ap  *dV  ■  t+AtR  -  /„  .r  „  Be  e  tdV 
ty  apy6  y6  aP  ty  t  cip  aP  ty  t  ap  aP 


(2.54) 


which  Is  a  nonlinear  equation  In  the  Incremental  displacements  u  . 

a 

Linearization  of  Eq.  (2.54)  may  be  obtained  by  assumlna  that  e 


aB 


e  0»  „  *  Ae  „  and  that  S„*Cn.e..  The  equation  then  becomes 

cip  oP  <tp  ap  aPy<5  y& 


/*  C  .  .  e  .  Ae  B  tdV  +  /,  .t  „  A  n  .  tdV  •  t+AtR  -  f¥  .  Ae  .  *dV 
ty  aPy6  y6  ap  ty  t  ap  aP  ty  t  aP  aP 


(2.55) 


The  electromagnetic  body  force  may  be  calculated  from  the  linearized 
Maxwell  equations  (2.34)  or  (2.35).  Whichever  set  of  equations  Is  used, 
one  has 


pfa  "  X  t^a  +  (t>i  *  +  (4  *  t^a  +  (4  X  (2'56) 

The  last  term  In  Eq.  (2.56)  Is  of  the  second  order  In  the  perturbed 
quantities  and  should  be  dropped  to  have  a  consistent  linearization. 
However,  In  some  problems  this  last  term  Is  not  small.  This  point  will 
be  discussed  In  the  next  section. 
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The  mechanical  boundary  condition,  when  referred  to  the  boundary 
surface  S  and  the  unit  normal  vector  In  the  configuration  £  ,  Is 

V  (?'57) 

When  second  order  theory  needs  to  be  used,  £  should  be  replaced  by  ji  In 
Eq.  (2.39). 

The  elasticity  material  tensor  and  the  specific  mass  p  are 

deformation  dependent  In  the  theory  of  finite  elasticity.  However,  for 

small  strain  and  finite  deformation  problems  for  conducting  metals,  the 

deformation  effects  will  be  very  small,  l.e.,  K  .  -v  6  .  In  these  cases. 

a  f  i  ^  a  i 

In  this  study,  the  constitutive  tensor  Is  taken  the  same  as  In  the  un¬ 
deformed  original  configuration  of  the  material. 

2.2.3  Summary  of  Equations 

The  linearized  equations  of  the  theory  of  maqnetoelastlclty  are 
summarized  below. 

Electromagnetic  Subsystem 

B  =  .8  +  b,  E  *  .E  +  e,  J  +  j 

Field  Equations 

(1)  small  deformation  cases 


1  *  0 

(2.58a) 

a  ,a 

«8y  (tEy  +  V’B  +  IT  ba  "  0 

(2.58b) 

-  b  -  ii  .1  *0 

(2.58c) 

a  By  Yt“  o  -a 

■  0 

a, a 

(2.58d) 

(11)  larqe  deformation  cases 

b  -  u  *B.  » 

a, a  a,B  t  B,a 
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(2.59a) 
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Constitutive  Fauatlons 


.<1  +  j  =  o  J.E  +  e  +  c  „  u  .B  ] 
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Boundary  Conditions 

•  [[  £  }  =  0 

x  ^  tl + « +  i  *  t? 1  3  0 

e  *  tr  fe/u0  3  -  o 
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(2.59d) 


(2.60) 


(2.61a) 

(2.61b) 

(2.61c) 

(2.61d) 


where  Is  the  unit  vector  normal  to  the  boundary  surface 


S  (f.  t)  -  0, 


N  "  SV  t 
a  a 


’P  ’P 


In  the  second  order  theory,  jjl  Is  replaced  by  ji  olven  by 

n .  ■  N  (6 .  +  N  N  e.  6,  -  u  . ) 

1  a  la  6  y  Py  la  a,1 

Mechanical  Subsystem 

Equations  of  Motion 

(1)  small  deformation  cases 


(2.62) 


(2.63) 
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pfa  "  X  t$*  +  M  X  +  U  x  £  )a 


(2.65a) 


(11)  large  deformation  cases  (U.L.  formulation) 
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Constitutive  Equations 


S  „  =  C  .  .  e  . 
aB  aByB  yS 


where 
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Boundary  Conditions 


(2.67) 


N«  •  I  t*''t5„ .  *  T„«  J  ■  0 


(2.68) 


Is  replaced  by  n^,  Eq.  (2.63),  in  the  second  order  theory. 


2.3  Magnetoelastlc  Plates 

The  linearized  equations  of  magnetoelasticity  presented  in  the  previous 
section  are  specialized  for  thin  plate  problems  In  this  section.  In 
addition,  some  simpllfylnn  assumptions  are  introduced. 

2.3.1  Literature  Survey 

There  are  two  types  of  studies  of  maonetoelastlc  plates  In  the 


literature.  The  first  one  deals  with  the  Interaction  of  maqnetlc  fields 


with  soft  ferromaanetlc  plates.  The  maqnetlc  force  and  body  couple  are 
aenerated  from  the  maqnetlzatlon  of  the  material  In  this  case.  Eddy 
currents  are  usually  not  considered.  The  second  one  deals  with  conduct- 
Ino  materials.  In  this  case  the  eddy  currents  can  he  Generated  from  both 
the  time  variation  of  the  maonetlc  field  and  the  motion  of  the  plates. 

Most  studies  of  this  type  consider  the  effect  of  motion  only.  A  few 
studies  consider  transient  maonetlc  fields  but  assume  that  the  electric 
and  maonetlc  fields  are  uncoupled  and  that  the  eddy  current  can  have  no 
reaction  to  the  applied  field.  For  the  analytic  treatment  of  the  subject 
In  both  types,  almost  all  authors  consider  stronq  and  uniform  static 
external  maqnetlc  fields. 

Panovko  and  Gubanova  (Ref.  3°)  first  discussed  the  buckllno  of  a 
cantilever  ferromaonetlc  plate  In  a  uniform  static  transverse  maqnetlc 
field  using  a  negative  foundation  type  of  maonetlc  force  distribution. 

Moon  and  Pao  (Ref.  40)  recognized  that  in  a  uniform  maonetic  field  there 
can  be  no  net  force  on  a  maonetlzed  body.  A  model  based  on  a  distribution 
of  maonetlc  couple  alonq  the  length  of  the  plate  was  built,  with  the  jump 
conditions  of  the  maonetic  field  calculated  on  the  surfaces  of  the  deform¬ 
ed  plate.  The  calculated  critical  maonetlc  fields  aaree  well  Qualitative¬ 
ly  with  the  experimental  data,  but  the  experimental  data  are  lower  than 
the  theoretical  values  by  a  factor  of  about  1.8.  Since  then,  a  number  of 
theoretical  and  experimental  studies  have  been  conducted  trylno  to  resolve 
this  discrepancy,  Includina  the  development  of  a  self-consistent  theory 
of  linear  ferroelastlc  contlnua  by  Pao  and  Yeh  (Ref.  41).  A  detailed 
review  of  many  of  these  papers  may  be  found  in  Moon  (Ref.  25). 

More  aeneral  boundary  value  problems  have  been  stualed  by  Van  de  Ven 
(Ref.  35)  and  Parkus  (Ref.  11).  In  Parkus  (Ref.  11),  rectangular  and 


circular  plates  of  an  elastic,  soft  ferromaonptlc  material  In  a  uniform, 
static  transverse  maqnetic  field  are  discussed.  The  Unearned  magneto- 
elastlc  theory  based  on  Mutter  and  Pao's  procedure  (Ref,  30)  is  used  for 
analysis.  The  equation  of  motion  of  the  plate  is  derived  by  assumlnq 
uniform  tangential  induced  electric  fiplds  and  uniform  normal  induced 
maqnetic  fields  across  the  thickness  of  the  plates.  The  classic  Klrchhoff 
assumption  is  used  for  the  kinematics  of  the  plate.  Eddy  current  Is  not 
considered,  and  the  effect  of  the  in-plane  deformation  on  the  flexural 
deformation  is  also  nealected.  The  influence  of  the  mannetic  field  on  the 
deformation  of  the  plate  is  modelled  as  a  surface  effect.  The  critical 
maqnetic  field  that  causes  the  buckling  of  a  clamped  circular  plate  has 
been  calculated.  The  resultina  expression  contains  the  3/?  power  of  the 
thickness  to  radius  ratio  for  an  axisvmnetrlc  bucklinq  mode,  the  same  as 
in  the  Moon  and  Pao  model  for  a  periodically  pinned  bpam-nlate  in  sinusoi¬ 
dal  deformation.  No  comparison  to  experimental  data  of  this  result  is 
reported  In  this  work.  However,  it  shows  unusual  completeness  of  deriva¬ 
tion,  startinq  from  the  dlpole-Amperlan  current  model  of  electrodynamics 
In  a  movlnq  body,  through  the  linearization  of  the  various  field  enuations 
and  boundary  conditions,  to  the  modellinq  of  the  specific  boundary  value 
problem  of  magnetoelastic  plates. 

Kaliskl  (Ref.  42)  and  Dunkin  and  Erinaen  (Ref. 18)  studied  the 
vibrations  and  wave  propaaations  in  ferromaqnetic  plates  under  the  action 
of  uniform,  static  magnetic  fields.  The  theories  Included  conductinq 
currents,  but  did  not  predict  any  Instability  phenomena.  Moon  and  Pao 
(Ref.  43)  in  another  paper  showed  that  the  freauency  of  vibration  of  the 
beam  plate  decreases  as  the  maanetlc  field  increases.  The  vibration  and 
stability  problems  are  then  connected  to  each  other,  and  the  critical 


magnetic  field  can  be  determined  from  the  relation  between  the  frequency 
and  the  field.  In  two-dimensional  ferromaanetlc  plates,  Moon  (Ref.  44) 
and  Srlnivasan  (Ref.  45)  found  that  the  maqnetlc  field  can  either  decrease 
or  Increase  the  natural  frequencies,  dependinq  on  Its  direction.  To 
achieve  bucklinq  and  a  decrease  In  frequency,  the  field  has  to  be  nearly 
parallel  to  the  plate  normal.  In  Moon  (Ref.  44),  the  bucklinq  theory  of 
beam-plates  has  also  been  extended  to  two-dimensional  ferroelastic  plates 
in  oblique  mannetic  fields. 

For  nonferrous  conductinq  materials,  Dunkin  and  Frinqen  (Ref.  18) 
discussed  the  vibration  of  an  infinite  beam-plate  in  a  static,  tanqentlal 
maqnetlc  field  and  showed  that  the  eddy  current  acts  as  a  damplno  force. 
Their  treatment  of  the  interaction  terms,  however,  is  not  complete. 
Ambartsumyan,  et  al.  (Ref.  35)  have  published  a  monopraph  (In  Russian)  on 
maqnetoelastlclty  of  thin  shells  and  plates  under  the  action  of  static 
maqnetlc  field.  In  a  recent  paper,  Ambartsumyan  (Ref.  46)  discussed  the 
oscillations  of  thin  conductinq  plates  In  static,  uniform  transverse 
maqnetlc  field.  The  linearized  maonetoelastlc  equations  (Ref.  36)  were 
used  for  analysis.  The  induced  tangential  electric  field  and  the  Induced 
normal  magnetic  field  are  assumed  uniform  across  the  thickness  of  the 
plate.  An  improved  bendlnq  theory ^of  plates  (Ref.  47)  is  used,  and  the 
shear  deformations  of  the  plate  considered.  The  electromaanetlc  boundary 
conditions  are  Imposed  on  the  surfaces  of  the  undeformed  plate.  The 
calculated  electromagnetic  body  force  is  in-plane  and  includes  the  effects 
of  both  the  rotation  of  the  plate  normal  and  the  shear  deformations.  The 
equations  of  motion  are  then  Integrated  and  averaged  over  the  plate's 
thickness.  The  equations  derived  are  applied  to  Infinitely  lonq  plates, 
simply  swpoorted  on  the  edges.  The  comparison  of  the  results  to  the 
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predictions  based  on  the  assumption  of  undeformablo  normals  shows  very 
different  behavior.  The  latter  Indicates  the  classical  phenomena  of 
decreasing  frequency  and  buckllno  with  Increasinn  field  in  all  modes. 

The  former  Indicates  dampina  without  oscillations  In  the  first  mode  for 
all  values  of  the  transverse  shear  effect.  For  the  hiqher  modes,  the 
deformable  normal  theory  shows  decreasing  frequencies  and  Increasinn  damp- 
inn  for  small  thickness-to-lenoth  ratios  at  the  benlnnlnn,  and  increasina 
frequencies  at  larner  thickness-to-lenath  ratios.  No  experimental  data 
are  reported  to  verify  these  results.  However,  it  clearly  Indicates  the 
importance  of  maonetic  shear  forces  In  such  problems. 

In  an  unpublished  report.  Moon  (Ref.  37)  presented  some  preliminary 
studies  of  the  vibrations  of  conduction  plates  carryinq  stronq  transport 
currents.  Fdelen's  nonlocal  variational  calculus  (Ref.  48)  was  applied 
to  the  Lagranqlan  function  of  the  coupled  electromechanical  system. 
Electrical  heatinq  was  not  considered.  The  results  indicate  stiffening 
of  the  system;  however,  further  studies  are  needed  for  confirmation. 

Moon  (Ref.  25)  has  conducted  experiments  on  circular  aluminum  plates 
levitated  above  a  nelqhborinn  coll  with  sinusoidal  currents.  Clrculatlna 
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currents  of  up  to  10  A  were  Induced  In  the  plate  tooether  with  the 
associated  Joule  heating.  To  minimize  the  thermal  effects,  resonance  and 
free  vibration  techniques  were  used  within  10  seconds  after  the  current 
was  turned  on.  The  chanaes  in  frequency  sauared  were  found  to  be  linearly 
related  to  the  souares  of  the  exciting  coll  currents.  The  frequencies  of 
the  diameter  modes  were  found  decreased  and  that  of  the  symmetrical  mode 
Increased,  However,  thermoplastic  effects  cannot  be  separated  from  the 
observed  channes  In  frequency.  Further  studies  have  been  abandoned  be¬ 
cause  of  the  inevitable  complications  of  Joule  heatinq.  Clearly,  all  the 
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maqnetoelastic  studies  of  conductinq  bodies  have  to  be  subjected  to  the 
same  limitation  (Ref.  25). 

2.3.2  Stream  Function  Method  for  Eddy  Current  Problem  on  a  Movino  Plate 
In  the  last  subsection  some  research  on  the  vibrations  of  thin 
conductinq  plates  under  the  action  of  a  uniform,  static  field  have  been 
reviewed.  The  discussions  there  serve  to  reveal  some  aspects  of  the 
physical  nature  of  the  problem.  In  this  subsection,  forced  vibrations 
of  plates  Induced  by  time-dependent,  nonuniform  maonetic  fields  will  he 
considered.  Eddy  currents  can  be  qenerated  by  both  the  motion  of  the 
plate  and  the  time  variation  of  the  maqnetlc  field.  Practical  non- 
uniform  fields  such  as  those  aenerated  by  current-carryino  coils  are  of 
special  Interest.  The  tanoential  component  of  such  a  field  will  exert 
a  normal  force  on  the  plate  and  cause  a  forced  vibration.  The  induced 
vibration  and  the  mutual  Interactions  between  the  fields  and  the  motion 
are  the  object  of  the  study. 

The  linearized  electromaonetic  field  relations,  Eos.  (2.58),  (2.60), 
and  (2.61)  will  be  used  for  analysis.  They  are  referred  to  the  confiaura- 
tlon  and  boundary  surfaces  of  the  undeformed  plate  in  the  small  deformation 
cases,  and  to  the  confiouration  and  boundary  surfaces  of  the  last  known 
position  of  the  plate  in  larae  deformation  cases.  The  unit  normal  vector 
to  the  boundary  surface  is  defined  by  Fo.  (2.62).  In  makino  such  a  choice, 
one  loses  the  couplina  terms  in  Eos.  (2.59).  An  iterative  procedure  will 
be  suqqested  for  the  simultaneous  satisfaction  of  both  the  plectromaonetlc 
and  mechanical  field  eauations.  This  procedure  is  basically  an  extension 
of  the  Iteration  scheme  for  the  out-of-balance  forces  widely  used  in 
geometrically  nonlinear  elastic  analysis. 
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Eos.  (2.64)  and  (2.66)  will  be  applied  to  the  mechanical  subsystem 
In  the  small  deformation  and  In  the  Incremental  laroe  deformation  analysis, 
respectively.  The  classical  Klrchhoff  assumption  of  the  undeformable 
normal  to  the  midsurface  Is  adapted.  Transverse  and  In-plane  shear  defor¬ 
mations  are  not  considered. 

The  followlnq  two  basic  assumptions  are  made  In  this  study: 

(1)  the  Induced  current  density  ^  Is  constant  across  the  thickness  of 
the  plate  and  is  parallel  to  the  midsurface  of  the  plate; 

(2)  the  normal  component  of  the  Induced  maanetlc  field  Is  constant 
across  the  thickness  of  the  plate. 

These  assumptions  differ  slightly  from  what  other  investigators  have  used 
In  that  one  neglects  the  normal  component  of  the  current  inside  the  plate. 

The  first  assumption  is  a  pood  approximation  when  the  penetration 
depth  of  the  magnetic  field  Is  large  compared  to  the  thickness  of  the 
plate,  since  the  boundary  condition  Eq.  (2.61d)  will  force  the  current 
to  flow  parallel  to  the  surface.  The  second  assumption  is  a  valid  one 
for  such  a  current  distribution,  since  the  induced  current  is  related  to 
the  normal  component  of  the  Induced  magnetic  field  through  Eg.  (2.58c). 

The  penetration  depth  (skin  depth)  6  is  large  when  the  freouency  of  the 
magnetic  field  is  low.  It  Is  expressed  by  the  foil  owl  no  formula  (Ref.  49) 

62  *=  JL  (2.6«) 

voes 

In  which  a)  Is  the  frequency  of  the  magnetic  field  and  <j  the  conductivity 
of  the  material . 

Except  for  these  two  assumptions,  the  applied  external  magnetic 
field  Is  assumed  unaffected  by  the  generation  of  the  Induced  eddy  current. 
That  Is,  the  power  supply  cf  the  external  driving  coll  Is  assumed  to  be 
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current  controlled  Instead  of  voltaqe  controlled.  The  generation  of  the 
eddy  currents  In  the  external  coll  thus  does  not  need  to  be  considered. 

With  the  above  assumption,  the  spatial  and  time  variations  of  the 
applied  field  can  be  calculated  directly  from  the  drlvlna  current  and 
geometry  of  the  external  coll.  Only  the  Induced  parts  of  the  electro¬ 
magnetic  field  variables  need  to  be  determined. 

To  solve  for  the  Induced  field  variables.  If  Eqs.  (2.58)  In  their 
present  form  are  used,  the  electromagnetic  field  eouatlons  In  reolons 
outside  the  plate  need  to  be  considered.  Poundarv  conditions  for  the 
Induced  field  at  Infinity  also  have  to  be  Imposed.  To  avoid  solving  such 
a  three-dimensional  problem,  the  current  density  Is  chosen  as  the  primary 
variable  In  this  study.  The  electromagnetic  field  problem  will  be  re¬ 
formulated  by  uslna  a  generalization  of  the  stream  function  method 
suggested  by  Moon  (Ref.  25).  Since  Eqs.  (2.58)  and  (2.60)  are  the  same 
as  Eqs.  (2.24)  and  (2.25),  respectively.  In  the  followlna  derivation  the 
notation  employed  Is  that  of  the  total  £,  £,  and  ^  In  the  present  con¬ 
figuration,  but  referred  to  the  reference  configuration. 

Since  the  current  density  ^  is  uniform  across  the  thickness,  one 
defines 

\  t2-70) 

where  h  Is  the  thickness  of  the  plate.  From  the  continuity  condition 
Eq.  (2.24d) : 

v*^  =  hv*^*0  (2.71 ) 

a  stream  function  ^  is  Introduced  for  the  current 

|  *  v  *  (ji  |l)  »  •  |l  x  7  ^  (2.72) 

where  Is  given  by  Fq.  (2.62).  Comparing  Eq.  (2.71)  to  (2.24c),  one 


finds  that  if'  Is  basically  the  normal  component  of  the  Induced  magnetic 


field  B‘, 


.  .  h  r1 

* '  %  N 


(2.73) 


The  total  magnetic  field  Is  separated  Into  the  applied  external  part 


|°  and  the  Induced  part  J^, 

« *  f  *  e1 


(2.74) 


The  Induced  part  Is  related  to  the  eddy  current  distribution,  under  the 


quasi-static  assumption,  through  the  Biot-Savart  law 
'  *  4’h  J* 


(2.75) 


Eg.  (2.24a)  Is  automatically  satisfied  here  where  £  is  at  the  mid-surface 
of  the  plate. 

The  stream  function  may  now  be  determined  by  substitution  of  Egs. 
(2.72),  (2.74)  and  (2.75)  Into  the  Faraday's  law,  Eo.  (2.24b), 


V  «  V  x  (i|i  )J)  *  .  ah  (J^°  +  ^)  +  nh  ?  *  (J  *  £°) 


(2.76) 


Only  the  normal  component  of  the  vector  eguatlon  will  be  used  because  the 
In-plane  current  Is  of  primary  Interest.  Note  also  that  the  last  term  In 
Eq.  (2.76)  has  been  linearized  by  dropplno  the  self-induced  field.  The 
variation  of  *  £°  across  the  thickness  contains  all  the  shear  deformation 
effects.  However,  since  the  flexural  motion  of  the  plate  Is  of  primary 
Interest  In  this  research,  this  term  will  refer  to  y  on  the  mid-surface 
of  the  plate  only.  In  small  deformation  cases,  It  Is  the  oriolnal  mid-plane 
of  the  plate.  In  the  large  deformation  cases.  It  Is  the  curved  mid-surface 
in  the  reference  configuration  of  the  deformed  plate. 


To  have  a  oeneral  description,  let  an  orthoaonal  curvilinear  coordinate 
system  r  be  Introduced  on  the  mid-surface,  Fla.  2.4,  with  unit  base  vectors 
£  (n=l  ,2,3) .  jg  Is  *he  same  unit  normal  vector  JJI  defined  by  Fa.  ( ? . 62 ) . 

One  has,  on  the  mid-surface, 

T  i  =  r  i  (c-j»  C  2 ) ,  F  p  =  f  2  ( C  *  F 2 )  ,  F  ^  =  F  3  ( F  i  *  ^  (2.77) 


The  Infinitesimal  line  element  on  the  mid-surface  Is  qlven  by 
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The  current  density  on  the  surface  can  be  written  in  the  form 

l  =  !1  $1  +  r2  ('l*  c2>  €? 

Then, 

’  •  1  ■  %  <a2?  V  i  ■» 


i  =  _L  *L_  ,  t _ L  £L_ 

1  ^22  ^^2  2  a-ji  3 c i 


(2.7P) 

(2.79 a) 
(2.79b) 


(2.81) 


(2.82) 


(2.83) 


The  normal  component  of  the  left  hand  side  of  Fq.  (2.76)  becomes 

e,  •  [V  x  V  x  (*  e  )]  - - 1 -  [|—  (!^.  |1-)  +|~  (111  It-)]  (2.84) 

^3  ^3  a^  ^22  ^^1  ^11  ^^1  9^2  ^22  ^^2 


Let  G  =  K  -  JC'»  R  =  I?!  s  k  -  JC‘1-  Note  that  ?/r3  b  *1 

where  is  the  gradient  and  operates  on  the  £'  variables.  One  has,  by 


Figure  2.4  Orthogonal  curvilinear  coordinate  system^ 
on  the  mid-surface  of  the  deformed  plate; 

C  is  the  Cartesian  coordinate  of  the  reference 

*va 

configuration  of  the  deformed  plate. 


substitution  Fa.  (? .1?)  Into  (?.75), 
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The  normal  component  of  which  Is 

6il  <S>  ■  IrarV  I?i  »  '  'i  ■  S> dv' 
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Note  that  because 
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one  has 
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G  =  (VT  ♦  *  !&)  $3  +  vi  ($3  *  V 
By  uslnq  the  followina  relations, 

/v  7  •  j  dV  =  /s  |  •  dS 
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(2.88) 


(2.8«) 


(2.90a) 


(2.90b) 
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one  has 
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The  second  surface  Inteoral  becomes 


A  ;S'  l  V1  ♦  *  S  dS< 


(?.O0c) 


(?.91) 


under  the  assumption  of  undeformahle  normals,  (j^  •  N)  •  =  0. 


The  volume  inteoral  contains  the  1/R  sinoularitv  and  diveroes  at  the  mid¬ 
surface  as  R  approaches  0,  However,  for  flate  plates,  •  N  =  1  and 
V1  ♦  *  %  ®  0  for  all  £* .  To.  (?.°1)  becomes 


BN  W  =  TT  *  W  +  57F  ;V  (S?  *  V  *  V1  *  $3  dS’ 


which  is  free  from  the  singular  nature  of  the  Riot-Savart  law  at  the  mid¬ 
plane  of  the  plate. 

It  Is  observed  that  the  v.  •  Q  term  in  Eq.  (?.Q1)  is  related 

I  % 

to  the  curvatures  and  their  derivatives  of  the  deformed  plate  and  is 
therefore  of  the  second  order  in  maonitude.  Furthermore,  thp  tvne  of 
sinoularitv  associated  with  it  has  been  reduced  by  the  order  of  one. 

One  must  assume  finite  deflections  but  small  rotations  in  nonlinear  prob¬ 
lems  to  neglect  the  last  term.  Eo.  (2.Q1)  then  takes  the  followino  form 
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The  time  rate  of  chanoe  of  which  Is 
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(2.94) 


The  first  three  terms  obviously  correspond  to  the  riold  body  state.  The 
other  terms  represent  the  perturbation  caused  by  the  relative  rotation, 
transverse  deflection,  and  stretchina  of  the  plate  elements.  The  riold 
body  terms  are  expected  to  have  the  dominant  effect.  The  coupllnq  terms 
will  be  discussed  In  the  next  subsection. 

For  riqld  flat  plates,  Fla.  2.5a,  one  has  (Ref. 25) 
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In  the  derivation  of  the  above  equation,  the  boundary  condition  on  the 
narrow  edqes  of  the  plate  Is  ♦  *  o.  For  a  rlqid  plate  Infinitely  lonq 
In  the  y-dlrectlon  and  of  finite  width  l  in  the  x-directlon,  Flq.  2.5b, 
with  fields  Independent  of  the  y-coordinate,  one  can  Intearate  the  second 
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term  In  Eq.  (2.95)  with  respect  to  v'  and  obtain 
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Eqs.  (2.95)  and  (2.96)  will  be  used  In  Chapter  Three  for  eddy  current 
calculations  on  Hold  plates. 

The  normal  component  of  the  last  term  In  Fo.  (2.76)  Is 
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where 

A]  =  u?  R3°  -  u 3  R?°,  A?  =  u3  ^  R3°  (2. OR) 

In  the  r  -  coordinate  system. 
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The  three-dimensional  transient  field  problem  thus  reduces  to  a  two- 
dimensional  one  on  the  mid-surface  of  the  plate.  To  the  above  equations, 
one  needs  to  add  the  boundary  conditions  for  ii'  described  in  the  followlna. 


Assume  a  simply  connected  plate.  At  the  boundary  curve  C,  Flo.  2.5a, 
the  current  density  I  should  be  tanoent: 


l  ■  *n  '  ° 

(2.QQ) 

or 

<v  *  x  V  *  $n  =  v  *  *  $t  c  0 

(2.100) 

By  Intearatlno  the  above  equation  alona  the  boundary, 

one  can  replace 

Eq.  (2.99)  by 

<i=0  on  C 

(2.101) 

2.3.3  Coupllnq  Effects 


The  kinematic  nature  of  the  Rlot-Savart  law  has  been  recoonized  In 
the  last  subsection.  Eq.  (2.94)  has  been  derived  based  on  the  assumption 


i  v 


of  nondeformable  normal  and  small  rotation  Qf  the  plate  element.  The 
divergence  property  of  the  Biot-Savart  law  on  the  mid-surface  has  been 
avoided  by  transforming  the  volume  intenral  Into  a  surface  intenral.  The 
normal  component  of  the  self-field  on  the  mid-surface  of  the  plate  has 
therefore  been  modelled  as  caused  by  a  distribution  of  sinoularities  on 
the  top  and  the  bottom  surfaces  of  the  plate. 

The  coupllna  terms  in  Fa.  (2.94),  however,  are  too  complicated.  The 
effects  of  the  motion  on  the  field  come  throuqh  the  relative  rotational 
and  translational  velocities  of  the  whole  plate  in  a  aloha!  way.  However, 
these  terms  correspond  to  the  interaction  between  the  motion  of  the  plate 
and  the  self  maqnetlc  field  and  are  of  the  second  order  In  the  theory  of 
the  coupled  problem.  One  therefore  drops  these  coupllna  terms  and  considers 
only  the  local  coupllna  effect,  i.e.,  the  last  term  in  Eq.  (2.7B),  which 
has  already  been  linearized.  The  coupllna  terms  in  Fq.  (2.94)  mlqht  be 
important  for  transport  current  problems,  however  (Ref.  37). 

With  the  above  linearization,  the  eddy  current  problem  is  formulated 
as  follows.  For  finite  deformation  incremental  analysis,  one  has 
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In  the  small  deformation  cases.  It  reduces  to 
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where  u,  v,  and  w  are  the  In-plane  x-dlrection  velocity,  the  in-plane 
y-direction  velocity,  and  the  transverse  velocity,  respectively.  In  the 
small  deformation  of  infinitely  Iona  plates,  it  is 
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The  effect  of  in-plane  deformation  has  been  included  in  Fos.  (2.10?)  to 
(2.104).  The  correspondina  equations  of  motion  are  provided  by  Fqs. 

(2.66)  to  (2.6 7)  for  the  finite  deformation  incremental  analysis.  For 
small  deformation  cases,  they  can  be  obtained  bv  apolvina  the  nondeformable 
normal  assumption  to  Fq.  (2.64).  In  the  present  study,  they  are  taken  to 
be  the  classical  equations  for  plates  sub.iected  to  both  in-plane  and  out- 
of-plane  loads,  and  are 
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For  small  deformations  of  Infinitely  lone  plates,  Fos.  (2.105)  reduce  to 
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where 


(2.107) 


r  =  2i  i  ,  _ 
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2.3.4.  Magnetic  Body  Forces  In  Plates 

Equations  (2.102)  to  (2.107)  will  he  the  basis  for  the  later  chanters. 
Their  derivations  have  been  based  on  several  assumptions.  Except  for  the 
undeformable  normal  of  the  plate  and  the  constant  eddy  current  density  and 
constant  normal  component  of  the  induced  mapnetlc  field  across  the  thick¬ 
ness  of  the  plate,  every  variable  has  been  referred  to  the  mid-surface  of 
the  plate.  Specifically,  the  magnetic  force  density  has  been  treated  as 
constant  across  the  thickness  and  equals  in  maqnltude  to  its  value  at  the 
mid-surface.  In  these  calculations,  the  equation  has  also  been  linearized 
by  dropping  the  effect  of  the  Induced  maqnetlc  field.  In  this  subsection, 
the  implications  of  these  assumptions  will  be  considered. 

The  assumptions  of  constant  current  and  constant  normal  component  of 
the  Induced  maanetic  field  across  the  thickness  averaqe  the  Induced  current 
distribution.  In  the  calculation  of  the  normal  force,  the  tangential 
component  of  the  applied  field  is  evaluated  at  the  mid-surface,  and  the 
contribution  from  the  induced  field  B^,  is  not  Included,  Figure  2. 6a. 

Since  the  tangential  component  of  the  induced  field  will  vary  approximately 
linearly  across  the  thickness  and  have  zero  value  at  the  mid-surface,  the 
normal  force  pf*  corresponding  to  it  will  have  the  same  variation  and  be 
self-egull Ibratlng.  The  result  for  the  normal  force  should  therefore  be 
a  good  representation  for  the  total  force  across  the  thickness. 

The  same  approach  is  used  for  the  calculation  of  the  in-plane  force. 
Again,  good  results  should  be  obtained  for  the  total  force  by  averaging 


across  the  thickness.  Here,  however,  the  variation  of  the  normal  component 
of  the  applied  field  B°  olves  a  variation  of  the  In-plane  force  across  the 
thickness,  which  creates  a  body  moment  distributed  over  the  whole  plate, 
Fiq.  2.6b.  This  effect  Is  not  Included  in  the  present  model  since  all  the 
variables  have  been  evaluated  on  the  mid-surface.  The  hodv  moment  can  be 
easily  Included,  however.  Moon  (Ref.  25)  has  derived  the  followlnn  equa¬ 
tion  for  the  plate  which  includes  such  an  effect,  but  neolects  the  in-plane 
force  in  the  plate. 


(2. 108) 
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When  the  plate  is  thin,  and  the  drivina  coil  far  away  from  the  plate,  the 
applied  magnetic  field  will  have  only  a  small  variation  across  the  thickness. 
The  magnetic  moment  in  such  cases  can  be  nealected  without  Introducing  much 
error. 

A  coupling  effect  not  considered  In  the  present  model  comes  from  the 

rotational  velocity  of  the  plate  normal.  Fin.  2.6  (c).  The  rotational 

velocity  gives  an  effective  tangential  current 
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in  a  normal  direction  across  the  thickness.  This  causes  a  tangential  body 


force  distribution 


Inteoratlna  over  the  thickness,  one  obtains  the  body  moment  distribution 
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which  is  opposite  to  the  rotational  velocity  and  thus  has  a  rotational 
damplno  effect.  This  Is  the  "classical"  coupllno  term  reported  by 
Ambartsumyan  et  al .  (Pef.  36).  To  consider  this  effect,  one  needs  to  add 
the  rotational  Inertia  force  term  In  the  equation  of  motion.  In  the 
present  studv,  one  Is  limited  to  the  stretchlna  and  the  transverse  flexural 
deflection  of  the  plate.  The  rotational  effect  Is  not  considered. 

A  visualization  of  the  maanetic  body  force  In  an  Infinitely  lono  plate 
may  be  obtained  from  Faraday's  concept  of  the  maanetic  llnp  of  force.  The 
maanetic  line  of  force  Is  the  maanetic  field  line  to  which  certain  amount 
of  electromotive  tension  has  been  applied.  The  force  depends  solely  upon 
the  strenqth  of  the  maanet  producina  It.  A  compressive  force  Is  also 
developed  between  the  field  lines  and  Is  at  a  rlnht  anale  to  them  (Pef.  50). 
When  the  line  of  force  Is  cut,  due  to  either  the  time-variation  of  the 
field  or  the  motion  of  the  plate,  exactly  the  same  amount  of  force  Is 
released.  The  mechanical  stress  must  then  be  developed  to  equilibrate 
the  above-mentioned  tension  and  pressure. 

In  Fla.  2.7 a,  the  applied  field  Is  varvlna  at  a  relatively  low  fre¬ 
quency.  Eddy  current  Is  Induced  In  the  plate  and  the  aenerated  self 
field  Is  added  to  the  applied  field.  The  total  maanetic  field  lines  are 
distorted  from  straloht  but  still  penetrate  throuoh  the  plate.  The  tension 
and  pressure  associated  with  the  maanetic  field  lines  are  exerted  on  the 
plate.  It  Is  to  this  case  that  the  theory  developed  In  this  chapter  Is 
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applied  under  the  two  basic  maanetoelastlc  assumotlons. 

When  the  frequency  of  the  applied  field  Is  hlnh,  the  qenerated  self 
field  completely  screens  the  applied  field  from  penetratlno  the  plate. 

This  situation  Is  shown  In  Flo.  2.7b.  The  skin  depth  In  this  case  Is 
limited  to  a  very  thin  layer  on  the  top  and  bottom  surfaces  of  the  plate, 
and  the  eddy  current  Is  basically  a  surface  current  on  the  plate.  The 
same  situation  may  be  obtained  at  low  frequencies  for  conductors  with  very 
larqe  conductivity  (a  %  -).  The  maanetlc  force  In  this  case  acts  like  a 
surface  force.  Boundary  conditions  Eos.  (2.13a)  and  (2.13b)  should  be 
used  for  the  electromaonetlc  field  conditions  In  such  problems.  The 
maanetlc  force  Is  usually  called  the  maanetlc  pressure  In  such  cases. 


Chapter  3 


EDDY  CURRENT  CALCULATION  ON  RIGID  CONDUCTORS 

The  qoverrrino  equations  derived  In  Chapter  Two  for  the  calculation 
of  eddy  currents  are  applied  to  rlold  conductlnq  plates  In  this  chapter 
as  a  first  step  In  the  study  of  the  coupllno  of  eddy  currents  and  defor¬ 
mations.  The  one-and  two-dimensional  equations  for  steady  state  eddy 
current  problems  are  nondlmenslonal Ized  and  solved  by  the  finite  element 
method.  The  role  of  the  nondlmenslonal  parameters  In  these  equations 
are  Investlqated.  Finite  element  transient  analysis  Is  also  applied  to 
the  one-dimensional  problems.  The  results  of  the  analysis  are  verified  by 
comparlnq  with  analytic  solutions  In  the  hioh  and  low  frequency  extreme 
cases  and  with  experimental  data.  The  capability,  Input,  and  output  of 
the  computer  proqrams  developed  are  briefly  described. 

A  literature  survey  of  the  eddy  current  calculation  on  conductlnq 
plates  Is  qlven  in  the  first  section.  The  Intimate  relation  between  the 
quasi-static  field  equations  and  the  circuit  equations  Is  the  basis  of 
many  studies  and  provides  another  point  of  view  to  the  stream  function 
formulation  of  the  eddy  current  calculation  used  In  this  work.  This 
subject  Is  discussed  In  the  second  section  of  this  chapter. 

3.1  Literature  Survey 

Because  of  the  Importance  of  the  eddy  current  calculation  In  the 
present  research,  a  brief  literature  survey  on  this  subject  Is  olven  In 
this  section.  The  purpose  of  the  survey  is  to  provide  a  backoround  for 
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the  comparison  of  the  stream  function  method  with  other  approaches  to  this 
problem.  Some  key  features  of  eddy  current  calculation  and  recent  advances 
are  also  hlqhllqhted. 

As  discussed  In  Chapter  Two,  eddy  currents  are  found  In  any  conductlnci 
material  which  Is  subjected  to  tlme-varyinq  maanetlc  field.  Its  existence 
Is  manifested  by  (1)  Induced  Joule  heatlna,  (?)  a  maonetlc  reaction  field, 
and  (3)  a  maonetlc  force  resultlno  from  the  Interaction  of  the  Induclno 
and  reaction  fields.  It  Is  possible  to  put  the  eddy  currents  to  oood  use. 
Examples  are  Induction  furnaces,  which  use  the  ohmic  losses  of  hloh  fre¬ 
quency  eddy  currents;  electromaanetlc  nondestructive  testlno,  which  uses 
the  maonetlc  reaction  field  of  eddy  current;  and  maonetlc  formlno  and 
levitation,  which  use  the  maonetlc  force  oenerated  by  eddv  currents.  It 
Is  sometimes  necessary  to  take  steps  to  reduce  the  effect  of  eddy  currents. 
An  example  is  the  need  to  minimize  losses  to  enable  the  cores  of  trans¬ 
formers  and  rotatlno  electrical  machines  to  carry  the  reoulred  maonetlc 
flux;  one  wav  this  Is  achieved  Is  by  lamlnatlno  the  cores.  Studies  of  the 
effects  of  eddy  currents  may  be  found  in  the  technical  literature  related 
to  all  of  these  appl Icatlons. 

Rlqorous  analytic  solutions  to  eddy  current  problems  can  seldomly  be 
obtained  and  are  usually  limited  to  infinitely  ex*0  Hed  conductors  or 
closed  conductor  shapes  such  as  spheres  or  r"  os  fanv  such  problems 

and  some  solutions  may  be  found  In  Smythe  (Ref.  51).  Conductors  with 
finite  lenqth  and  more  aeneral  shape  are  more  difficult  to  handle.  The 
solutions  are  usually  obtained  by  either  analytic  technloues  such  as 
Iterative  approximation  and  Rayleioh-Reltz  method,  or  by  numerical  means 
such  as  finite  difference,  finite  element,  and  boundary  Inteoral  equation 
methods.  The  choice  of  the  solution  method  depends  to  a  laroe  extent  on 


the  formulation  of  the  problem.  The  maonetlc  enerov  and  variational 
principle  are  the  basis  of  many  studies.  The  circuit  analoov  provides 
another  model  more  Intuitively  appeallnn  to  electrical  ennlneers. 

The  analytic  technloues  have  the  advantaae  that  the  problem  parameters 
appear  as  variables  in  the  solutions.  The  effpct  of  alterina  one  or  more 
such  parameters  can  therefore  be  fairly  readily  appreciated.  The 
numerical  methods  have  wider  application  but  the  dlsadvantaoe  that  the 
parameters  are  concealed  In  the  numerical  results  produced  from  a  olven 
set  of  data.  Dimensional  analvsls  in  this  case  can  be  used  to  reduce  the 
excessive  computation  with  numerous  sets  of  data.  Some  illustrations  of 
thp  application  of  dimensional  analysis  in  the  present  study  are  olven  In 
sections  3  and  a  of  this  chapter,  A  brief  introduction  to  the  use  of 
dimensional  analysis  In  eddy  current  problems  Is  olven  In  Stoll  (Ref.  5?); 
this  Is  one  of  the  few  references  which  deals  with  this  specific  application 
of  dimensional  analysis. 

Rooks  and  monooraphs  addresslno  the  analysis  of  eddy  currents  arp 
scarce  In  the  literature.  Two  monooraphs  on  this  subject  arp  Stoll  (Ref. 
53)  and  Lammeraner  and  Stasl  (Ref.  53).  The  problems  discussed  are  mostly 
two-dimensional  ones  solved  bv  analytical  methods  or  finite  difference 
technloues.  The  application  of  finite  element  and  boundary  Intearal 
equation  methods  In  this  area  has  recently  been  developino  rapidly.  A  new 
monoaraph  edited  bv  Chari  and  Silvester  (Ref.  5A)  contains  many  discussions 
of  these  methods  In  eddy  current  problems. 

The  eddy  current  problems  studied  most  In  the  literature  are  the  two- 
dimensional  ones.  Taklna  rectanoular  coordinates  to  Illustrate  these, 
one  may  have  problems  In  which  (1)  the  currents  flow  In  the  z-dlrectlon 
only,  and  Jz  ,  Hx  ,  and  Hy  are  functions  of  x  and  y;  and  (?)  thp  maonetlc 


field  has  a  slnale  component  In  the  z-dlrectlon  only,  and  H  ,  J  ,  and  J 
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are  functions  of  x  and  y.  In  the  first  case.  It  Is  convenient  to  use  the 

maonetlc  vector  potential  which  Is  parallel  to  the  current  density  to 

formulate  the  problem.  The  maonetlc  vector  potential  £  Is  defined  as 

?  *  £  =  £  (3.1a) 

v  •  £  ■=  0  (3.1b) 

In  case  (1),  A  ■*  A  =0  and  A  Is  related  to  current  density  throuoh 
XV  z 

(Ref.  52) 

\  (V  E  f  f  £rJz  ln  Ijc-jc' I  d  (area,)  (3-?) 

area 

A  two-dimensional  diffusion  equation  may  then  be  derived  for  A  with 
appropriate  boundary  conditions. 

Silvester  and  Haslam  (Ref.  55)  first  presented  a  finite  element 
Galerkln  analysis  of  eddy  current  fields  In  maonetotellurlc  problems  usino 
such  a  formulation.  Chari  (Ref.  56)  showed  a  similar  formulation  usino  a 
variational  approach  for  the  eddy  current  problem  In  maonetlc  structures. 
Since  the  maonetlc  vector  potential  must  be  calculated  both  inside  and 
outside  the  conductor,  the  discretization  of  space  must  be  extended  to 
reolons  sufficiently  far  away  from  the  conductor.  Chari's  study  concerns 
the  hlqhly  permeable  iron  part  of  the  electrical  machines  where  slonlflcant 
skin  effects  and  ohmic  loss  are  developed.  The  eddy  current  problem  on 
Infinitely  lonq  plates  In  the  present  research  belonos  to  the  same  class 
of  two-dimensional  problems  but  Is  directed  to  conductfno  materials  with 
small  relative  permeability.  For  such  material  the  skin  effect  Is  not 
well  developed  In  the  low  ana  Intermediate  frequency  ranoe.  The  stream 
function  method  used  has  the  advantaoe  of  llmltino  the  discretization  to 
the  conductor  body  only,  at  the  pxnense  of  solvina  a  system  of  linear 


equations  with  a  full  coefficient  matrix. 

For  the  case  In  which  the  maonetlc  field  has  only  a  slnale  component 
H  ,  manv  authors  use  the  electrical  vector  potential  T  to  formulate  the 
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problem.  The  electrical  vector  potential  T  Is  defined  such  that 

v*T-^  (3.3) 

For  strictly  two-dimensional  problems,  T  has  only  one  nonzero  component 

T  ,  which  equals  H  .  A  two-dimensional  diffusion  equation  for  T  can 
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aciain  be  derived.  The  problem  may  then  be  solved  either  by  finite  element 
method  throuqh  a  variational  formulation  (Ref.  57),  or  by  boundary 
Intearal  equation  method  throuoh  a  transformation  of  the  differential 
equation  Into  an Inteoral  equation  (Ref.  58). 

Silvester  (Ref.  5^)  developed  the  useful  concept  of  eddy  current  modes 
from  the  eloenfunctlon  expansion  of  the  solution  for  this  class  of  problems. 
The  orthoaonal Ity  of  the  elaenfunctlons  enables  an  equivalent  circuit 
representation  of  the  eddy  current  distribution  to  be  developed.  Fach 
eloenfunctlon,  or  mode,  forms  a  separate  R-L  circuit  and  there  is  no 
maqnetlc  coupllno  between  the  modal  circuits.  In  a  subseouent  paper  (Ref. 
60),  Silvester  extends  the  modal  network  theory  to  Infinitely  Iona  flat 
conductlnq  plates.  The  current  density  Jz  Is  used  to  formulate  the  prob¬ 
lem,  and  consideration  Is  limited  to  the  conductor  Itself.  The  plate  Is 
subdivided  Into  N  parallel  strips  of  equal  width.  Fach  strip  Is  assumed 
to  carry  a  uniformly  distributed  current  at  all  times.  The  thickness  of 
the  plate  comes  Into  the  calculation  of  the  resistance  R  of  the  strips. 

The  Inductances  are  calculated  uslnq  the  qeometrlc  mean  distance  theory 
of  linear  conductors.  The  solutions  are  then  found  bv  sol  vino  the  matrix 
circuit  equations.  This  eddy  current  modal  theory  has  later  been  applied 
to  the  studies  of  many  othpr  eddy  current  problems  (Refs.  61,  6?). 
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For  the  case  In  which  the  pddv  current  flows  in  ’-direction  only, 
Silvester  (Ref.  63)  developed  an  inteqral  eouation  essentially  of  the 
Fredholm  type  which  treats  current  density  directly.  This  eouation  has 
been  solved  with  various  methods  by  Schaffer  and  Randevet  (Ref.  64), 
bv  Silvester  (Ref.  65),  by  Gopinath  and  Silvester  (Ref.  66),  and  hv 
Silvester,  Wona  and  Burke  (Ref.  61).  Trowbrldoe  (Ref.  54)  and  Blddlecombe 
et  al.  (Ref.  67)  also  developed  an  intearal  equation  formulation  for  eddv 
currents  usino  the  maonetic  vector  potential. 

For  two-dimensional  eddv  current  problems  in  cylindrical  coordinate 
systems,  the  only  nonzero  components  of  current  and  the  maonetic  vector 
potential  are  in  the  ci rcumferential  direction.  A  two-dimpnsi ona 1 
diffusion  eouation  can  be  constructed  for  the  maonetic  vector  potential. 
Donea,  Giuliani  and  Philippe  (Ref.  6P)  presented  a  finite  element  formula¬ 
tion  usinq  a  variational  method.  Becker  and  Pillsburv  (Pef.  ?A) ,  and 
Miya,  et  al.  (Ref.  2P)  solved  it  by  usino  a  finite  element  Galerkin 
technique.  Both  studies  consider  counled  electromannetic-mechanical 
effects.  Many  electromaonetic  NOT  studies  also  consider  problems  with 
cylindrical  symmetry.  Most  of  them  use  a  variational,  or  enerov  functional, 
formulation.  Fxamples  are  Palanisamy  and  Lord  (Ref.  6R),  Nphl  and  Pemerdash 
(Ref.  7D),  and  Hwano  and  Lord  (Pef.  71). 

All  the  two-dimensional  problems  discussed  assume  the  conductors  to 
be  infinitely  lono.  For  conductor  plates  with  finite  lenoth,  the  eddv 
current  problem  becomes  three-dimensional.  To  simplify  the  problem,  it  is 
usually  assumed  that  the  skin  depth  is  larqe  compared  to  the  thickness  so 
that  the  current  density  can  be  taken  as  constant  across  the  thickness. 

The  network  analoov,  in  many  respects  eoulvalent  to  eddv  current  modal 
theory,  plays  an  Important  role  In  such  problems.  This  analoov  will  be 


discussed  in  more  detail  in  the  subseouent  section  of  this  chanter. 

Turner  (Refs.  78,  73,  74)  has  developed  cnmouter  proorams  for  two-  and 
three-dimensional  transient  eddy  currents.  The  conduction  materials  are 
represented  hv  a  network  of  current-carryi no  line  elements  in  order  that 
Maxwell's  field  conations  mav  be  replaced  hv  KirchofC's  circuit  rules.  Loon 
currents  are  used  as  the  variables  in  the  analysis  so  that  the  Kirchoff's 
mode  eouatlons  are  satisfied  automatical! y.  Thin  conducting  platps  and 
shells  are  treated  as  two-dimensional  problems  bv  such  a  network  model. 

Ko'-earl  and  Suzuki  (Ref.  75)  developed  a  finite  element  circuit  method 
for  the  analysis  of  transient  eddy  currents  on  thin  conduction  plates  and 
shells.  The  electric  vector  potential  is  used  to  formulate  the  problem. 
Unlike  the  two-dimensional  problems  such  as  lono  prismatic  bars,  the 
electric  vector  potential  now  is  not  eoual  to  the  mannetic  field  but  differs 

by  an  additional  scalar  potential  !■,  i.e. 

B  r  u  T  +  v  ifi  (7.4) 

'V  “V 

A  set  of  circuit  eouations  are  formulated  in  this  method  from  the  balance 
of  the  mannetic  enernv  and  Joule  loss  of  the  system.  Eddy  current  modal 
analysis  is  then  performed  throuob  an  eioenvalue  analysis  of  the  derived 
finite  element  circuit  eouations.  This  formulation  oives  a  set  of  linear 
equatir  'ith  fully  symmetric  matrices.  The  effect  of  the  thickness  comes 
into  the  .alculation  of  the  resistance  matrix.  The  inductance  and  resist¬ 
ance  matrices  in  ehis  method  are  all  calculated  by  numerical  inteoration. 

Compter  and  Hamels  (Ref.  76)  studied  a  coupled  mechanical-eddv  current 
problem  for  a  movable  conduction  disk  usinn  a  network  circuit  model. 
Carpenter  (Pefs.  77,  78)  described  a  finite  element  network  model  and  its 
application  to  eddy  currents  in  thin  plates.  Thus  most  pddv  current  studies 
for  thin  conduction  plates  use  a  circuit  model  for  analysis.  The  applica- 


t i on  of  variational  methods  Is  difficult  because  the  dissipation  term 
affects  the  variational  formulation.  Hammond  and  Penman  (Ref.  7R)  studied 
the  upper  and  lower  bounds  of  the  resistances  and  Inductances  of  thick  and 
thin  conductina  plates  uslno  a  variational  approach,  but  the  problems  thev 
treated  are  for  the  classical  skin  effect  of  current  carryino  conductors, 
not  the  Induction  problems  treated  In  the  present  research.  De  Me.v  (Ref. 
80)  presented  a  variational,  Inteoral  eouatlon  study  of  eddy  currents  In 
plates  uslna  a  stream  function  representation.  But  he  neolects  the 
reaction  of  the  eddy  currents  on  the  maonetlc  field  and  therefore  only 
obtains  a  local  solution  of  the  problem. 

3.2  Quasl-statlc  Field  Equations  and  Circuit  Equations 

As  discussed  In  the  last  section,  there  Is  an  Intimate  relation 
between  the  circuit  model  and  the  various  formulations  for  eddy  currents 
In  thin  plates.  This  relation  Is  examined  in  this  section  throuah  the 
correlations  between  the  quasl-statlc  Maxwell's  equations  and  the  rules 
of  circuit  theory.  The  comparison  provides  a  point  of  view  alternative 
to  the  stream  function  method  used  in  this  study. 

The  scalar  potential  4  and  the  vector  potential  £  are  first  Introduced 
for  time  varylnq  fields.  They  are  defined  as  (Ref.  32) 

£  =  v  x  £  (3.1a) 

<3-5> 

and 

(Lorentz  condition)  (3.6) 

These  relations  come  from  Equations  (2.24).  The  electromaqnetlc  field  can 
be  completely  described  by 
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V?  k  '  11  c  "  "  ll4 

n 


(3.7a) 


2  x  3  .  o 

V  iji  -  u  £  - S-  4>  =  -  — 

3t 


(3.7b) 


Neolectinq  retardation  in  the  quasi-static  case,  one  obtains  the  oeneral 
solution 


fi,fc/VHdV 


<J>  =  r  —  p  dV 
v  4  tt  c  V  r 


(3.Ra) 


(3.8b) 


The  relations  between  quasi-static  field  equations  and  circuit  equations 
will  be  discussed  in  terms  of  these  two  potentials. 

The  circuit  theory  includes  the  two  followino  Klrchoff's  laws: 

(1)  the  current  law:  At  any  DOint  in  a  circuit,  the  current  flow  is 
equal  to  the  rate  at  which  charge  is  passino  throuoh  a  cross-section 
of  the  conductor  at  that  point; 

(2)  the  voltaqe  law:  Around  a  closed  circuit,  the  alnebralc  sum  of  the 
potential  difference  Is  zero. 

These  two  laws  may  be  put  into  the  following  equation  form  when  there  Is 
no  capacitor  in  the  circuit,  Fiq.  3.1,  which  is  the  case  in  the  present 
study  since  eddy  current  arclnq  is  not  considered. 

(3.9a) 


1  !in  =  1  W 
dl 


(3.96) 


e  '  IR  +  L  Ht 

In  which  I  is  the  total  current  across  the  cross  section  of  the  conductor, 
e  is  the  applied  electromotive  force;  R  is  the  resistance  per  unit  lenqth; 
and  L  the  Inductance  of  the  circuit. 


!R 


Figure  3.1  Metallic  conductor  wire  showing  path  of 
Integration  lying  on  surface  and  the 
equivalent  RL  circuit. 


Figure  3.2  Mutual  Inductance  between  circuits. 

Variation  of  jJ'dv'  Induces  an  emf  In  d£ 


Eq.  (3.9a)  comes  from  the  equation  of  continuity  of  current  anrl  Is 
automatically  satisfied  bv  the  potentials  jA  and  $  defined  before.  The 
relation  between  Fq.  (3.9b)  and  the  ouasl-static  field  equations  Is  dis¬ 
cussed  below. 

A  circuit  Is  defined  as  any  closed  path.  For  the  one  shown  In  Flo. 

3.1,  one  first  separates  the  total  electric  field  £  In  the  conductor  Into 
two  parts:  E’  Is  that  due  to  charaes  and  currents  defined  by  Fq.  (3.5), 
and  is  the  impressed  field  due  to  enerqy  sources  such  as  the  external 
excltlnq  coll. 

Then 

f  ■  i  -  S' 

■U*hk*v*  (3-11? 

Inteoratlno  around  the  circuit  alona  a  nath  on  the  surface  of  the  conductor 
one  has 

)  f  •  dt  -  |  y  i  ■  H  *  (  ft  t  •  "i  *  f  •  «4  I3-’?) 

The  term  on  the  left-hand  side  of  the  above  equation  Is  the  applied  electro¬ 
motive  force,  which  Is  desianated  e  as  In  Fq.  (3.9b).  The  last  term  Is  zero 
since  without  retardation  one  has 

(i  v<p  •  d£  =  |  y|-  dt  -  j  dt  =  o  (3.13) 

The  second  term  on  the  rlqht-hand  side  of  Eq.  (3.12)  may  be  written  as 
^  |  £  •  dt  for  a  rlqld  stationary  conductor,  in  which  £  Is  determined  by 
Eq.  (3.8a)  in  the  quasi-static  case.  Assume  that  the  conductor  Is  a  wire 
with  small  cross-sectional  dimensions  so  that  dV’  In  Eq.  (3.8a)  may  be 


replaced  by  I  d£'  where  I  Is  the  total  current  across  the  cross  section  of 
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the  wire.  One  has 

£  =  $7  f  f  d*' 


and 

d 

HF 


£  •  d* 


dl 

3F 


In  which 


L  =  £- 
4  TT 


|  f f  <>r 


dt 


(3.14) 

(  .  .5) 


is  the  inductance  of  the  circuit  and  Is  a  function  only  of  the  oeometrv  of 
the  circuit.  Another  way  of  defininq  L  is  the  followlno: 

LI  ■  f  t  •  1* 

•  U  7  -  4  •  d  i 

’'st-U  I3-1*) 

or 

L  -}/s  5  •  d  {  (3.17) 

In  which  (3  Is  the  total  maqnetlc  Induction. 

The  first  term  on  the  rlqht-hand  side  of  Eq.  (3.12)  Involves  J /? 

'V 

which  represents  the  total  electric  field  £  at  the  surface  of  the  conductor. 

For  sinusoidal  time  variations  Is  related  to  the  total  current  I  by 

4  l  (3.18) 

c 

in  which  Ac  Is  a  complex  number  havlnq  the  dimensions  of  area,  but  not 

actually  equal  to  the  cross-section  area  of  the  conductor.  It  Is  complex 

because  the  current  density  at  the  surface  of  a  conductor  Is  not  In  phase 

with  the  total  current  when  a  skin  effect  exists.  In  such  cases  1/oA  , 

c 

deslqnated  by  Z,  Is  called  the  Impedance  (per  unit  lenqth)  of  the  conductor. 

When  the  frequency  is  sufficiently  low  that  the  skin  effect  Is  not 
well  developed,  the  current  Is  essentially  uniform  over  the  entire  cross- 
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section  of  the  conductor.  In  such  cases  Ac  Is  the  actual  cross-sectional 
area  A  of  the  wire  and 


"  7  i  ■  “l  ‘  |  3T  d<  ‘  IR 

In  which 


R  = 


1 

0  off 


dt 


(3.19) 


(3.20) 


Is  the  usual  resistance  per  unit  lenqth  of  the  wire. 

The  ouasl-static  field  equations  are  therefore  related  to  the  equili¬ 
brium  equation  (3.9b)  for  a  closed  circuit.  In  the  Induction  problem,  the 
emf  e  Is  Generated  by  the  varylnq  maonetlc  field  due  to  the  excltinq  coll, 
Fiq.  3.2.  Assume  that  the  Induced  current  In  the  excltinq  coll  b  is  zero, 
l.e.,  t:iat  1°  Is  voltaoe-controlled  or  that  the  excltinq  coll  has  Infinite 
conductivity.  One  has 

«  *  -  f  h  1°  ■  st  >5-21 

in  which  jA°  Is  the  component  of  vector  potential  due  to  the  maanetic  field 
of  the  excltinq  coll  and  Is  olven  by 

«°  *  St  \  ?  >i0  dVh 


Substltutlna  Ea.  (3.2?)  into  Ea.  (3.21),  one  has 


e  =  -  M 


dl° 

W 


In  which 


dt 


(3.22) 


(3.23) 


(3.24) 


Is  the  mutual  Inductance  between  the  colls  a  and  b. 
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The  relationships  established  above  for  a  sinole  discrete  current  loop 
use  exactly  the  same  assumptions  applied  in  the  present  work  for  a  thin 
conducting  plate,  l.e.,  the  skin  effect  Is  not  slanificant  and  the  current 
Is  constant  across  the  thickness  of  the  plate.  A  further  similarity  occurs 
between  the  circuit  equation  that  describes  the  discrete  circuit  and  that 
for  the  eddy  current  modes  on  the  conducting  plates.  The  eddy  current 
modes  are  certain  distribution  patterns  of  the  current  density  that  are 
electrically  and  maanetlcally  uncoupled  from  each  other.  The  modal  distri¬ 
bution  can  be  obtained  by  decomposina  the  total  current  distribution  Into 
components  orthonormal  to  both  resistance  and  inductance  operators.  The 
derivation  Is  briefly  described  below. 

The  law  of  conservation  of  magnetic  energy,  or  the  Poyntlng's  theorem. 


for  the  Induction  problems  states  that  (Ref.  5?) 

hii  dU 

dll  .  u  _  m 

+  w  ‘  3T 


(3.25) 


In  which  U  Is  the  total  maqnetlc  energy  due  to  the  current  distribution, 

W  Is  the  total  ohmic  loss,  and  U  the  mutual  maqnetic  eneray  between  the 

m 

external  maanetlc  field  and  the  induced  current.  Assumlno  constant  current 
across  the  thickness  of  the  plate,  one  has  (Ref.  40) 


H_//l 

8ir  S  S  rP0 


4p  '  4o 


(3.26a) 


W  =  /  n  ^  ^  dS  (3.26b) 

Um  =  /  ^  dS  (3.26c) 

In  which  n  Is  the  area  electrical  resistivity,  £°  Is  the  magnetic  vector 
potential  of  the  externally  applied  field,  and  ^  Is  the  total  induced 
current  density  across  the  thickness  of  the  plate. 
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By  uslna  an  electric  vector  potential  T  for  the  current  density, 
and  the  finite  element  technique  Kamearl  and  Suzuki  (Ref. 75)  obtained  the 
followlno  coupled  circuit  eouatlons  for  the  plate 

[L]  ^  (VI  +  [R]  {VI  =  ^  {el  (3.27) 

In  which  [L]  and  [R]  are  the  time- Independent  Inductance  and  resistance 
matrices,  respectively.  Both  [Ll  and  [R]  are  positive-definite  and  fully- 
populated;  however,  for  a  uniform  mesh  they  are  also  symmetric.  Apnlvlna 
the  aenerallzed  eloenvalue  analysis  to  the  coupled  system  of  equations, 
Kamearl  and  Suzuki  obtained  the  eddy  current  modes  (FI,  n  =  1,2,  ....  N, 


where  N  Is  the  total  number  of  the  nodal  points.  Since 

[L]  (En>  -  *n  [Rl  (Enl  (3.2ft) 

one  has 

<VT 

<E/  tLl  (En)  -  »„  (W»> 

The  modes  are  therefore  electrically  and  maanetlcally  uncoupled.  By 

expandlnq  the  potential  T  by  the  eddy  current  mode, 

N 

{T)  =  Z  Pn  (Enl  (3.30) 

n*=l 


one  obtains  the  uncoupled  equation  for  each  mode 

n*l ,2,  ....  N 

II  uv  II 

In  which 


Pn  +  Xn  Pn  "  Ht  en 


e  =  { E  }  {el 
n  n 


(3.31) 


(3.32) 


The  modal  equations  (3.31)  are  electrically  and  maonetlcally  uncoupled  to 
each  other  and  are  linked  maqnetlcally  to  the  drivlnq  field  throuqh  Eq. 


(3.3?).  Fa.  (3.31)  Is  completely  equivalent  to  the  circuit  eouation  (3.%) 
with  resistance,  inductance,  and  decay  time  constant  1,  Xp,  and  x  ,  respec¬ 
tively.  The  full  analoov  between  the  circuit  theory  and  the  ouasi-static 
description  of  eddy  currents  on  thin  plates  is  thus  completed. 

3.3  Eddy  Currents  on  Infinitely  Long  Plates 

The  finite  element  Galerkin  technioue  is  applied  to  solve  steady 
state  and  transient  eddy  current  equations  for  infinitely  lono  conductinq 
plates.  The  formulations  and  the  Fortran  codes  developed  are  described 
in  subsection  3.3.1.  To  facilitate  parameter  studies  for  this  problem, 
the  steady  state  eddy  current  equation  is  nondimensional Ized.  A  study  of 
the  relevant  nondimensional  parameters  is  presented  in  subsection  3.3.?, 
toaether  with  other  examples  which  provide  analytical  and  experimental 
verifications  of  the  numerical  results. 


3.3.1.  Finite  Element  Formulations  for  Steady  State  and  Transient  Problems 
The  transient  eddy  current  equation  is  aiven  by  Ea.  (2.104),  with¬ 
out  the  couplino  terms. 


an 


1 
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HI 

n 


dX '  = 


(3.33) 


in  which  X  Is  used  to  denote  the  position  in  physical  units,  and  l  is  the 
total  width  or  span. 


3. 3.1. a  Formulation  and  Implementation  for  Steady  State  Problems 

For  steady-state  harmonic  currents  in  the  plate,  Eo.  ( 3 . ^ 3 )  may 


be  nondimensional Ized  into 


hijuii  Jilin  fin  rm 


V 1  ■  V  * 


11 


-  iR*  +  l^f  fl  - -‘-g--^  dr.  =  1RR(x)  (3.34) 

dx  "  J0  U-x)  +  j 

in  which  x  and  £  are  distances  across  the  width  nondimenslonal ized  with 

respect  to  the  thickness  h;  *  =  <tl°e^u)t;  R°  =  (u/a)8(x)  I°e^u,t ;  the 

2 

maonetic  Reynolds  number  R  -  much  ,  which  is  related  to  the  skin  depth  6 
throuoh  R  *  ?(h/<s)^;  Ic  is  the  current  in  a  nearbv  exciter  coil;  and  ! 
is  the  nondlmensionali zed  width. 

The  finite  element  (FFl  Galerkin  method  is  used  to  solve  the  inteoro 
differential  equation  (3.’dl.  t  is  approximated  nloballv  and  locally  by 
piecewise  linear  models 
fi 

$  =  r  M.  (3.35a) 

k*l 


=  T.  N*  *  (3.35b) 

k*l 

in  which  G  is  the  total  number  of  noda  i  points,  the  superscript  E  denotes 
the  E-th  element,  and  are  the  qlobal  interpolation  functions  qenerated 
from  the  local  linear  element  shape  functions  NE.  The  linear  aloebraic 
equations  for  each  element  are 


2 

-  T 

k»l 


Sjk*k 


2 

1 

k=l 


PEL'tL 
jk  k 


+  i 


G 

7. 

k=l 


°Jk*k 


in  which 


(3.36) 
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dNE  dNE 
dx 


W  dx  ;  Pik  3  R 


Q- 


R 


jk  77 


Mk(e)Wj(r.)  dC  ;  re  =  R 


f  NEBdx 


(3.37) 


the  welohtino  functions  W^  used  in  calculatlna  the  inteoral  term  is 

Wfc.  (f)  =  |  -  1...  -T  dx  (3.38) 

-1  ^(r-x)‘+T 

Both  Eos.  (3.37)  and  (3.38)  are  intearated  analytically.  The  resultinq 
element  matrices  are  shown  in  Appendix  A  for  both  the  local  and  nonlocal 
terms.  The  assembled  olohal  matrix  is  complex,  oenerallv  nonsymmetrlc, 
and  fully  populated.  The  usual  limited  connectivity  (banded  matrix) 
finite  element  representation  is  lost  due  to  the  nonlocal  nature  of  the 
problem. 

When  the  frequency  of  the  external  maonetic  field  is  low,  the  effect 
of  the  self  field  is  small.  The  solution  in  this  case  will  approach  the 
local  solution  obtained  bv  dropoinn  the  second  and  third  terms  on  the  left 
hand  side  of  Eq.  (3.34).  When  the  frequencies  are  hiah,  the  distribution 
of  the  eddy  current  is  primarily  influenced  by  the  self  field.  The 
solution  will  approach  the  one  obtained  by  droooino  the  first  term  on  the 
left  of  Eq.  (3.34).  The  solution  so  obtained  is  an  asymptotic  value  for 
the  hiqh  frequency  limit,  and  Is  called  the  "imaoe  solution"  in  the 
1 iterature. 

A  Fortran  proqram  EDDY 1  has  been  developed  based  on  Eqs.  (3. 36 )-( 3. 38) . 
It  is  desioned  to  perform  the  followino  calculations: 

(1)  Steady  state  local  solutions  of  stream  function,  current, 
temperature,  and  pressure, 

(2)  Steady  state  nonlocal  solutions  of  stream  function,  current, 
temperature,  and  pressure,  and 

(3)  Steady  state  imaoe  solution  of  stream  function,  current, 
temperature,  and  pressure. 


i 


/') 

To  opfform  the  calculations  necessary  to  construct  a  spectrum,  EDDY  1  performs 
the  above  calculations  for  any  number  of  maonetlc  Reynolds  numbers  In  one 
run.  Multiple  load  cases  (up  to  in)  are  allowed.  Fach  load  case  will  he 
solved  for  all  the  maonetic  Reynolds  numbers  specified  for  that  run. 

Presently,  EDOY 1  can  handle  uniform  magnetic  field  cases  and  cases 
where  the  fields  are  qenerated  bv  any  number  of  conductor  wires  parallel 
to  the  plate.  The  strenoths  and  senses  of  the  currents  and  the  positions 
of  the  wires  may  be  different.  Space  has  also  been  left  to  include  other 
types  of  exciting  fields  in  the  prooram. 

Program  EDDY1  takes  the  geometric  Information  of  the  plate  and  magnetic 
field  source  information  as  input.  The  total  numbers  of  nodal  points,  load 
cases,  parameter  sets,  and  the  plottina  option  of  the  output  must  he 
specified  on  the  master  input  card.  Nodal  point  information  mav  be  opnerated 
for  portions  of  the  mesh  that  are  uniform. 

The  type  of  source  of  maonetic  field  is  specified.  No  other  informa¬ 
tion  Is  needed  for  a  uniform  maonetic  field.  For  the  parallel  wire  cases, 
the  positions,  senses,  and  the  relative  strenoths  of  the  currents  In  each 
wire  must  be  Input  subsequently. 

All  the  Input  data  must  be  nondimensionallzed  accordlno  to  the  con¬ 
ventions  presented  In  the  formulation. 

Values  of  the  stream  function,  eddy  current,  temperature  Induced  In  a 
half-cycle  of  the  current,  and  time-averaoed  maonetic  pressure  exerted  on 
the  plate  are  produced  as  output  of  the  prooram.  The  stream  function  and 
eddy  current  are  calculated  In  complex  form.  The  stream  function  values 
are  calculated  at  the  nodal  points  of  the  mesh,  while  the  current,  tempera¬ 
ture,  and  pressure  are  evaluated  at  the  centroids  of  the  elements.  The 
stream  function  and  eddy  current  are  olven  in  complex  form,  and  the  modulus 


and  phase  anale  are  calculated  for  the  eddy  current.  Values  of  the  modulus 
and  phase  annle  for  different  freouencles  may  be  used  to  nenerate  the 
spectrum  of  the  current  any  any  point  on  the  plate.  These  curves  and  the 
spectrum  of  the  excltlno  current  may  be  used  to  calculate  the  transient 
current  at  the  point  by  fast  Fourier  transform  technloues. 

All  the  outputs  of  the  prooram  may  be  plotted  uslnn  a  printer-plotter 
subroutine.  This  capability  Is  optional  and  can  be  snecif1ed  In  the  Input. 

In  this  prooram  the  local  solutions  are  printed  out  and  plotted 
parallel  to  the  nonlocal  solutions  for  comparison.  When  therp  are  several 
load  cases,  the  results  are  printed  out  in  orouns  In  the  seouence  o f  the 
load  numbers.  In  the  multiple  parameter  (freouencv)  cases,  this  is  done 
for  each  Reynolds  number.  Title  lines  will  be  printed  for  each  parameter 


to  distlnoulsh  the  different  orouns  of  output.  Imaoe  solutions,  which  are 
Independent  o *  the  Reynolds  number,  are  printed  out  last  in  the  seouence 
of  load  cases. 


3 . 3 . 1 . b  Formulation  and  Implementation  for  Transient  Problems 

The  Halerkin  finite  element  (FF)  version  of  the  transient  eouation 
(3.33)  Is  formulated  in  a  wav  similar  to  the  steady  state  problem.  It  is 
Inteorated  with  respect  to  time  by  an  implicit  scheme.  ti>  Is  appro  lima  ted 
In  physical  units  oloballv  and  locally  by 


'!•  *  F  fY 

k  =  l 


E  F 

*  =  T.  N  i 

k=l  k 


(3.39a) 


(3.39b) 


The  linear  equations  for  each  element  are 
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where  the  weighting  function  Is: 


(3.40) 


(3.41) 


(X* ) 


t  h  (X) 


(3.42) 


Equations  (3.41)  and  (3.42)  are  again  Integrated  analytically,  and  the 
resulting  expressions  are  those  given  In  Appendix  A.  The  global  matrix 
form  of  Eq.  (3.40)  Is 


[A]  {*}  -  [S]  {*}  •  {R>  (3.43) 

In  which 

[A]  «  [01  -  [P]  (3.44) 

Again,  [A]  Is  full  and  generally  nonsymmetrlc. 

The  transient  equations  (3.43)  are  Integrated  using  the  following 
scheme 

4e  w  -  <’-•>  tsn 

-  (^  [A]  +  e  [S]>  (♦)t  +  e  (R)t  +  O-e)  (R>t+At  (3.45) 

e»  0.5  has  been  used  In  all  the  problems  analyzed. 


A  Fortran  program  EDDY1T  has  been  developed  based  on  Eqs.  (3.41 )- (3.45) 
and  has  been  built  Into  a  more  general  transient  analysis  program  EDDYBEAM 
for  the  coupled  electromagnetic-mechanical  effects.  The  purely  transient 
one-dimensional  eddy  current  analyzer  EDDY1T  may  be  recovered  from  It  by 
specifying  the  Input  parameter  LINK  to  be  zero.  EDOY1T  Is  designed  to 
perform  the  following  calculations  for  pure  eddy  current  problems: 

(1)  transient  local  solutions  of  stream  function,  current, 
temperature,  and  magnetic  force,  and 

(2)  transient  nonlocal  solutions  of  stream  function,  current, 
temperature,  and  magnetic  force. 

EDDY1T  handles  one  load  case  In  one  run  only.  The  two  types  of 
magnetic  fields  Included  are  a  time-varying  uniform  magnetic  field  and  a 
field  generated  by  any  number  of  conductor  wires  parallel  to  the  plate. 

The  strength  and  sense  of  the  currents  and  the  positions  of  the  wires  may 
be  different,  but  the  time  variations  of  all  the  currents  must  be  the 
same.  The  time  function  may  either  be  analytic  and  calculated  at  each 
time  step,  or  numerical  and  read  In  from  the  Input  data.  Presently  a 
half-sine  function,  a  continual  sinusoidal  function,  and  a  ramp  function 
of  time  have  been  Included  In  E0DY1T.  Flexibility  has  been  left  for  the 
user  to  create  additional  modules  If  other  analytic  forms  of  the  time 
function  are  preferred.  The  temperature  calculated  Is  the  Instantaneous 
temperature  rise;  the  effect  of  heat  conduction  Is  not  Included. 

Geometric  Information  of  the  plate,  source  Information  of  the  magnetic 
field,  relevant  physical  constants,  time  Integration  parameters,  and  the 
numerical  Information  (If  any)  describing  the  time  function  of  the  field 
must  be  Input  to  the  program  EDDY1T.  The  units  of  all  the  Input  values  must 
be  consistent.  The  rationalized  MKS  system  (Glorgl  system)  has  been  used 
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for  all  the  problems  In  the  present  research. 

Mesh  generation  and  Incremental  plotting  capabilities  are  also  Includ¬ 
ed  in  EDDYlT.TIme  history  plots  at  selected  points  and  spatial  variation 
plots  at  selected  time-steps  for  both  current  and  magnetic  force  may  be 
produced  In  the  same  scale  for  comparison.  Appropriate  mesh  generation 
and  plotting  Information  must  be  Input  If  these  capabilities  are  used. 

Values  of  the  stream  function,  eddy  current,  magnetic  force,  and 
temperature  at  each  time  step  are  produced  as  output  of  the  program.  The 
stream  function  values  are  calculated  at  the  nodal  points  of  the  mesh. 

The  eddy  current  and  temperature  are  evaluated  at  the  centroids  of  the 
elements  and  are  constant  within  each  element.  The  magnetic  forces  are 
calculated  at  the  Integration  oolnts  of  each  element  and  vary  within  each 
element  with  the  tangential  component  of  the  applied  field.  The  outputs 
may  be  plotted  using  a  library  subroutine.  A  short  form  of  the  output  Is 
provided  as  an  option  If  long  time  solutions  are  needed.  Local  and  non¬ 
local  solutions  are  handled  as  different  jobs  and  must  be  executed  separate¬ 
ly.  The  type  of  analysis  (local  or  nonlocal)  must  be  soeclfled  In  the 
Input. 

3.3.2  Numerical  Results  and  Experimental  Verifications  for  Infinitely 

Long  Plates 

The  programs  EDDY1  and  EDDY1T  have  been  applied  to  a  number  of  problems 
to  test  and  verify  the  programs'  validity  and  to  demonstrate  their  utility 
for  problems  of  scientific  and  practical  Interest.  A  portion  of  the 
verification  process  Is  the  comparison  of  the  computed  results  with  experi¬ 
mental  measurements.  A  summary  of  these  numerical  results  Is  given  below. 


3. 3. 2. a  Steady  State  Analysis  by  Propram  EDDY1 

The  first  test  of  the  numerical  techniques  discussed  In  this 
section  Is  the  calculation  of  Induced  currents  In  a  lonp  rectangular  plate 
by  a  two-dimensional  maqnetlc  field  uslnq  program  EDDY1, Figure  3.3.  Near 
the  center  of  the  plate.  Section  C-C  In  Figure  3.3,  the  end  effects  of  the 
plate  are  negligible,  and  the  currents  are  one-dimensional  or  are  parallel 
to  the  long  edges  of  the  plate.  The  external  magnetic  field  used  In  the 
test  cases  was  that  due  to  parallel  current  filaments  above  the  plate  and 
a  uniform  time  dependent  magnetic  field.  In  addition  to  the  Induced 
current  distribution  across  the  plate,  the  Induced  temperatures  and  maqnetlc 
pressure  distributions  were  calculated. 

The  following  objectives  were  met  with  the  program  EDDY1. 

a)  Comparison  of  low  magnetic  Reynolds  number  (R)  results  with 
direct  quadrature  of  local  theory. 

b)  Comparison  of  hloh  R  results  with  the  Imaoe  method. 

c)  Demonstration  of  the  Importance  of  nonlocal  theory  for 
moderate  frequencies  or  R, 

d)  Comparison  of  finite  element  calculations  with  experimental 
Infrared  measurements  reported  In  Refs.  101  and  102. 

e)  Combination  of  the  finite  element  (FE)  results  for  different 
field  frequencies  with  the  fast  Fourier  transform  (FFT)  to 
predict  dynamic  currents  and  pressures. 

f)  Demonstration  of  the  Importance  of  the  edge  effect  In 
Increasing  the  current  density  and  magnetic  pressure. 

g)  Calculation  of  the  magnetic  forces  on  a  Iona  plate  due  to 
a  tilted  coll. 
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h)  Examined  the  effect  of  a  banded  matrix  as  an  approximation 
to  the  full  matrix  or  complete  nonlocal  solution. 

1)  Demonstration  of  the  size  effect  of  the  width  on  the  Induced 
current  through  a  different  nondlmenslonal  magnetic  Reynolds 
number  Rw. 

Each  of  these  will  be  described  separately  below. 

a)  Comparison  of  FEM  with  Direct  Quadrature  for  R«l.  When  R«1 
one  may  drop  the  nonlocal  terms  In  Eg.  (3.34)  and  obtain  the  local 
theory  for  the  stream  function,  i.e. 

^£=1RB(x)  (3.46) 

dx^ 

The  solution  for  the  one-dimensional  case  can  be  found  by  direct  Integration 
of  a  given  B(x).  Comparison  of  the  low  Reynolds  number  solution  for  the 
finite  element  method  (FEM)  and  direct  quadrature  Is  shown  In  Figure  3.4 
for  a  pair  of  current  ^laments  centered  above  a  long  conducting  plate.  The 
results  show  excellent  agreement. 

b)  Comparison  with  the  Image  Method.  When  the  frequency  Is  high,  i.e., 
R»l,  the  nonlocal  effects  act  to  prevent  the  magnetic  field  from  penetrat¬ 
ing  the  plate.  In  this  limit  the  solution  can  be  approximated  by  consider¬ 
ing  an  Image  coll  below  the  plate  which  serves  to  cancel  out  the  normal 
component  of  the  total  magnetic  field  (Figure  3.5).  The  results  of  the 
one-dimensional  FEM  with  the  Image  coll  calculation  are  shown  In  Figure  3.6 
and  again  show  excellent  agreement. 

c)  Nonlocal  vs.  Local  Theory  for  Eddy  Currents.  A  comparison  of  the 
local  and  nonlocal  theories  for  eddy  current  distribution  Is  shown  In 
Figure  3.7  for  different  R  or  frequencies.  It  Is  clear  that  even  though 
the  skin  depth  may  be  several  times  the  thickness  (R  *  0.1,  6  *  4.5h),  the 
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FEM  and  Direct  Quadrature  Solutions 
■  a  pair  of  current  filaments  centered 
conducting  plate  (R  =  0.01). 
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Figure  3.5  The  image  method  for  high  R. 
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Figure  3.6  Comparison  of  FEM  and  image  method  solutions 
for  high  R  for  a  pair  of  current  filaments 
centered  above  a  long  conducting  plate  (R  =  5). 
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nonlocal  distribution  differs  significantly  from  the  local  theory.  The 
Importance  of  the  nonlocal  effects  Is  not  recoqnlzed  by  some  eddy  current 
calculations  which  neglect  the  reaction  field  of  the  Induced  current 
(Ref.  80). 

A  more  fundamental  difference  between  the  nonlocal  and  local  effects 
concern  the  magnetic  pressure.  When  the  nonlocal  or  self  field  of  the 
Induced  currents  Is  neqlected,  the  average  force  on  the  plate  for  a  sinus¬ 
oidal  current  Is  zero.  The  nonlocal  or  self  field  effects  are  necessary 
to  obtain  a  time-averaged  force  on  the  plate  between  the  sources  of  the 
external  field  (such  as  colls  or  current  filaments)  and  the  currents  In 
the  plate. 

The  phase  difference  between  the  drivlnq  maqnetlc  field  and  the 
Induced  eddy  current  Is  90  degrees  In  the  local  theory,  as  Is  evidenced 
by  Eq.  (3.46).  The  phase  difference  Is  zero  In  the  Imaoe  solution.  The 
Induced  current  densities  In  Figures  (3.4)  and  (3.6)  therefore  have  the 
same  phase  across  the  width  of  the  plate.  For  the  Intermediate  case,  the 
eddy  current  density  has  different  phase  angles  at  different  locations  of 
the  plate.  Figure  3.7  shows  the  modulus,  or  the  maximum  value,  of  the 
current  density  across  the  plate. 

d)  Comparison  of  FE  Calculations  with  Experimental  Measurements. 

An  Important  feature  of  this  research  program  has  been  the  experimental 
verification  of  the  numerical  calculations.  To  check  the  calculations, 
measurements  of  one  of  several  physical  quantities  along  the  surface  of 
the  plate  must  be  made  such  as  electric  or  magnetic  fields.  Induced 
temperature,  or  magnetic  pressure.  Measurement  of  electric  or  magnetic 
fields  associated  with  the  Induced  currents  In  the  plate  involve  taklnq 
data  of  a  sufficient  number  of  points  to  map  out  the  current  distribution. 
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It  has  the  advantage  that  the  time  hibtory  of  the  electromaanetlc  field 
variable*  at  buch  points  may  be  more  easily  determined.  Such  a  technique 
will  be  used  in  the  experimental  verification  of  the  results  of  the 
transient  analysis.  To  verify  the  current  distribution,  an  Infrared 
scanning  technique  which  senses  the  Incremental  temperature  distribution  Is 
used  (Ref. 102).  This  technique  is  based  on  the  fact  that  for  small  times 
after  the  Induced  currents  are  generated,  the  heat  conduction  can  be 
neglected  and  the  Induced  temperature  Is  proportional  to  the  square  of 
the  current  density.  The  heat  equation  Is  given  by 


- kV2J  +  c  ^  (3.47) 

In  which  <  Is  the  thermal  conductivity,  and  c  is  the  heat  capacity  of  the 
plate.  Under  appropriate  conditions  on  the  gradient  of  temperature  and 
time,  one  may  neglect  the  first  term  and  write 


T 


ft 


J2dt 


(3.48) 


o 

To  simulate  a  one-dimensional  problem,  a  rectangular  multi  turn 

induction  coll  is  placed  parallel  to  a  flat  stainless  steel  plate  and 

induced  temperatures  are  measured  across  the  plate  under  the  middle  of  the 

coil.  Figure  3.8.  The  two-dimensional  Infrared  scanner  used  is  described 

2 

In  Ref.  102,  An  Infrared  thermogram  of  the  temperature  or  J  pattern  in 
the  plate  due  to  a  pulsed  rectangular  coil  Is  shown  in  Figure  3.9.  The 
Induction  coil  has  damped  sinusoidal  currents.  Comparison  of  the  qualita¬ 
tive  behavior  of  the  measured  temperature  and  the  calculated  temperature 
distribution  alonq  the  plate  uslnq  an  effective  frequency  Is  shown  in 
Figures  3.10  and  3.11  for  two  different  frequencies.  Both  calculated  and 


Figure  3.8  Method  of  detecting  eddy  current  flow  patterns 
in  thin  structures. 


Nonlocal, 
steady  state 


Comparison  of  f"C  solution  and  measured  temperature 
for  a  pair  of  current  filaments  nearer  one  edge  of 
a  long  conducting  plate  (R  =  0.074). 


measured  data  have  been  normalized.  Excellent  qualitative  aoreement  can 
be  seen  for  the  coll  centered  above  the  plate,  Figure  3.10,  and  for  the 
coll  near  the  edge  of  the  plate.  Figure  3.11. 

Quantitative  aqreement  Is  difficult  because  the  calculations  are  for 
a  harmonic  excitation  field  and  the  experiment  uses  pulsed  current  In  the 
coll.  This  problem  Is  solved  by  calculating  the  induced  current  density, 
as  a  function  of  frequency,  and  uslna  a  fast  Fourier  transform  (FFT)  to 
calculate  the  Induced  current  as  a  function  of  time,  as  shown  In  Fiaure 
3.12.  This  data  Is  then  Integrated  using  Eg.  (3.48)  above  to  calculate 
the  temperature  as  a  function  of  time  at  the  point  of  maximum  temperature 
In  the  plate.  Using  the  maximum  temperature  In  time  alona  with  the  calcu¬ 
lated  distribution  as  In  Figure  3.12,  a  quantitative  comoarlson  of  calcu¬ 
lated  and  measured  temperature  were  made  as  shown  In  Fiaure  3. IP.  The 
measured  value:  are  about  20%  below  the  calculated  values.  A  comparison 
of  the  FFT  result  and  a  subsequent  transient  analysis  result  for  the  same 
problem  Is  given  In  sub-subsection  3.3.2.b.  In  light  of  this  later 
comoarlson,  the  difference  between  the  experimental  result  and  the  FFT 
analysis  appears  to  arise  partly  from  the  rough  calculation  In  the  FFT, 
l.e.,  the  cutoff  frequency  used  may  be  too  low  and  the  number  of  terms 
used  In  the  finite  sum  approximation  of  the  Fourier  Intearal  may  be  In¬ 
sufficient.  The  crudeness  of  the  FFT  result  Is  also  evidenced  by  the  late 
arrival  of  the  peak  value  of  the  induced  current  In  Fiaure  3.12.  In  the 
FFT  analysis  of  the  present  problem  a  half-sine  pulse  of  the  drlvlna  current 
Is  used,  while  In  the  experiment  the  real  driving  current  has  a  damped 
half-sine  variation.  The  major  cause  of  the  difference  between  the  calcu¬ 
lated  peak  temperature  and  the  measured  one  Is  probably  heat  conduction 
the  effect  of  which  Is  not  Included  In  the  calculation.  The  Infrared 


Figure  3.12  Transient  exciting  and  induced 
currents  as  functions  of  time. 


picture  was  taken  approximately  160  ms  after  the  pulse  was  applied.  This 
may  explain  the  loss  of  the  sharp  peak  of  the  temperature  curve  under  the 
wire.  In  sum,  the  FFT  result  obtained  for  this  problem  Is  satisfactory 
considering  the  approximations  Involved. 

e)  FE  and  FTT  Methods  for  Pulsed  Maonetlc  Fields.  As  discussed 
above,  the  distribution  of  current  In  the  plate,  as  well  as  Induced 
pressure  and  temperatures  are  for  harmonic  excitation.  For  pulsed  or 
transient  excitation  fields  jl0(t),  one  may  decompose  the  field  Into  Its 
spectral  components 


V<“> '  k 


B(t)  e1uitdt 


(3.49) 


If  the  Induced  current  Is  calculated  as  a  function  of  frequency  J(u>),  for 
an  excitation  Boe"1a)t,  the  time  variation  of  current  can  be  found  from  the 
Integral 

m 

J(t)  J(<o)  f(u)e"i“td(D  (3.50) 

The  function  J(uj)  Is  found  by  calculatlna  the  Induced  current  for 
various  frequencies  uslno  the  nonlocal  theory  as  shown  In  Flqure  3.13. 

A  polynomial  Is  then  fitted  to  these  points  over  the  slonlflcant  frequency 
domain  of  f((u).  The  Fourier  Inteoral  Is  then  approximated  by  a  finite 
sum  and  the  summation  carried  out  uslno  a  fast  Fourier  transform  aloorlthm. 
The  resulting  time  history  of  current  In  the  plate  Is  shown  In  Flaure  3.12. 
From  this  history  the  maonetlc  force  and  Induced  temperature  histories  In 
the  plate  can  be  calculated.  A  detailed  discussion  of  the  FFT  method  In 
the  transient  solution  of  the  diffusion  equation  has  been  given  In  Chapter 
8  of  Reference  54. 


Figure  3.13  Spectrum  for  current  at  x  =  -0.62A  for 
local  and  nonlocal  solutions. 


f)  Effect  of  Edges  on  Induced  Eddy  Currents.  One  of  the  Immediate 
results  of  Infrared  scan  experiments  on  Induced  eddy  currents  In  plates 

Is  the  observation  that  hlah  current  densities  will  occur  when  the  excita¬ 
tion  field  Is  concentrated  near  the  edges  of  the  plates.  This  can  be  seen 
In  Figure  3.14a,b,  where  a  coll  Is  moved  closer  to  the  edge  of  a  rect¬ 
angular  plate.  The  results  of  the  one-dimensional  program  also  show  a 
distinct  edge  effect  for  both  the  paired  current  filament  Induction, 

Figure  3.4,  and  the  uniform  field  Induction,  Floure  3.15. 

g)  Pressure  Distributions  Due  to  Tilted  Induction  Coils.  Both 
magnetic  forming  devices  and  magnetic  levitation  devices  use  current  carry¬ 
ing  colls  near  plate-like  conductors.  The  effect  of  tlltlno  the  induction 
coll  has  been  calculated  as  shown  In  Figures  3.16  and  3.17.  The  Increased 
magnetic  pressure  under  the  filament  close  to  the  plate  produces  a  moment 
on  the  plate  and  a  restoring  moment  on  the  coll.  The  effect  of  lateral 
movement  of  the  coll  can  also  be  seen  In  Figure  3.17,  and  the  effect  of 
lateral  movement  on  the  restorlna  moment  can  be  calculated. 

h)  Effect  of  Reduced  Matrix  Band  on  Nonlocal  Solutions.  One 

manifestation  of  the  FE  Implementation  of  the  nonlocal  theory  Is  that  the 

algebraic  FE  eguatlons  become  full  rather  than  banded.  Nevertheless,  the 

nonlocal  effects  which  cause  this  loss  of  bandlno  are  proportional  to 
2 

1/r  In  which  r  Is  the  distance  from  a  self-field  source  point  on  the 
sheet  conductor  to  the  field  point  In  auestlon,  Eo.  (3.34).  One  method 
of  restricting  the  nonlocal  effect  at  any  field  point  on  the  sheet  would 
be  to  exclude  source  points  at  distances  from  the  field  point  greater 
than  some  cut-off  value  of  r.  For  a  uniformly  spaced  finite  element  mesh, 
this  would  yield  a  banded  matrix,  but  the  size  of  the  band  would  now  be 
determined  by  the  cut-off  distance.  The  choice  of  an  appropriate  cut-off 


Figure  3.14  Photographs  showing  effects  of  edges 
on  induced  eddy  currents. 


Figure  3.16  Magnetic  pressures  due  to  a  tilted  pair  of  current 
filaments  centered  above  a  long  conducting  plate 


Magnetic  pressures  due  to  a  tilted  pair  of  current 
filaments  nearer  one  edge  of  a  long  conducting 
plate  (R  =  0.07,  tilt  =  10  deg). 


distance  will  be  dependent  upon  the  particular  problem.  The  solution  to 
the  restricted  nonlocal  problem  with  an  appropriate  cutoff  distance  would 
be  neqllglbly  different  from  the  full  nonlocal  solution,  and  considerable 
savlnqs  In  computation  would  be  achieved  due  to  a  reduction  In  the  non¬ 
local  Integrations  required,  Eq.  (3.34). 

As  a  preliminary  test  of  this  strategy,  the  expedient  of  merely 
reducing  the  bandwidth  of  the  FE  equations  was  utilized  by  neqlectlnn  all 
terms  outside  of  selected  bandwldths.  For  the  nonuniform  mesh  necessary 
for  the  eddy  currents  Induced  by  a  pair  of  external  current  filaments, 
this  procedure  Is  not  equivalent  to  the  selection  of  a  cutoff  distance. 

(A  finer  mesh  Is  required  under  each  filament).  Nevertheless,  as  Illus¬ 
trated  In  Flqure  3.18,  the  reduction  In  bandwidth  provides  a  variation 
In  the  solution  that  approaches  the  full  nonlocal  solution.  These  results 
provide  evidence  that  reduced  nonlocal  solutions  will  be  efficient  and 
useful  In  some  situations.  Further  work  will  be  needed  to  develop  criteria 
for  the  selection  of  appropriate  cutoff  distances  and  to  Implement  It  in 
the  computer  prooram.  For  the  efficient  solution  of  laroe,  dense  matrix 
systems,  Lachat  and  Watson  (Ref.  81)  presented  a  similar  ad  hoc  approach 
for  the  numerical  treatment  of  boundary  inteoral  equations  In  three- 
dimensional  elastostatlcs.  Jenq  and  Wexler  (Ref.  82)  described  a  successive 
element  Iterative  scheme  In  their  studies  of  three-dimensional  static  field 
problems  using  a  boundary  Integral  equation  formulation. 

1)  Dimensional  Analysis  of  the  Size  Effect  of  the  Plate.  The 

numerical  results  obtained  from  proqram  EDDY1  presented  so  far  have  been 

based  on  a  single  nondlmenslonal  parameter:  the  maonetic  Reynolds  number  R. 

R  Is  defined  as  R  *  wyoh  and  Is  related  to  the  skin  depth  6  through 
2 

R  a  2(h/s)  ,  as  discussed  in  connection  with  Eq.  (3.34).  For  infinite 

i 


Figure  3.18  Effect  of  reduced  matrix  band  on  nonlocal  solutions  for  a 
pair  of  current  filaments  centered  above  a  long  conducting 
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plates,  except  for  the  geometric  parameters  characterizing  the  positions 
of  the  excitinq  wires,  this  is  the  only  possible  nondimensional  parameter 
characterizing  the  plate.  For  plates  with  finite  width  t,  the  size  effect 
of  the  width  (span)  comes  in  an  another  nondimensional  parameter  charac¬ 
terizing  the  effects  of  both  the  thickness  and  the  width  of  the  plate  Is 
possible.  The  new  nondimensional  parameter  Is  designated  by  Rw  and  is 
defined  and  related  to  R  as 

R  *  uuoht  =  2h t/6^  =  Rt/h  (3.51) 

w 

The  role  of  Rw  is  illustrated  through  the  one-dimensional  plates  with 
varyinn  width/thickness  ratio  in  a  uniform,  sinusoidally  varyino  maqnetic 
field.  The  analyses  with  EDDY1  have  been  performed  by  Hara  (Ref.  83). 

Fioure  3.19  shows  the  maximum  induced  current  at  the  edoes  of  the 
plate  versus  the  Reynolds  number  R  usina  t/h  as  a  parameter.  Figure  3.20 
shows  that  the  different  curves  can  be  collapsed  to  a  s inale  curve  by 
using  the  new  nondimensional  parameter  Rw. 

Figure  3.21  shows  the  sum  of  the  moduli  of  the  Induced  currents 
through  the  width  of  the  plate  versus  R  for  different  t/h  values.  Again, 
the  different  curves  can  be  mapped  to  one  by  usina  the  new  parameter  R 

n 

as  shown  in  Figure  3.22. 

To  determine  the  total  torque  and  force  on  a  conductor,  one  needs 
to  know  the  real  current  distribution.  Both  the  real  and  the  imaginary 
parts  of  the  induced  currents  therefore  need  to  be  considered  (Ref.  33). 
Figure  3.23  shows  the  real  and  Imaginary  parts  of  the  total  Induced 
current  as  functions  of  the  Reynolds  number  R  for  two  different  t/h 
values.  The  real  part  a'  is  the  total  current  across  the  width  In  phase 
with  the  driving  magnetic  field,  the  imaginary  part  a"  is  the  total 
current  90  degrees  out-of-phase  with  the  driving  field.  These  curves 


/aI/B  at  edge  of  plate 


Figure  3.19  Maximum  normalized  current  density  at  the  edge 
of  a  plate  subject  to  a  uniform  transverse  field 
as  a  function  of  Reynolds  number  R  =  2(h/6)z. 


Total  induced  current  at  ait  =  2 w 


are  plotted  In  Figure  3.24  uslnn  the  new  parameter  R  .  Aqaln,  qood 
correlation  Is  obtained  althouph  some  differences  occur  at  high  Reynolds 
number  cases. 

3.3.2.b  Transient  Analysis  by  Propram  EDDY1T 

The  transient  analysis  propram  EDDY1T  has  been  verified  analytically 
and  experimentally  In  various  cases.  The  followlno  Is  a  summary  of  some 
of  the  results: 

a)  Comparison  of  the  local  solutions  with  the  analytic  results. 

b)  Comparison  of  the  steady  state  nonlocal  solution  with  the 
ED0Y1  results. 

c)  Comparison  of  the  transient  solution  with  the  FFT  and 
Infrared  experimental  results  of  Ref.  102. 

d)  Comparison  of  the  transient  solution  and  experimental 
measurements  of  eddy  currents  obtained  with  a  search  coll  (Ref.  101). 

All  the  cases  analyzed  use  the  physical  parameters  of  the  problem  In  the 
rationalized  MKS  system.  Each  of  these  problems  Is  described  In  more 
detail  below. 

a)  Comparison  of  the  local  Solutions  with  the  Analytic  Results.  The 
local  solutions  of  the  two-wire  problem  studied  previously  with  the  program 
EDDYl(Sect1on  3. 3. 2. a;  Figure  3.10)are  calculated  and  presented  In  Floure 
3.25.  Since  the  reaction  field  Is  neglected,  the  solution  may  be  obtained 
by  direct  Integration  and  Is  the  same  for  each  cycle  of  excitation.  A 
stainless  steel  plate  Is  used  In  this  problem.  The  FFM  solutions  obtained 
with  two  different  time-step  sizes  do  not  differ  significantly  from  the 
analytic  results. 


Figure  3.24  Averaged  curr 


Figure  3.25  Comparison  of  theoretical  and  finite  element 
method  (FEM)  predictions  of  local  transient 
eddy  current  effects  for  the  long  plate  in 
Figure  3.12  subject  to  a  sinusoidally  vary¬ 
ing  magnetic  field . 
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b)  Comparison  of  the  Steady  State  Nonlocal  Solutions  with  the  EDDYl 
Results.  The  nonlocal  solution  of  the  same  problem  with  a  continual 
sinusoidal  excltlnq  current  Is  calculated  by  proaram  EDDY1T.  Fiaure  3.26 
shows  that  the  Induced  current  at  point  6  becomes  nearly  steady  state 
after  two  cycles  of  excitation.  The  maanetlc  Reynolds  number  In  this 
problem  Is  R  ■  0.079.  The  modulus  and  phase  anole  of  the  induced  current 
when  It  becomes  steady  state  are  the  same  as  predicted  In  Fiaure  3.13. 

In  Fiaure  3.26b,  the  maanetlc  oushinq  force  is  shown,  and  the  fact  that 
Its  averaqe  value  Is  nonzero  Is  clearly  evident. 

The  size  of  the  time  steD  used  In  this  oroblem  Is  At  =  0.041  ms,  or 
one  twentieth  of  the  period  of  the  drlvlno  current.  Flqure  3.27a  shows 
the  eddy  current  distribution  at  the  53-rd  time  step  when  the  Induced 
current  at  point  G  is  almost  at  Its  maximum.  Floure  3.27b  shows  the 
current  distribution  at  the  57-th  time  step  when  the  Induced  current  at 
point  G  is  nearly  zero.  Fiaure  3.27c  shows  the  current  distribution  at 
the  60-th  time  step  (<ut  *  6ir)  when  the  drlvlnq  current  in  the  wire  is 
zero.  Fiqure  3.27d  shows  the  current  distribution  at  the  65-th  time  step 
{wt  ■  6»yw)  when  the  drivina  current  in  the  wire  reaches  Its  maximum  value. 
The  correspondlnp  maanetlc  force  distributions  at  these  particular  time 
steps  are  shown  In  Fiaure  3.28. 

c)  Comparison  of  the  Transient  Solution  with  the  FFT  and  Infrared 
Experimental  Results.  The  transient  eddy  current  problem  discussed  in 
Flqures  3.10,  3.12,  and  3.13  Is  re-analyzed  uslna  proaram  EDDYlT.  The 
results  are  shown  In  Flqures  3.29  and  3.30.  The  Induced  current  density 
calculated  at  point  G  underneath  the  wire  Is  presented  In  Fiaure  3.2R 
and  compared  with  the  FFT  result  from  Floure  3.12.  It  Is  seen  that  the 
two  results  have  nearly  the  same  maonltude  but  the  one  calculated  by  EDPY1T 


Force  (xIO  N/m)  Current  density  (x  l03omp/m) 


induced  tempercture,  ru 


shows  more  accurate  time  characteristics.  This  Is  evidenced  by  the  late 
arrival  of  the  second  peak  of  the  FFT  result  which  should  occur  physically 
at  the  time  when  the  drlvlnq  current  becomes  zero.  The  correct  arrival 
time  of  the  second  peak  has  been  predicted  by  the  EDPY1T  result. 

The  Induced  temperature  rise  Is  shown  In  Floure  3.30  and  compared 
with  the  Infrared  experimental  results.  The  Infrared  picture  Is  taken 
approximately  160  ms  after  the  pulse  Is  applied,  as  mentioned  In  connec¬ 
tion  with  Flqure  3.10.  The  prediction  of  the  FFT  result  for  the  tempera¬ 
ture  Is  also  shown  for  comparison. 

d)  Comparison  of  the  Transient  Solution  and  Experimental  Measure¬ 
ments  of  Eddy  Current  Obtained  with  a  Search  Coll.  A  search  coll 
technique  has  been  used  to  measure  the  eddy  current  Induced  on  an 
aluminum  plate.  The  arranaement  of  the  experiment  is  described  schemetl- 
cally  In  Flqure  3.31.  To  simulate  a  one-dimensional  plate,  the  measure¬ 
ments  were  taken  alono  the  central  line  of  a  rectanoular  plate  with 
lenqth/width  ratio  approximately  eoual  to  4.2  (Ref.  101). 

To  measure  the  real  Induced  current,  the  search  coll  readlnos  at 
points  A  and  B  were  first  taken  usino  a  plastic  dummy  plate.  The  measure¬ 
ments  were  then  made  at  points  A  and  R  for  the  conduction  plate.  The 
induced  current  densities  are  obtained  by  subtractlna  the  readlnqs  obtain¬ 
ed  for  the  plastic  plate  from  the  readlnqs  for  the  real  plate.  The  search 
coll  readlnos  obtained  at  point  A  underneath  the  wire  are  shown  in  Floure 
3.32. 

The  Induced  current  densities  at  points  A  and  P  so  determined  are 
compared  with  the  numerical  results  obtained  by  usino  prooram  EDDY1T  In 
Flqure  3.33.  The  qualitative  aqreement  of  the  results  are  aood  at  both 
points.  Quantitatively  some  differences  exist.  These  quantitative 


Oscilloscope 


Search  coil 


d  =  32  mm,  h  =  1 1  mm 


t 

Figure  3.31  Arrangement  of  the  search  coil  experiment  for 

a  long  plate,  (a)  Isometric  view,  (b)  Plan  view. 
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Figure  3.33  Numerical  and  search  coil  experimental  results  of  the 
induced  current  density  for  the  plate  in  Figure  3.31: 
(a)  at  point  B,  (b)  under  the  wire  at  point  A. 
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differences  probably  arise  laroely  from  the  approximations  Involved  In  the 
experiment.  Namely,  each  of  the  two  search  colls  used  has  a  finite  lenqth 
of  about  10  mm  and  finite  cross  sectional  dimension  of  about  10  mm.  The 
dimension  of  the  excltlno  wire  Itself  Is  about  ?  mm.  The  measurements 
obtained  are  therefore  averaged  values  over  the  10  mm  lenoth  of  the  search 
coll.  One  example  of  this  averaging  effect  in  the  experimental  results 
can  be  seen  from  Flqure  3.33.  The  two  experimental  curves  have  nearly 
the  same  peak  value  at  points  A  and  B,  while  In  the  Infrared  results  and 
in  all  the  numerical  calculations.  Figure  3.34,  the  Induced  current  under 
the  wire  Is  definitely  higher. 

3.4  Eddy  Current  on  Finite  Plates 

The  finite  element  Oalerkin  formulation  and  the  Fortran  code  EDDY2  for 
two-dimensional  problems  of  steady  state,  harmonic  currents  In  flat  plates 
are  presented  In  Section  3.4.1.  The  treatment  of  thp  nonlocal  Inteoral 
term  Is  described  In  some  detail.  Some  numerical  results  and  their  verifi¬ 
cations  are  presented  In  Section  3.4.?. 


3.4.1  Finite  Element  Formulation  for  Steady  State  Problems 

The  two-dimensional  eddy  current  equation  Is  obtained  from  Fq. 
(2.103)  by  dropping  all  the  coupling  terms.  For  steady  state,  harmonic 
currents  In  a  flat  plate,  the  eouatlon  may  be  nondlmensional Ized  Into  the 
following  form. 
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In  which  the  coordinates  are  nondlmensional Ized  with  respect  to  half  the 
thickness  (■>}■);  1 1>  =  {■^■He^u’t;  B°  *  (y)ee^10^;  the  maanetlc  Reynolds  number 


Is  R 


,  which  is  related  to  skin  depth  £  throuoh  R 
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B  Is  the  reference  maonetlc  field. 

The  finite  element  GalerMn  method  Is  used  to  solve  Fo.  (3.51).  q  Is 
approximated  Globally  and  locally  by 
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in  which  G  is  the  total  number  of  nodal  points,  F  denotes  the  Fth  element, 
are  the  quadratic  olobal  Interpolation  functions  oenerated  from  the 
local  element  shape  functions  N^.  Six-node  trianoular  elements  are  used 
here.  The  local  element  shape  functions  are  all  nuadratlc  In  this  case. 
The  element  alqebralc  equations  are 
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The  weiqhtlno  function  W;j  Is  alven  by 


(3.53) 


(3.54a) 
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A  change  of  the  sequence  of  Integrations  In  the  Galerkin  formulation  of 
the  Integral  term  In  Eq.  (3.51)  has  been  applied  to  obtain  the  weighting 
function  W*j.  Because  of  the  numerical  difficulty  associated  with  the 
sharp  variation  of  the  kernal  function  in  Eq.  (3.55),  the  welqhtino 
function  Is  calculated  analytically  within  each  element  and  numerically 
outside  the  element.  Eq.  (3.54c)  is  then  Inteorated  entirely  by  numerical 
quadrature. 

The  formulation  in  Eqs.  (3.54c)  and  (3,55)  is  in  terms  of  the  qlobal 
coordinates  through  the  global  position  of  the  source  point  U,n).  To 
simplify  the  Integration,  each  element  Is  first  mapped  onto  a  standard 
triangle  through  a  rigid  body  translation  and  rotation.  A  linear  change 
of  scales  of  the  axes  Is  then  performed  so  that  the  standard  triangle  Is 
independent  of  the  coordinates  of  the  nodes  of  the  element  being  mapped, 
and  the  coordinates  of  the  source  point  appear  only  as  parameters  in  the 
analytical  Integration.  The  same  form  of  the  Integrated  result  may  then 
be  used  for  all  the  elements  In  the  finite  element  mesh. 

Since  N*j  Is  quadratic,  six  basic  Integrations  with  numerators  1,  x, 

2  2 

y,  x  ,  xy,  and  y  in  Eq.  (3.55)  need  to  be  performed.  These  Integrated 
results  are  used  to  form  the  weighting  function.  The  six  Integrations, 
the  weighting  function  w|j,  and  the  element  matrices  [KE]  and  [PE]  are 

detailed  In  Appendix  B. 

The  resulting  global  matrix  Is  again  complex,  nonsyrrmetrlc,  and  fully 
populated.  Local,  nonlocal,  and  Image  solutions  can  be  calculated  .lust  as 
In  the  one-dimensional  case.  Eddy  currents  can  be  calculated  through 
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numerical  differentiation,  and  maanetlc  pressure  and  temperature  evaluated 
after  the  stream  function  is  solved. 

A  Fortran  program  EDDY2  has  been  developed  based  on  the  formulation  of 
Eqs.  (3.53)  -  (3.55).  It  calculates  the  local  and  nonlocal  solutions  of 
stream  function,  eddy  current,  temperature,  and  pressure.  As  of  this 
writing,  the  image  solution  and  two-dimensional  graphic  output  capabilities 
remain  to  be  Implemented.  Uniform  magnetic  field  and  fields  due  to  any 
number  of  magnetic  dipoles  can  be  handled.  Magnetic  fields  aenerated  from 
other  types  of  colls  of  Interest  can  be  added  readily.  The  capability  of 
analyzing  for  multiple  frequencies  also  remains  to  be  implemented. 

The  geometry  of  the  plate  and  the  description  of  the  external  magnetic 
field  are  the  two  basic  forms  of  data  needed  by  Drooram  EDDY2.  The  total 
numbers  of  nodal  points,  load  cases,  and  elements  need  to  be  specified. 
Coordinate  and  boundary  condition  must  be  given  for  each  Input  node. 
Intermediate  nodes  may  be  generated  for  any  groups  of  nodal  points  that 
are  uniformly  spaced.  Element  information  may  also  be  generated.  Although 
only  six-node  triangles  are  Included  In  the  present  version  of  E0DY2,  the 
program  has  been  structured  so  that  other  types  of  elements  may  be  added. 
Element  group  Information  and  the  master  card  for  each  element  group 
therefore  need  to  be  Input  too. 

The  program  allows  for  different  orders  of  numerical  integrations. 

Six-  and  seven-point  formula  are  now  provided.  The  order  may  be  specified 
on  the  master  element  qroup  card. 

Presently  two  types  of  magnetic  fields  may  be  analyzed:  the  uniform 
field  and  the  magnetic  field  due  to  any  number  of  dipoles.  For  uniform 
field  no  other  Information  Is  needed.  For  dipole  field  the  positions  and 
the  relative  strengths  and  senses  of  each  magnetic  dipole  must  be  given 


subsequently. 

Values  of  the  stream  function,  eddy  current,  temperature  Induced  in 
a  half  cycle  of  the  current,  and  time-averaged  magnetic  pressure  are 
produced  as  output.  The  stream  function  is  calculated  at  the  nodes  of 
the  finite  element  mesh.  Current,  temperature,  and  pressure  are  evaluated 
at  the  centroid  of  each  element.  The  stream  function  and  current  are 
calculated  In  complex  form.  The  modulus  and  phase  angle  of  the  current 
are  evaluated  In  the  Interest  of  spectral  analysis  for  the  calculation  of 
transient  currents. 

3.4.2  Numerical  Results  and  Experimental  Verifications  for  Finite  Plates 
The  two-dimensional  finite  element  code  EDDY2  has  been  verified  for 
a  limited  number  of  problems  involving  uniform  applied  fields.  The 
numerical  results  from  EDDY2  have  been  checked  with  both  analytic  solutions 
and  numerical  calculations  obtained  from  other  methods.  Infrared  measure¬ 
ments  are  also  used  to  verify  the  results  qualitatively.  More  extensive 
verifications  are  needed  to  check  the  full  capabilities,  but  comparison 
results  for  fields  other  than  uniform  are  scarce  except  for  low  Reynolds 
number  (local)  cases. 

The  following  objectives  have  been  met  with  the  program  EDDY2  up  to 
this  writing. 

a)  Comparison  of  low  Reynolds  number  results  with  the  analytic 
solutions, 

b)  Comparison  of  low  and  Intermediate  Reynolds  number  results 
for  long  plate  with  the  numerical  results  from  program  EDDY1. 

c)  Comparison  of  low  Reynolds  number  results  for  a  square 
plate  with  edge  crack  with  the  BIEM  numerical  results  and 
Infrared  measurements  (Ref.  102). 


d)  Dimensional  analysis  of  the  size  effect  of  the  width  of  a 


square  plate. 

Each  of  these  is  described  separately  below. 

a)  Comparison  of  Low  Reynolds  Number  Results  with  the  Analytic 
Solutions.  In  the  low  Reynolds  number  limit  one  may  drop  the  nonlocal 
terms  In  Eq.  (3.51)  and  obtain  the  local  theory  for  the  stream  function 
in  the  two-dimensional  case, 

v2*  =  12irRe(x,y)  (3.56) 

The  solution  of  the  above  Poisson  eouation  is  a  common  practice  in  many 
physical  problems.  For  the  uniform  normal  external  magnetic  field,  the 
problem  is  analogous  to  the  torsion  of  a  shaft  in  the  theory  of  elasticity, 
for  example.  For  a  square  and  rectangular  plate  In  the  low  Reynolds 
number  limit, the  EDDY2  solution  has  been  checked  with  the  series  solution 
for  torsion  of  rectangular  shafts.  The  agreement  between  the  stream 
function  and  current  density  in  one  case  and  the  stress  function  and  shear 
stress  in  the  other  is  excellent. 

b)  Comparison  of  Low  and  Intermediate  Reynolds  Number  Results  for 
Long  Plate  with  the  Numerical  Results  from  Program  EDDYl.  For  a  long 
rectangular  plate  excited  by  a  harmonic  uniform  field  with  the  Reynolds 
number  R  =  0.0012,EDDY2  gives  the  stream  function  contour  shown  in  Figure 
3.35.  Across  the  middle  section  of  the  Iona  plate,  the  problem  is 
essentially  one-dimensional  and  may  be  approximated  by  an  analysis  with 
EDDYl.  The  close  correlation  between  the  two  calculated  induced  currents 
across  that  section  from  EDDY2  and  EDDYl  are  shown  in  Figure  3.36. 

For  the  intermediate  Reynolds  number  case  (R  =  0.2  in  EDDYl),  the 
comparison  between  the  stream  functions  calculated  from  EDDYl  and  EDDY2 
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Figure  3.36  Comparison  of  one-  and  two-dimensional  solutions 
of  currents  induced  across  the  middle  of  a  long 
rectangular  plate  by  a  harmonic  uniform  field 
(R  =  0.03  in  EDDY1 ) . 


are  shown  in  Fiaure  3.37.  The  EDDY!  results  are  obtained  with  ten  linear 
elements  over  the  half-width  of  the  plate,  while  the  EDDY2  results  are 
from  a  sinqle  quadratic  element  over  the  same  distance, the  aoreement 
between  the  two  should  be  therefore  reasonable  considerino  the  difference 
in  the  discretization. 

c)  Comparison  of  Low  Reynolds  Number  Results  for  a  Souare  Plate 

with  Edae  Crack  with  the  BIEM  Numerical  Results  and  Infrared  Measurements. 

The  flow  of  induced  current  around  an  edqe  crack  In  a  souare  plate  is 

calculated  using  the  two-dimensional  code  in  the  laroe  skin-depth  limit 

for  a  few  notch  widths.  A  uniform  applied  mannetic  field  is  assumed  in 

these  calculations.  The  stream  function  contours  for  one  case  are  shown 

in  Figure  3.38.  These  contours  show  the  flow  of  current  around  the  crack 

tip.  Contours  of  constant  temperature  are  also  shown  in  Figure  3.39. 

The  infrared  measurements  of  the  induced  ohmic  heating  on  the  same  olate 

but  induced  by  a  wire-source  field  are  shown  in  Fioure  3.40.  The  tempera- 
2 

ture  or  J  hot  spot  near  the  edqe  of  the  crack  or  slit  is  clearly  shown. 

The  same  hot  soot  is  predicted  in  the  FE  results,  Figure  3.30,  f0r  the 
uniform  maanetic  field. 

Quantitative  comparison  of  the  numerical  and  experimental  results  is 
difficult  in  this  problem  because  of  the  difficulty  in  creatino  a  uniform 
field  within  a  sufficiently  larae  region  with  the  present  experimental 
facilities.  However,  the  FE  results  have  been  compared  with  the  numerical 
results  obtained  from  a  boundary  integral  equation  (R IE )  code  developed  by 
Mukherjee  (Ref.  84),  This  code  has  oriqinally  developed  to  solve  for 
stresses  around  a  crack.  The  comparison  of  the  stream  functions  calculated 
is  shown  In  Figure.  3.41.  Close  aoreement  is  obtained  between  the  two 
numerical  results. 
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Figure  3.37  Eddy  current  '•.i.rcan  function  on  the  nid-lcngth 
crocs-nect ion  of  a  long  rectangular  plate 
( length :width  =  4:1,  R  =  0.2). 


Figure  3.38  Stream  function  contours  for  a  notched  plate 
excited  by  a  harmonic  uniform  field  (notch- 
width  =  2h,  R  =  0.001). 
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Figure  3.39  Isotherms  for  a  notched  plate  excited  by  a 
harmonic  uniform  field  (notchwidth  =  2h, 

R  =  0.001). 


- ~ 
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Figure  3.40  Infrared  thermogram  showing  hot  spots  due  to  eddy  current 
flow:  (upper)  plate  without  crack,  (lower  left)  hot  spot 
due  to  flow  around  a  crack,  (lower  right)  color  quantized 
hot  spot  at  top  of  crack  (magnified). 


Comparison  of  eddy  current  stream  lines  for  a  notched 
square  plate  in  a  uniform  field  computed  by  boundary 
integral  equation  method  (BIEM)  and  finite  element 
method  (FEM).  (Notchwidth  =  2h,  notch  length  =  40h, 
plate  length  =  120h). 
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The  B IEM  code  used  In  this  problem  basically  solves  the  Poisson 
equation  (3.56).  Compared  to  the  prooram  EDDY2  In  the  low  Reynolds  number 
case.  It  has  the  advantaae  that  only  boundary  elements  need  be  used.  This 
much  reduces  the  number  of  eauatlons  to  be  solved.  But  the  extension  of 
the  BIE  method  to  Include  the  nonlocal  effect  cannot  be  easily  achieved. 

d)  Dimensional  Analysis  of  the  Size  Effect  of  the  Width  of  a 
Square  Plate.  For  a  souare  plate  with  different  width/thickness  (D/h) 
ratios,  the  program  EDDY2  was  run  by  Hara  (Ref.  83).  By  using 
the  original  definition  of  the  maanetic  Reynolds  number  and  the  FE  mesh 
shown  in  Figure  3.42,  one  obtains  the  relations  between  the  modulus  of 
the  current  at  the  centroidal  points  of  the  elements  versus  the  parameter 
R  for  different  D/h  values.  The  exciting  maanetic  field  is  assumed  uniform 
here.  The  moduli  of  the  current  at  all  centroid  points  of  the  elements 
therefore  have  the  same  value. 

By  using  a  different  definition  of  the  Peynolds  number. 


rI  =  $  R  =  h  (3.57) 

6 

one  obtains  the  curves  shown  in  Figure  3.43.  At  small  R*  values  these 
three  curves  correlate  to  each  other  well,  but  the  differences  begin  to 
grow  as  the  R*  values  Increase. 

By  usina  the  same  curves  shown  in  Fiqure  3.42,  but  another  definition 
of  the  magnetic  Reynolds  number  R^, 


R11 


(3.58) 


one  obtains  the  curves  shown  In  Fiaure  3.43,  which  show  a  much  better 
correlation  for  the  high  R*1  values.  However,  the  aoreement  is  still 
not  as  close  as  that  obtained  for  the  one-dimensional  plate.  Figure  3.20. 


m 


•M 


Modulus  of  the  induced  current  density  at  the  element 
centroids  vs.  the  modified  magnetic  Reynolds  number 
R*  =  (D/h)R“  Dh/Unfi2)  for  a  square  plate  subject  to 
a  uniform  transverse  field. 
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This  rouah  parameter  study  shows  that  the  number  of  nondlmensional 
parameters  for  the  finite  plate  is  more  than  one,  and  that  they  are 
probably  related  to  each  other  In  some  fashion  that  cannot  he  revealed 
by  the  simple  chanaes  of  parameter,  Ea.  (3.57)  and  (3.58).  More  detailed 
dimensional  analysis  is  suaaested  to  study  the  size  effect  of  a  two- 
dimensional  plate. 


magnetic  Reynolds  number  Rll  =  R(D/h)^  =  (D/6)^/4tt  for  a  squ 
a  uniform  transverse  field. 


Chapter  4 


LINEAR  VIBRATION  OF  INFINITELY  LONG  MAGNETOELASTIC  PLATES 

The  magnetically  induced  vibrations  of  infinitely  lonq  conducting 
plates  are  considered  in  this  chapter.  The  plates  have  constant  span  £ 
alona  their  Infinite  length.  The  exciting  maanetic  fields  considered 
are  invariant  in  space  in  the  length  direction.  The  eauations  governlno 
the  coupled  system  in  this  case  are  Equations  (2.104)  and  (2.106).  The 
finite  element  formulation  used  for  analysis  is  presented  in  Section  4.1. 
The  staqqered  transient  analysis  is  applied  to  the  two  coupled  equations. 
The  computational  procedure  is  described  In  Section  4.2.  Numerical 
results  and  experimental  verifications  are  presented  in  Section  a. 3. 

The  various  coupling  effects  are  discussed  in  Section  4. a. 

The  study  described  in  this  chapter  has  been  limited  so  far  to  small 
transverse  flexural  motion  of  plates  subjected  to  pulsed  maqnetlc  fields. 
The  in-plane  force  effects  are  not  included.  The  total  coupling  effects 
and  the  influence  of  static,  uniform  magnetic  fields  on  the  free  vibrations 
of  the  plate  are  left  for  future  investigations.  The  couplino  effect 
through  rotational  motion  is  also  left  for  further  studies. 

4.1  Basic  Equations  and  Finite  Element  Formulation 

The  governing  eauations  for  Infinitely  lonq  conductina  plates  subject¬ 
ed  to  excitinq  magnetic  fields  invariant  in  space  in  the  lenath  direction 
are  derived  in  Chapter  Two  and  are  Egs.  (2.104)  and  (2.106).  In  this 
chapter  the  coupling  effects  in  the  induced  vibrations  are  limited  to  the 
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transverse  motion  of  the  plate.  The  governinn  eouations  for  this  case 
are  the  following  specialized  versions  of  Eg.  (2.104)  and  (?.106b): 


T? 


at 


(x-x' ) 


(x1) 


dx‘ 


h‘ 

T 


=  tir-lx  Bx°n 

and 

4  2 

n  3  w  .  .  3  w  .  3t|)  Do 

+Ph^"-Bx 


(4.1) 


(4.2) 


In  both  equations  the  couplino  effects  appear  on  the  right-hand  side  only. 

The  same  finite  element  formulation  procedure  as  described  in  Chapter 
Three  is  applied  to  Eq,  (4.1)  and  results  in 


<  It  [flI  •  n-e)  ts]  )  wtMt 


=  (  [A]  +  0  [s]  }  {^)t  +  e  ( R } t  +  (1-e)  {R)t+At 


where  { R >  is  oiven  by 


:3T-57<It  >;>!<"■ 


(4.3) 


(4.4) 


All  the  other  finite  element  matrix  and  vectors  are  the  same  as  defined  in 
Section  3.3. 

The  transverse  deflection  w  in  En.  (4.2)  is  approximated  by  the  usual 


cubic  model 
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w  =  T.  w.  =  lCEJ  (wl  (4.5) 

k=l  K  K 

where 

{w}T  =  lw1  e1  w 2  e?j  (4.6) 

and 

CE  =  1  +  2£3  -  3C2  CE  =  -x  (5  -  l)2 

CE  =  3^2  -  2e3  CE  =  -x  U2  -  £)  (4.7) 

with 
E  =  x/lE 

The  following  matrix  equation  of  motion  Is  then  obtained 

[M]  (w)  +  [K]  {w}  =  {F}  (4.R) 

in  which 

(FI  »  |  (C)  ||  B°  dx  (4.9) 

Usina  +  he  Newmark  inteqration  scheme  (Ref.  TOO),  one  has 

where 

[K1  =  aQ  [Ml  +  [K]  (4.11a) 

<fw  ■  ,F,t*st  *  mi  {a0  (w)t  *  a?  rt)t  *  >3  <w>t)  H.'ib) 


with 
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a3  *  k  ■ 1f 

6  >  0.50  a  >  0.25  (0.5  +  <5)? 

The  trapezoidal  rule  Is  used  for  all  the  studies  with  a  =  and  6  ®  ^  . 


4.2  Staggered  Transient  Analysis  and  Computational  Procedure 

The  eguations  (4.1)  and  (4.2)  may  he  solved  simultaneously  in  principle. 
In  that  case,  the  discretized  eguations  will  have  the  followino  form 


[A]  {*}  -  [s]  HI  ■  {R>  -  [D]  {w} 


[M]  {wl  +  [K]  {wl  =  [P]  {*} 


where 

(R1  --  oh 

[D]  =  oh 


>  TT  dx 

E  9g° 

{NE}  i||— j  B°  dx  +  oh  f  (NE>  lCEj  dx 


[P]  .  }  (CE)  l|Si|  Bj  dx 
E 


Define  a  state  vector  {A}  with 
{a}T  »  i*  wj 

Eg.  (4.12)  may  then  be  put  into  the  following  form 
[M]  (A)  +  [D]  (A)  +  [K]  (A)  =  {R} 


where 


[Ml  = 

0  0 

[D]  = 

A  D 

[K]  = 

-s  O 

0  M 

» 

0  0 

» 

-P  K 

(4.12a) 

(4.12b) 

(4.13a) 

(4.13b) 

(4.13c) 

(A. 14) 

(4.15) 
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{R}T  ■=  [R  01 

Eq.  (4.15)  may  be  Intearated  directly  to  oive  the  solution  simultaneously. 

There  are  several  disadvantages  of  the  simultaneous  solution  scheme 
however.  Besides  the  inefficiency  of  the  solution  of  a  laroer,  unsymmetrlc 
matrix  equation,  the  same  finite  element  mesh  must  be  used  for  both  sub¬ 
systems  for  the  treatment  of  the  couplino  terms.  Since  a  finer  mesh 
usually  has  to  be  used  for  the  magnetic  subsystem,  this  will  result  in  an 
excessive  number  of  structural  elements  for  the  coupled  problems.  The 
other  consideration  is  the  different  time  characteristics  of  the  maanetic 
and  mechanical  subsystems.  Since  it  is  the  structural  response  rather 
than  the  stress  wave  motion  In  the  structure  that  is  of  interest  in  this 
study,  one  would  like  to  know  the  range  of  the  frequency  of  the  field  with¬ 
in  which  the  mutual  interactions  are  Important.  At  hiaher  frequency  of 
the  drlvinq  field,  one  may  orefer  different  time  steps  for  the  two  sub¬ 
systems  because  the  maanetic  force  will  be  more  like  an  impulse  to  the 
structure.  For  these  reasons  a  staaaered  solution  scheme  is  attractive. 
Since  the  two  governina  eauations  are  only  weakly  coupled,  this  approach 
Is  feasible  and  is  adanted  in  this  work.  Different  hut  conformable  meshes 
are  used  for  the  two  subsystems.  The  solution  state  of  the  coupled  problem 
is  advanced  by  sequentially  executino  Eauations  (4.3)  and  (4.10).  The 
Interaction  terms  appear  in  each  as  an  external  force  effect.  Temporal 
predictors  may  be  used  to  calculate  these  force  terms  more  accurately 
(Ref.  85).  In  the  present  study,  no  predictor  is  used.  The  accuracy  is 
improved  by  usina  a  smaller  time  step. 

A  Fortran  program  EDDYBFAM  has  been  developed  for  the  linear  vibration 
of  infinitely  lono  plates.  The  computational  scheme  employed  is  basically 
a  combination  of  the  transient  eddv  current  prooram  FDDY1T  and  the  struc- 


tural  analysis  prooram  BEAM.  Different  levels  of  couplino  of  these  two 
proarams  may  be  specified  in  the  input.  The  functions  of  these  oroqrams, 
which  use  different  but  conformable  finite  element  meshes,  is  shown  in 
Fiqure  4.1 . 

The  first  phase  of  the  proqram  is  the  initiation  of  the  two  sub- 
proorams  EDDY1T  and  BEAM  and  the  specification  of  the  type  of  couplino. 

If  the  parameter  LINK=0,  the  prooram  provides  the  transient  eddy  current 
calculation  only.  If  LINK=1 ,  the  mapnetic  force  calculated  on  the  riaid 
plate  Is  used  in  the  structural  analysis  prooram.  The  VxB  effect  on  the 
eddy  current  will  not  be  considered.  If  LINK=?,  two-way  interactions  are 
Included  and  the  solution  corresnonds  to  the  coupled  problem  referred  to 
the  initial  undeformed  position  of  the  plate.  For  LINK=3,  the  vector  f R> 
in  Fq.  (4.4)  is  calculated  on  the  deformed  position  of  the  nlate  at  the 
previous  time  step  usino  the  unit  normal  vector  of  the  undeformed  olate. 

The  maqnetic  force  calculated  is  still  applied  to  the  orioinal  nositlon  of 
the  plate  for  the  analysis  of  the  structural  response.  A  limited  follower 
force  nature  of  the  maonetic  force  is  then  obtained  in  the  analysis.  The 
results  obtained  from  the  different  levels  of  couplino  between  the  two 
subsystems  serve  to  indicate  the  relative  importance  of  the  various  inter¬ 
action  effects. 

The  major  computational  task  is  performed  by  the  prooram  FDDY1T.  As 
shown  in  Fiqure  4.1,  this  computation  consists  of  the  time  intearation  of 
the  eddy  current  diffusion  eouation.  At  each  time  step,  the  induced  current 
1^  and  maqnetic  force  Fz  are  calculated  at  each  inteqration  point.  The 
maqnetic  forces  are  then  transferred  to  the  structural  analysis  prooram 
BEAM.  The  maqnetic  forces  are  inteorated  numerically  and  the  effective 
consistent  nodal  force  vector  calculated  to  oive  the  displacement  at  the 


EDDYBEAM 


Indicate  type  of  analysis 

LINK  =0,  1,2,3: 
Initiate  EDDY1T,  BEAM 


EDDY1T 


Figure  4.1  Computational 

Scheme  of  EDDYBEAM 
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new  time  step.  The  Vx£  values  are  then  evaluated  at  the  inteoration  points 
bv  interpolation.  If  the  chanqe  of  distance  between  the  plate  and  coil  is 
considered  (LINK=3),  the  calculation  of  VxR  is  performed  at  the  new  position 
of  the  plate.  The  same  alternative  is  applied  to  the  calculation  of  the 
rate  of  chanoe  of  the  external  maonetic  field  £°.  These  calculated 
values  are  then  used  to  beain  the  next  time  step. 

In  addition  to  the  Input  Information  described  in  sub-subsection  3 . 3 . 1 . b 
necessary  for  the  transient  eddy  current  calculation,  qeometrlc  and  material 
attributes  for  the  structural  problem  must  be  Input  If  LINK  is  not  zero. 

Since  a  finer  mesh  usually  is  needed  for  an  accurate  calculation  of  the 
induced  current  and  force,  different  but  conformable  meshes  may  be  used  for 
the  maanetlc  and  mechanical  subsystems.  The  conformahil Ity  of  the  two 
different  finite  element  meshes  is  achieved  by  dividlno  each  structural 
element  Into  several  eoual-lenath  field  elements.  The  number  O'f  divisions 
may  be  different  for  each  structural  element  and  must  be  input  as  part  of 
the  mesh  Information. 

The  same  size  of  time  step  is  used  for  the  calculation  of  both  the 
maanetlc  and  mechanical  subsystems.  A  lumped  mass  matrix  is  used  for  the 
mechanical  problem.  Initial  conditions  and  mechanical  loads  may  be  input 
to  the  proqram.  Mesh  generation  and  plotting  capabilities  for  the  time 
history  and  spatial  variation  of  displacement  and  bendlna  moment  are 
provided  as  options. 

The  output  from  the  ED0Y1T  portion  of  the  nrooram  Is  described  In 
Chapter  Three.  In  addition,  displacement  velocity,  and  acceleration  at 
each  structural  node  together  with  the  bendlno  moment  at  specified  points 
are  produced  as  output  of  the  subproaram  REAM.  These  are  printed  out  after 
the  output  of  the  subproaram  E0DY1T  for  each  time  step.  Displacements  and 


bending  moments  may  be  plotted  uslno  a  library  subroutine  If  the  plotting 
option  Is  chosen. 


4.3  Numerical  Results  and  Experimental  Verifications 

The  program  EOPYBFAM  has  been  applied  to  the  mapnetical ly  induced 
vibrations  of  aluminum  and  stainless  steel  (SS  304)  plates.  The  properties 
of  these  two  materials  are  listed  below. 


Material 

o  (Amp/Volt.m) 

v  (N/Amp?) 

p  (N.Sec?/m4) 

E  (N/m?)  v 

aluminum 

3.80  x  107 

An  x  10”7 

?.7  x  103 

7.0  x  1010  0.34 

SS  304 

1.30  x  IQ6 

4w  x  10"7 

8.0  x  103 

2.0  x  1011  0.30 

The  study  so  far  has  been  limited  to  single  pulsed  currents.  To  have 
sufficient  Interaction  time,  the  time  characteristics  of  the  electro¬ 
magnetic  and  mechanical  subsystems  are  compared  below  for  cantilever  plates. 
The  information  is  used  to  desion  meaningful  numerical  and  experimental 
tests. 

Let  f  =  frequency  of  vibration  of  the  nth  mode  of  the  plate,  and 
f^  =  freouency  of  the  driving  current  In  the  coll. 

Then 


In  which  E  Is  the  effective  Young's  modulus  for  long  plates.  Let  the  skin 
depth  6  be  N  times  the  thickness  of  the  plate  h. 


2  2 

N  n 


1 

T 


t‘C1 


fE  -  1 

f - n-J 

1,  no  fr  lv 

Matching  the  two  freouencies,  f  =  f^,  one  obtains 


t 


2 


-J-  ^?K  ?  no  /I  h3 

2  ,r?  n  '  0 


(4.17) 


(4.  IP) 


The  lenoth  and  thickness  of  the  nlate  for  different  combinations  of  mode  and 
skin  depth/thickness  ratios  are  shown  in  Table  4.1.  The  followino  values 
are  used  in  the  calculation  of  the  tabled  values  (Ref.  86). 


n 

1  2 

3  4  5  6  7 

8 

K 

n  i 

1.875  4.684 

7.865  10.006  14.137  17.270  20.420 

23.562 

The  frequencies  and  half-periods  of  the  oscillations  for  both  the  maanetic 
field  and  the  structure  correspondino  to  the  values  in  Table  4.1  are  calcu¬ 
lated  in  Table  4.2.  The  problem  tested  is  determined  by  choosing  suitable 
lenoth  and  thickness  so  that  the  frequency  of  the  driving  current  matches 
that  of  the  vibration  beam  for  a  selected  mode. 

The  transient  eddy  current  prooram  EDDY1T  has  been  verified  in  Chanter 
Three.  The  structural  analysis  prooram  BEAM  has  also  been  tested  for 
several  dynamic  problems.  One  problem  tested  Is  the  free  vibration  of  a 
simply  supported  beam  subjected  to  a  sinusoidal  initial  velocity  shown  In 
Figure  4.2.  With  four  beam  elements  for  half  of  the  beam  and  intearation 
by  the  trapezoidal  rule  version  of  the  Newmark  method,  the  results  show 
qood  agreement  with  the  theoretical  solution  for  both  choices  of  time 
step,  althouah  for  At  =  T/20  some  period  elonoation  (%  IT)  is  apparent. 

The  coupled  problems  studied  have  the  oeneral  arranoement  shown  in 
Fiaure  4.3  with  one  exciting  wire.  One  problem  analyzed  is  for  an  aluminum 
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Table  4.1  lenoth  and  Thickness  for  Different  Combinations 
of  Mode  Shane  (n)  and  Skin  Depth  (N  =  s/h) 

Aluminum  (At)  :  no  =  47.78  Sec/m  ,  JEJT  -  5000  m/Sec 


SS  304  :  no  =  1.8  Sec/m'  ,  »T7^  =  5000  m/Sec 


3,  a  7 
Q.  83 
If  .48 

f .  04 
15.13 

25.32 


8.45 

21.15 

35.30 


7.85 
10. ff 
32.  on 

12.0Q 
30. 2f 
50. f4 


If  .00 
*2.30 
0.78 


11.78 

?o.4Q 

40.34 

18.13 

45.40 

76.00 


25.3 

63.4 
106.17 


15.70 
30. 
65.,  . 

24.18 

60.53 

101.?o 


33.7 

84.5 

141.55 


7.63 
-  '.14 

82.24 

30.22 

75.66 

126.61 


176.04 


Table  4.2  Frequency  and  Mil f-perlod  of  Oscillation  fnr  Different 
Combinations  of  Thickness  (h)  and  Skin  Depth  (N  =  f/h) 


Aluminum  (AO 

:  iio  =  4  7 

.78  Sec/m2 

h 

_ 

f*  (1/sec)  /  td+ 

(msec) 

(mm. ) 

1 

2 

3 

4 

5 

1.0 

f 

6.66xl03 

1 .67x1 03 

7.40x1 0? 

4.16xl02 

2.66x10" 

0.075 

0.300 

0.675 

1.201 

1 .876 

1.5 

f 

2.  <36x10' 

7.40xl02 

3.2Qxl02 

1 .85x1 02 

1 ,18xl0? 

*d 

0.160 

0.675 

1.52 

2.702 

4.27 

2.0 

m 

1 ,67xl03 

4.16xl02 

1 .85xl02 

1.04X10? 

6.66X181 

0.300 

1.200 

2.702 

4.803 

7.505 

SS  304  :  v>n  =  1.8  m/Sec2 


1 

4.44x10’’ 

0.001 

1.11x10' 

0.005 

2.83x10* 

0.018 


*  f  = 


f 


E 


"  f'*'  (1/sec)  /  td+  (msec) 

2  3  4 

T.  llxl 05  4.94x10*  ?.78xina 

0.005  0.010  0.018 

2.78x10*  1 .23x10*  6.04x10* 

0.018  0.041  0.072 

7 . 07x1 03  3 . 14x1 03  1  .77x1 025 

0.071  0.159  0.283 


*d  =  ^ 


pulse  duration 


5 


1.78x10 


4 


0.028 

- , 

4.44xlOJ 

0.113 

1. 1 3x1  O' 
0.442 


Note:  f  can  be  made  eoual  to  f  by  adjustlna  the  lenoth  t. 
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plate  with  t  =  427  min,  a  =  33?  mm,  h  «=  2.1  mm,  and  o  =  9.1  nrn.  The  excit- 
ino  wire  Is  purposely  placed  above  the  nodal  point  of  the  second  vibration 
mode  of  the  cantilever  nlate.  Because  of  the  much  sharper  variation  of 
the  electromaqnetic  field  variables,  more  field  elements  than  beam  elements 
are  used  In  the  analysis.  The  meshes  of  the  two  suhsvstems  arp  made  con¬ 
formable  by  dividino  each  beam  elempnt  into  several  equal-lennth  field 
elements  as  shown  in  Fioure  4.4a.  Thp  freounncies  and  periods  of  the  f1rst 
three  modes  of  the  plate  in  this  problem  are: 


1  N 

1 

2 

3  1 

f  ( 1 /sec) 

o.51 

59.60 

156.87 

T  (sec) 

0.1052 

0.0168 

0.0060 

A  half-sine  pulse  with  3  ms  pulse  duration  and  a  500  Amp  drivino  current 
has  been  chosen.  The  pulse  duration  is  made  eoual  to  half  the  period  o' 
the  third  mode  of  vibration.  The  results  of  the  analysis  are  shown 
Fiaures  4,5  and  4.6  and  will  be  discussed  in  the  next  section. 

Experimental  verification  has  been  attempted  on  a  finite  lennth  bea1” 
plate  with  the  same  cross  sectional  properties  and  the  arrannement  of  the 
excltino  wire.  This  attempt  failed,  however,  because  of  the  small  lenoth- 
to-width  ratio  of  the  plate  used  in  the  experiment.  The  size  effect  on 
the  eddy  current  distribution  was  investioated  experimentally  (Ref.  101). 
The  affected  distance  from  the  ends  of  the  plate  was  found  to  be  about 
eoual  to  the  width  of  the  nlate.  A  new  conf iouration  that  can  be  better 
simulated  by  the  one-dimensional  analysis  was  then  desloned.  The  numerical 
and  experimental  results  for  this  second  conf iouration  are  described  below. 
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Figure  4.5  Nonlocal  predictions  of  eddy  current  density,  magnetic 
force,  and  transverse  displacement  at  t  =  0.9  ms  for 
cantilever  plate  with  single-wire  exciting  coil. 

Figure  4.4(a). 


Displacement  (mm)  Force  (N/m)  Current  density  (xIO 


Figure  4.6  (a),  (b)  Induced  current  and  force  under  the  wire  during 

the  pulse  duration,  (c)  vertical  displacement  at  node  6 
in  Figure  4.4(a)  during  1/4  the  period  of  the  first  mode. 


The  new  plate  has  a  lenoth/width  ratio  about  4.2  and  Is  shown  In 
Figure  3.31.  The  dimensions  of  the  new  conf iouration  are  ?  =  14Q  mm, 
a  =  117  rim,  h  =  1  mm,  and  e  =  4,r  mm.  The  transient  eddy  current  calcula¬ 
tion  is  first  verified  for  this  new  configuration  using  the  search  coil 
technique  for  e  =  11  mm.  The  results  have  been  presented  In  subsection 
3.3.2.  For  the  coupled  problem,  the  finite  element  meshes  shown  In 
Figure  4.4b  with  e  =  4#$  mm  are  used  for  analysis.  In  the  experiment, 
two  low  inductance  strain  oaoes  are  applied  at  the  position  of  node  3  In 
Figure  4.4b.  The  time  history  of  the  bending  moment  recorded  is  compared 
to  the  calculated  result  as  an  experimental  verification. 

The  freguencles  and  periods  of  the  first  four  modes  of  an  infinitely 
long  cantilever  plate  with  this  thicknpss  are 


N 

1 

f  (1/sec) 

T  (sec) 

37.6 

26. 6x10' 3 

In  the  free  vibration  test  of  the  actual  fini te-lenoth  plate,  the  same 
frequencies  and  periods  have  been  observed. 

Three  different  pulsed  currents  are  used  to  study  the  induced  motions 
of  the  plate.  The  study  of  pulsed  problems  has  been  motivated  by  potential 
applications  in  magnetic  forming,  pulsed  electric  maonets,  and  structural 
problems  in  hiqh  energy  devices.  The  durations  and  currents  of  the  driving 
pulses  are 

Pulse  No.  1  2  3 

I  (Amps.)  2700  7000  8500 


td  (sec.) 


5.0x10' 


2.5x10' 


1.2x10' 


167 


The  actual  pulsed  currents  are  Input  numerically  into  the  prooram  FDDYBFAM. 
The  calculated  bending  moments  at  the  position  of  the  strain  cage  are 
compared  with  the  measured  values  for  each  case  separately.  Floures  4. 7, 
4.9,  and  4.11  show  the  calculated  bendina  moments  for  the  first,  second, 
and  third  pulse,  respectively.  Floures  4.8,  4.10,  and  4.1?  show  the 
recorded  voltage  readlnos  from  the  strain  aaqes  for  each  of  the  pulses 
used.  These  voltaoe  readlnos  have  been  calibrated  and  compared  with  the 
numerical  results.  The  discussions  of  these  numerical  and  experimental 
studies  of  the  coupled  problems  are  oiven  in  the  next  section. 

4.4  Coupling  and  Nonlocal  Fffects 

Because  of  the  many  geometric,  material,  and  time  parameters  involved, 
the  dimensional  analysis  of  the  transient  coupled  Droblem  is  more  difficult 
than  that  of  the  steady-state  problem  and  has  not  been  attempted  yet. 
However,  from  the  limited  number  of  problems  studied  in  the  last  section, 
some  conclusions  about  the  various  coupling  and  nonlocal  effects  may  be 
drawn. 

The  results  shown  in  Fioure  4.6  for  the  first  problem  Indicate  that 
the  velocity  and  displacement  of  the  conducting  plate  both  have  a  dampino 
effect  on  the  magnetic  and  mechanical  subsystems.  The  magnitudes  of  the 
current,  force,  and  displacement  have  been  reduced  and  their  phases  changed, 
as  can  be  seen  In  Figures  4.5  and  4.6.  The  damping  also  reduces  the  energy 
transfer  from  the  magnetic  subsystem  to  the  mechanical  subsystem.  The 
efficiency  of  the  energy  transfer  may  be  obtained  by  calculating  the  power 
supplied  by  the  current  source  and  the  total  kinetic  and  strain  energies 
of  the  plate.  An  Indication  of  the  influences  of  various  coupling  effects 
on  the  energy  transfer  may  be  seen  from  Figure  4,6c  which  shows  the  long- 
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I - 1  Pulse  duration 

I - 1  3rd  mode  period 


Figure  4.7  Comparison  of  the  numerical  and  experimental 
results  of  the  bending  moment  at  node  3, 
Figure  4.4(b),  for  the  5  ms  pulse. 
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Figure  4.9  Comparison  of  the  numerical  and  experimental  results 
of  the  bending  moment  at  node  3,  Figure  4.4(b),  for 
the  2.5  ms  pulse. 
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e  4.12  Low  inductance  strain  gage  readings  for  the 
1.2  ms  pulse  for  two  different  time  scales. 
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time  solutions  of  the  displacement  at  node  6  (x  =  33(1  mm).  When  the 
displacement  reaches  its  maximum  value,  the  eneroy  density  of  the 
mechanical  subsystem  consists  entirely  of  strain  enemy.  The  pffects  of 
the  various  couplinos  in  reducino  the  total  enerny  density  at  that  noint 
is  obvious  from  the  fioure. 

The  nap  distance  between  the  wire  and  the  mid-plane  of  the  plate  Is 
approximately  four  times  the  thickness.  The  skin  depth  of  a  steady  state 
current  with  the  same  neriod  will  be  three  times  the  thickness  of  the 
plate.  The  small  nap  distance  and  the  nonlocal  effects  tend  to  concentrate 
the  mannetic  force  at  the  nodal  point  of  the  second  vibration  mode.  As  a 
result,  the  Induced  motion  of  the  platp  consists  mainly  of  the  first  and 
the  third  modes.  The  superposition  of  the  third  mode  on  the  first  mode  Is 
clearly  shown  In  Fiqure  4.6c. 

The  nonlocal  effect  is  still  sinniflcant  in  this  problem.  The  non¬ 
zero  averane  pushino  force  can  he  easily  seen  from  Flaure  4.6b.  A  pull  inn 
force  is  developed  close  to  the  end  of  the  Dulse.  Since  the  time  variation 
of  the  maanetlc  field  determines  the  induced  current  and  force,  it  is 
possible  to  deslon  a  drlvina  current  to  Generate  a  laroe  pullino  force. 
Instead  of  a  pushino  one,  on  the  conductor. 

Three  different  drivinq  currents  have  been  applied  in  the  second 
problem.  Since  the  position  of  the  strain  oaoe  is  far  away  from  the  wire, 
the  effects  of  the  various  kinds  of  couplinq  are  indlstinouishable  and 
only  one  curve  is  plotted  for  each  case.  Sood  quantitative  aoreement  has 
been  obtained  in  all  cases.  The  quantitative  difference  is  small  in  the 
lono  pulse  duration  case  and  larqpr  in  the  shorter  pulse  duration  cases. 

The  differences  may  be  partly  caused  by  the  flexibility  of  the  clamped 
end  support.  The  correlation  between  the  numerical  and  experimental  results 


is  considered  acceptable  in  an  experiment  of  this  kind  because  of  the 
end  effects  of  the  flnite-lenqth  plate,  the  neolect  of  the  structural 
and  rotational  maanetlc  damoino,  the  finite  size  of  the  strain  aaoe  and 
the  excitino  wire,  and  the  precision  of  the  experimental  measurements. 
These  comparisons  verify  the  model  and  the  numerical  results  of  the 
coupled  problem. 

The  lenqth  of  the  pulse  duration  has  an  obvious  effect  on  the  mode 
of  the  Induced  vibration.  In  the  5  msec  pulse  duration  case,  mainly  the 
first  vibration  mode  is  excited  as  can  be  seen  from  the  osclllascope 
pictures,  Fiqure  4.8.  As  the  pulse  duration  decreases,  the  participation 
of  the  hioher  modes  Increases  as  shown  in  Floures  4. in  and  4.12  durino 
half  the  period  of  the  first  mode  of  vibration.  Note  that  in  Fioure  4.11 
some  hlqher  frequency  peaks  in  the  experimental  result  cannot  be  represent 
ed  in  the  numerical  result.  This  is  because  the  same  finite  element  mesh 
is  used  for  all  three  pulsed  currents.  To  represent  those  hinher  freoupnc 
peaks  more  accurately,  more  beam  elements  would  need  to  be  included  In  the 
finite  element  mesh, 

A  schematic  comparison  of  the  pulse  duration  and  the  third  mode 
period  has  been  included  In  Fioures  4.7,  4.9,  and  4.11  to  demonstrate 
this  effect.  More  hioher  frequency  modes  will  be  excited  as  the  pulse 
duration  further  decreases.  The  maonetic  force  in  this  case  will  act 
more  like  an  impact.  The  nonlocal  effect  will  also  be  more  important  as 
the  short  pulse  incU/des  more  hioh  freouency  components.  The  nonzero 
averaoed  pushlno  force  from  the  ronlocal  pffect  durino  the  pulse  duration 
serves  as  the  Impact  force  in  such  a  case. 

A  quantitative  description  of  the  pffect  of  pulse  duration  on  the 
induced  vibration  Is  difficult  to  draw  from  these  studies  presented.  Onp 


reason  Is  that  the  current  Is  vol tane-control 1 ed  and  a  different  amount 
of  enemy  has  been  applied  for  each  pulsed  current.  To  studv  the  effect 
of  nulse  duration  on  the  distribution  of  enerov  transferred  into  different, 
vibrational  modes,  the  same  electromagnetic  enemy  should  be  input  to  the 
whole  svstem,  and  the  mutual  inductance  between  the  plate  and  wire 
considered.  A  parameter  studv  of  the  effect  of  pulse  duration  can  then 
be  conducted  throuah  modal  analysis  of  the  vibration  plate. 

An  estimation  of  the  effect  of  the  difference  between  the  time 
characteristics  o*  the  two  subsystems  on  the  Induced  vibration  can  be 
obtained  in  the  followlno  sense.  xhe  magnetic  *orce  mav  be  Integrated 
during  the  pulse  duration  and  applied  on  the  plate  as  an  Impact  force. 

The  mechanica1  epnnnv  ’ "  th,-s  case  w-'l1  be  eoual  to  the  Hnetlc  energy 
in  the  eouivalent  inlfia1  velocity  problem.  Thp  mechanical  enerov  in  the 
coupled  problem  ^av  bP  calculated  *rom  thp  kinetic  and  elastic  pneroips 
of  thp  plate  after  ‘hp  pulse  has  bppn  applied.  Since  Jn  the  very  short 
pulse  case  the  induced  eddv  current  will  screen  thp  applied  magnetic  field 
out,  the  maonptic  feme  will  act  more  life  a  naonetic  pressure  as  the  puIsp 
duration  decreases.  The  comparison  hetwepn  the  two  mechanical  pnpmies 
thus  calculated  versu'  the  ratio  of  Dulse  duration  to  various  fundamental 
periods  of  vihratior  mav  indicate  the  effect  of  nulsp  duration  on  the 
induced  vibrations.  A  more  rigorous  studv  of  such  effects  should  be  based 
on  thp  enerov  conversion  between  thp  magnetic  and  mechanical  subsystems 
based  on  the  same  ener  ■’v  input,  as  described  In  the  last,  paragraph. 

Thp  studv  of  the  coupled  problems  in  this  chanter  h,  s  beer  limited 
to  the  case  of  one  exciting  wlm.  Thp  various  coinlino  and  nonlocal 
effects  are  of  course  problem  dependent  and  should  bp  Investigated  for 
other  types  of  applied  field,  Manv  different  types  pf  magnetic  field. 
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however,  can  be  Generated  from  superpositions  of  the  one-wire  fields. 
The  study  of  the  one-wire  problem  thus  serves  as  a  prototype  for  future 
studies. 


1 

i 


I 


* 


Chapter  5 


NONLINEAR  VIBRATION  OF  INFINITELY  LONG  MAGNFTOELASTIC  PLATES 

The  magnetically  Induced  vibrations  of  Infinitely  lono  conductlna 
plates  with  finite  deflections  but  small  rotations  are  considered  In  this 
chapter.  Beside  the  Interest  In  studyina  the  motion-dependent  nature  of 
magnetic  forces,  such  studies  may  find  applications  In  maonetlc  formlno 
or  In  devices  like  mannetodynamlc  circuit  breakers.  Moreover,  Geometric 
nonllnearltles  must  be  Included  for  Investloatlons  of  structural 
stability. 

The  plate  and  the  excltlno  maonetlc  field  treated  have  the  same 
oeneral  arrangement  as  In  Chapter  Four.  The  eouatlons  oovernlno  the 
coupled  system  now  are  Eouatlons  (2.102)  and  (2.6f)  to  (2.67).  The  basic 
equations  and  the  finite  element  formulation  are  presented  In  Section  5.1. 

An  Incremental-Iterative  procedure  Is  suggested  In  Section  5.2  for  the 
coupled  nonlinear  transient  system.  The  Implementation  of  the  procedure 
In  this  study,  however,  met  convergence  problems  because  of  the  interac¬ 
tion  between  the  time  Integration  operator  and  the  approximate  solution  pro¬ 
cedure  for  the  equations  of  motion.  These  difficulties  are  discussed  in 
Section  5.3  together  with  some  successful  results  obtained  for  static, 
nonlinear  elastic  problems. 

5.1  Basic  Equations  and  Finite  Element  Formulation 

The  linearized  equations  of  motion  oovernlno  the  nonlinear  mechanical 


system  are 
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/,  f  0  .  e  .  6e  .  *dV  +  /  6  n  tdV  *  t+AtR  -  /,  «e  tdV 

ty  ctRy6  aP  *y<5  ty  t  aP  aP  ty  t  aP  aB 

(2. 66a) 


t+AtR  *»  /.  t+Att  5  u  tdA  +  /  %  (t+Atf  .  u  )  6  u  tdV  (2. 66b) 

a  a  ty  a  a  a  ' 


A  nonlinear  beam  element  developed  by  Bathe  and  Bolourchl  (Ref.  87)  Is 
used  for  the  updated  Lagranqlan  (UL)  finite  element  formulation.  The 
components  of  displacement  Increment  are  approximated  by 

6  1 

u.  -  z  h?  U|f  1  *  1,2  (5.1) 

1  k=l  k  K 

In  which  the  h^  are  the  interpolation  function  for  the  kth  node  In  the 
element  and  the  uk  are  the  nodal  point  displacement  Increments  In  the  local 
axes  at  time  t.  Equation  (5.1)  Is  used  to  evaluate  the  linear  and  nonlinear 
parts  of  the  strain  tensor  In  Eo.  (2.66)  for  each  element.  The  formulation 
procedure  results  In  the  following  Incremental  finite  element  stiffness 
equations 


([\]  +  (*%])  (u>  *  (t+AtR)  -  (V)  (5.2) 

In  which  are  linear  and  nonlinear  strain  Incremental  stiff¬ 

ness  matrices  for  the  conflauratlon  at  time  t,  (t+AtR)  Is  the  vector  of 
Inertial  effects  and  externally  applied  nodal  loads  at  time  t+At;  and 
{^F}  Is  the  vector  of  nodal  point  forces  equivalent  to  the  element  stresses 
at  time  t. 

The  linear  stiffness  matrix  [^K^]  Is  computed  as 

[\]  -  ft  [^L^tC]  *dV  (5.3) 

In  which  [tBL]  Is  the  linear  strain-displacement  transformation  matrix  for 
the  configuration  at  time  t,  and  [C]  Is  the  material  property  matrix,  [^KjJ 
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Is  the  same  as  the  usual  linear  beam  stiffness  matrix  (Ref.  88). 

The  nonlinear  stiffness  matrix  [*<^1  accounts  for  the  nonlinear 
effect  of  the  stresses  at  time  t  and  Is  calculated  by 

[  V  ■  \  'V'Vi  rV> t<lv  (5-4) 

In  which  t*BNL]  Is  the  nonlinear  strain-displacement  transformation  matrix, 
and  [tt]  Is  the  matrix  of  Cauchy  stresses  In  the  conflouratlon  at  time  t. 

The  vector  of  nodal  point  forcea(tF)  accounts  for  the  linear  effects 
of  the  stresses  at  time  t  and  Is  computed  by 

{tf}  s  \  ft\]T{t*T}  tdV  (5*5) 

A 

In  which  {^xl  Is  the  vector  of  Cauchy  stresses  In  the  conflouratlon  at 
time  t. 

The  vector  of  externally  applied  nodal  loads  {t+AtR>  for  the  conflaura- 
tlon  at  time  t+At  Is  obtained  by  the  finite  element  evaluation  of  Eo.  (2.44) 
uslnq  Eq.  (2.66b).  The  deformation  dependent  nature  of  the  force  Is  pro¬ 
posed  to  be  treated  In  an  Iterative  way  In  the  computational  procedure. 
Details  of  the  various  element  matrices  and  vectors  are  presented  In 
Appendix  C. 

A  moving  curvilinear  coordinate  c  Is  Introduced  on  the  mid-surface  of 
the  plate  for  the  eddy  current  problem.  Figure  5.1.  As  a  first  step  In  the 
study  of  the  geometrically  nonlinear  coupled  problem,  the  self-field  effect 
Is  neglected.  Eq.  (2.102)  In  this  case  becomes  for  one-dimensional 
problems: 

.2,  , 

H a  oh  +  k \] 


(5.6) 


a**— xmani  ii  tit  iiinaa 


fj»  unite 
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Figure  5.1  Fixed  and  moving  coordinate  systems  on 
the  mid-surface  of  a  cantilever  plate. 


Figure  5.2  Configurations  of  a  cantilever  plate  at 
three  consecutive  Iteration  steps. 
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In  which  r,  denotes  the  unit  normal  vector  on  the  mid-surface  of  the 
deformed  plate. 

Using  the  Galerkln  formulation  and  the  linear  shape  function  approxi¬ 
mation 

r  2  r 

<T  =  I  „  (5.7) 

k=l  k  K 

one  may  discretize  Eg.  (5.6)  Into 


in  which 


(5.8) 


(5.9a) 


-  «i.  |  [6°  ‘  fj-  <«)z1  de  (5.%) 

After  the  stream  function  Is  solved,  the  magnetic  body  force  Is  calculated 
from 


t  t+At., 

p 


f  =  .  11  B°  e.  +  R?  n 
%  ?r  n  A-t  3c  t  ^ 


(5.10) 


This  expression  Is  used  in  Eg.  (2.6fb)  for  the  calculation  of  the  externally 
applied  load  vector  (t+AtR). 


5.2  Staggered  Transient  Analysis  and  Computational  Procedure 

Except  for  the  deformation- and  velocity-dependent  magnetic  force,  the 
solution  of  Eq.  (5.2)  may  follow  the  standard  Incremental  Iterative  pseudo¬ 
force  procedure  (Ref.  87).  With  an  Implicit  time  Integration  scheme,  this 
Involves  a  calculation  of  the  unbalanced  force  and  an  iterative  procedure 
to  achieve  the  equilibrium  among  the  Internal  element  stresses,  the 
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externally  applied  nodal  point  loads,  and  the  Inertia  force  at  time  t+at 
(Ref.  89).  For  deformation- Independent  loads,  such  a  procedure  Is  stralqht- 
forward.  For  some  types  of  follower  forces,  such  as  pressure  loadlnns,  a 
load  stiffness  matrix  mav  he  added  to  the  tanoent  stiffness  matrix  to 
account  for  the  deformation-dependent  nature  of  the  force  (Ref.  *>0).  For 
the  deformation-  and  velocity-  dependent  magnetic  force,  this  approach  Is 
not  applicable.  In  this  work,  an  updated  calculation  of  the  maqnetlc 
force  within  each  pseudo-force  Iteration  step  is  sugoested  to  achieve  the 
simultaneous  satisfaction  of  both  the  eddy  current  eauatlon  and  the  eoua- 
tions  of  motion.  The  procedure  for  this  calculation  Is  described  below. 

Consider  the  state  of  motion  of  the  plate  at.-f-th  Iteration  shown  In 
Figure  5.2.  The  maonetic  force  has  been  calculated  based  on  the  motion  of 
the  plate  at  the  j-lst  Iterative  step  and  applied  to  the  conflnuratlon  of 
the  plate  at  that  step.  If  the  unbalanced  force  calculated  Is  such  that  a 
prescribed  convergence  tolerance  is  violated  (Ref.  89),  the  eddy  current 
equation  Is  solved  aqain  based  on  the  motion  and  position  at  the  ,1-th 
iterative  step.  The  unbalanced  force  is  recalculated  with  this  newly 
calculated  maonetic  force  and  applied  to  the  configuration  at  the  j-th 
Iterative  step.  The  same  nonlinear  stiffness  matrix  [^K^]  and  nodal  force 
vector  {tF}  are  employed  In  each  iteration.  The  procedure  Is  continued 
until  the  tolerance  Is  satisfied,  l.e.,  until  the  error  In  equilibrium  or 
the  change  In  Incremental  enerqy  Is  acceptably  small  (Ref.  8Q). 

The  trapezoidal  rule  version  of  the  Newmark  method  (Ref.  100)  Is  used 
to  Inteqrate  the  eouatlons  of  motion.  A  ouasl-Newton  method  Is  used  for 
the  Iterative  procedure  (Ref.  89).  A  force  and/or  energy  tolerance  Is 
used  for  convergence  checks  (Ref.  89).  The  Incremental  Iterative  computa¬ 
tional  procedure  Is  summarized  In  Table  5.1. 
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Table  5.1  Summary  of  Incremental -Iterative  Inteoratlon  Procedure 


Initial  Calculations 

1.  Form  the  linear  stiffness  matrix  [K^]  and  mass  matrix  fM];  Initialize 
{%},  {°u},  {°u};  form  the  field  matrix  [s]. 

2.  Calculate  the  time-intepratlon  constants  and  converqence-tolerance 

constants.  Newmark  method:  e*1.0,  a-0.5,  oe0.25 

a0  =  l/(«At2)  a^  =  A/(aAt)  a?  =  l/(riAt)  a3  =  1 /-2  *  1 

a^  =  5/n  -1  a^  =  At( {/a  -?)/2  a^  =  ao  a^  =  -a,, 

a?  =  -a3  aq  =  At(l-A)  alf)  *  SAt. 

Force  tolerance  :  ftol  3  0.1  (unless  otherwise  noted) 

Eneray  tolerance  :  etol  =  0.001  (unless  otherwise  noted) 

3.  Form  the  effective  linear  coefficient  matrix  [K]  =  [K^]  +  ao  [M], 


For  Each  Time  Step 

1.  Calculate  the  maanetlc  field  at  the  new  position  and  form  the  field 
vector  {t+AtP};  calculate  the  vector  of  nodal  point  forces  {*F}. 

2.  Solve  the  eddy  current  eouatlon  (5.8)  and  calculate  the  Induced 


current  and  force  on  the  plate. 


3.  Integrate  and  transfer  the  maanetlc  force  to  the  structural  prooram 
and  form  the  load  vector  (t+AtR). 

4.  Update  [K]  for  nonlinear  stiffness  effects  to  obtain  ftK]  * 

[K1  +  and  trlanoularlze  (*£}:  [**]  «  f L)  fn]  n_]T. 

5.  Form  the  effective  load  vector  {t+AtP.}: 

{t+AtRl  *  {tR}  +  e  ({t+AtR>  ,  {tR})  +  [M]  (#1  +  a?  { tu> )  -  {V}  . 
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6.  Solve  for  the  displacement  increment  {u}. 

7.  If  required,  iterate  for  eciulllbrlum  as  follows: 

(a)  Initialize:  j=0,  {u^}=  {ul 

(b)  j  =  i+1 

(c)  Interpolate  the  new  displacement  and  velocity  and  transfer  to 
eddy  current  program, 

(d)  calculate  the  magnetic  field  at  the  new  position  at  the  same 
time  step  and  solve  for  the  maonetic  force, 

(e)  Integrate  and  transfer  the  maonetic  force  to  the  structural 
prooram  to  form  the  unbalanced  load  vector  {a  Pvu1 


(f)  solve  for  the  j-th  correction  to  the  displacement  increment 
{A  UW>  and  update  the  displacement  increment 
{u(j)}  *  +  {A  u(j)} 

(o)  check  iteration  converoence:  ||  (a  t+Atq(  i) }  ||  ^  j\\  { R}  ||  ^  <  ftol , 

and/or  {A  u(^}T{a  )}  j  {A  uO  )}T{t+At£(l )}  <  eto1  f 

(h)  if  convergence  {u}  *  {u^l  and  go  to  8; 

If  no  convergence  and  limit  on  number  of  Iterations  not  exceeded, 
go  to  (b);  otherwise,  stop  the  program  with  a  messaoe  to  restart 
using  a  smaller  time  step  size. 


8,  Calculate  new  displacements,  velocities,  and  accelerations 
{t+Atu}  =  ag  {ul  +  a?  { tul  +  a8  {ul 

(t+Atul  =  (*u>  ♦  a,  {K}  ♦  a10  (t+Atul 


{t+Atu>  .  {tu}  +  {u} 

9.  If  required  number  of  time  steps  have  been  calculated,  stop;  otherwise, 
go  to  1. 


186 


A  Fortran  program  BEAMNL  has  been  developed  for  the  mechanical  part 
of  the  problem  based  on  the  formulation  and  computational  procedure 
presented.  Some  geometrical 1y  nonlinear  static  problems  have  been  solved 
successfully  uslno  the  proaram  BEAMNL,  but  the  Iteration  procedure  for 
geometrically  nonlinear  dynamic  problems  experienced  some  converoence 
difficulties.  This  converoence  problem  and  some  limited  results  obtained 
are  discussed  and  presented  in  the  next  section.  The  present  version  of 
the  program  BEAMNL  is  briefly  described  below. 

BEAMNL  Is  designed  to  analyze  linear  and  oeometrically  nonlinear 
static  and  dynamic  problems  of  a  linear  elastic  beam  plate.  The  type  of 
analysis  may  be  specified  by  the  input  parameter  NANTY.  For  nonlinear 
problems,  the  Newton  Cotes  quadrature  formula  is  used  to  calculate  the 
element  stress  vectors  {*F}.  The  order  of  the  quadrature  formula  may  be 
specified  in  the  input  data.  One  load  case  for  a  sinole  type  of  analysis 
is  handled  at  each  run. 

Geometric,  material,  and  load  information  need  to  be  input  to  the 
program.  The  number  of  Incremental  solutions  and  number  of  Iterations 
allowed  must  also  be  specified  for  the  nonlinear  problems.  Time  inteora- 
tion  constants  and  convergence  parameters  also  need  to  be  input.  The 
initial  conditions  may  be  specified,  if  any,  (input  data  IGIVN^O)  for  the 
dynamic  problems. 

Displacements,  velocities,  and  accelerations  are  produced  as  output 
of  the  program.  The  number  of  iterations  at  each  incremental  step  Is 
also  printed  out  for  reference.  No  plotting  capabilities  are  incorporated 
yet. 


The  linear  part  and  nonlinear  static  part  of  the  prooram  RFAMNL  have 
been  tested  and  verified  throuqh  comparisons  with  analytical  and  experimen¬ 
tal  results.  One  of  the  nonlinear  static  problems  studied  is  the  larqe 
deformation  analysis  of  a  cantilever  subjected  to  an  end  moment,  shown  in 
Figure  5.3  (Ref.  87).  Five  beam  elements  are  used  for  this  problem  with 
Newton  Cotes  formula  of  order  three  in  both  the  lenath  and  the  thickness 
direction  of  the  element.  Seventy  incremental  steps  are  used  with  eouili- 
brium  iteration.  The  force  and  eneray  tolerance  used  are  0.1  and  0.001, 
respectively.  The  calculated  response  compares  well  with  the  analytical 
solution  within  the  range  of  the  moment  parameter  n  shown  In  the  fiqure. 

Rood  agreement  with  the  analytical  result  Is  expected  for  even  larqer 
values  of  n.  Note  that  In  this  problem  there  is  no  couplfm  between  the 
axial  force  and  the  flexural  deformation  since  It  is  a  Dure  bendinq  problem. 

The  second  problem  studied  is  the  larqe  deformation  analysis  of  a 
cantilever  subjected  to  an  end  load,  Fioure  5.4.  Seven  beam  elements  are 
used  with  the  mesh  finer  In  the  region  closer  to  the  support.  The  deflections 
of  the  beam  at  different  load  levels  are  shown  In  Fiqure  5.4.  The  curva¬ 
ture  at  the  root  of  the  cantilever  beam  Is  shown  as  a  function  of  load  In 
Fiqure  5.5,  The  theoretical  and  experimental  results  are  from  an  unpub¬ 
lished  study  by  Pao  and  Moon  of  Cornell  University.  The  FF  analysis 
predicts  deflections  accurately  for  tip  deflections  up  to  about  60?-  of 
the  cantilever  span,  corresponding  to  a  loadinq  parameter  f  of  nearly  2.0. 

The  curvature  Is  accurate  up  to  about  f  *  4.0.  Note  that  In  the  analytic 
results  the  beam  Is  assumed  Inextensible,  while  In  the  FE  results  the 
effect  of  axial  deformation  Is  included. 


Displacement  ratios,  u/i,  v/Jt,  <f>/2-rr 
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L  z  1 00  in 

I  =  0.01042  in 

A  =  0.5  in2 

E:  1 .2  x  I04psi 

M  =  Concentrated  end 
movement 


Moment  parameter  17  =  MJt/2rrEI 


Figure  5.3  Moment  deflection  curve  of  an  elastic  cantilever 
plate  loaded  by  a  moment  at  its  tip. 


—  Theoretical 
a  FEM  (BEAMNL) 


t  *  x  y  r/ei 


4  Comparison  of  F£  and  theoretical  predictions  of  the 
nonlinear  elastic  bending  deflections  of  a  tip-loaded 
cantilever  beam. 


s 


Figure  5.5  Comparison  of  theoretical,  experimental,  and  FE 
results  for  the  curvature  at  the  root  of  the 
cantilever  beam  shown  in  Figure  5.4. 
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Strono  couplinq  between  the  axial  force  and  the  flexure  deformation 
exists  in  this  problem.  When  the  loading  parameter  f  becomes  oreater  than 
?.5,  the  behavior  of  the  deflected  beam  becomes  more  like  an  arch.  Slow 
convergence  of  the  iteration  scheme  occurs  for  this  situation.  At  laroer 
values  of  f,  convergence  fails  with  the  present  quasi -Newton  method. 
Stronger  convergence  tolerances  with  ftol  =  1x1  f)"*  and  etol  =  lxlfl"fi 
have  been  used  for  this  structural  stiffenino  problem,  and  a  small  time 
step  is  needed  for  convergence  (Ref.  R9). 

Geometrically  nonlinear  analysis  has  been  attempted  for  the  laroe 
amplitude  free  vibrations  of  a  simply  supported  beam.  The  purpose  of  the 
analysis  Is  to  establish  the  frequency-amplitude  relationship  of  the  beam 
and  to  compare  it  with  the  analytical  result  (Refs.°l,  9?).  This  study  Is 
not  successful  because  the  Iteration  scheme  falls  to  converae  within  each 
time  step.  The  numerical  experience  shows  that  the  lack  of  converoence 
arises  from  the  hi qh  frequency  axial  vibrations  of  the  beam  plement.  The 
larqe  axial  forces  developed  from  these  axial  vibrations  cause  an  increase 
of  the  unbalanced  force  during  the  Iterative  procedure.  Attempts  have 
been  made  to  separate  the  inertia  force  effect  and  the  effect  of  the  linear 
approximation  to  the  nonlinear  equations  in  the  calculations  of  the  ur - 
balanced  forces,  but  with  the  present  quasi -Newton  scheme  all  the  attempts 
failed  because  of  the  intrinsic  Interaction  of  the  time  intearatlon  method 
with  the  approximate  solution  procedure  for  the  eauatlons  of  motion. 

A  quasi-Newton  method  and  structural  (i.e.,  one-dimensional)  beam 
elements  are  used  in  the  nonlinear  dynamic  analysis  In  this  work.  A 
literature  survey  reveals  that  for  transient  nonlinear  analysis  all  the 
works  surveyed  (Refs.  Q3  -  96)  whicl.  use  continuum  (i.e.,  two-dimensional 
plane  stress)  elements  adopt  either  a  quasi -Newton  or  the  full  Newton 
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method  for  Iteration.  Transient  nonlinear  analyses  uslna  structural  elements 
usually  employ  explicit  time  inteoratlon  schemes  (Ref.  °7,  98).  Studies 
which  use  structural  elements  tociether  with  an  implicit  Inteoratlon  scheme 
and  quasi -Newton  method  for  Iteration  have  not  been  found. 

In  transient  nonlinear  analysis  the  dynamic  properties  of  different 
types  of  elements  have  a  definite  effect  on  the  success  of  a  particular 
solution  method  used.  Efforts  to  use  Implicit  Inteoratlon  schemes  with  a 
quasi-Newton  method  for  Iteration  appear  to  be  futile.  It  is  suqaested 
that  the  full  Newton  or  the  BFGS  method  (Ref.  99)  be  used  for  the  Iterative 
procedure  In  the  future. 


rfvw-v 
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Chapter  6 
CONCLUSION 


This  thesis  has  presented  an  Integrated  study  of  the  magnetically 
Induced  vibrations  of  nonferrous  conducting  plates.  The  continuum 
mechanics  basis,  the  stream  function  modelling  of  the  coupled  problem, 
the  Intermediate  eddy  current  calculations  on  rigid  flat  plates,  and 
the  linear  vibrations  of  Infinitely  long  plates  have  been  discussed  In 
some  detail.  Some  preliminary  studies  of  the  coupled  nonlinear  problems 
of  Infinitely  long  plates  have  also  been  presented.  A  summary  of  the 
work  presented  In  this  thesis  Is  given  below.  Conclusions  from  the 
present  study  are  then  drawn.  Some  suggestions  for  further  research 
are  made  at  the  end  of  the  chapter. 

6.1  Summary 

The  thesis  contains  three  main  phases:  theoretical  modelling  of  the 
magnetoelastic  plate,  finite  element  eddy  current  calculations  on  rigid 
plates,  and  finite  element  numerical  studies  of  coupled  magnetoelastic 
problems.  The  numerical  results  obtained  have  been  compared  with  analytic 
and  experimental  results  at  each  stage  of  the  study.  The  physics  and 
continuum  mechanics  background  of  the  problem  has  been  presented  In 
Chapter  Two.  The  magnetic  force-field  method  has  been  chosen  to  formulate 
the  coupled  problems.  The  magnetic  force-energy  method  (Ref.  25)  has  not 
been  emphasized,  although  the  energy  transfer  between  the  magnetic  and 
mechanical  subsystems  has  been  briefly  mentioned  In  Chapter  Two,  and  the 
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magnetic  energy-circuit  method  for  eddy  current  calculation  discussed  In 
Chapter  Three. 

The  modelling  of  plate  problems  by  the  stream  function  method  has 
been  presented  In  Section  2.3.  All  the  simplifying  assumptions  have  been 
carefully  examined  and  the  limitations  of  the  stream  function  method  In 
the  large  deformation  cases  analyzed.  The  governing  equations  for  the 
magnetic  and  mechanical  subsystems  have  been  derived  and  the  various 
coupling  effects  discussed.  The  assumption  of  constant  magnetic  body 
force  across  the  thickness  of  the  plate  and  the  neglect  of  the  magnetic 
damping  moment  have  also  been  discussed.  Both  linear  and  nonlinear  prob¬ 
lems  have  been  studied.  All  the  equations  have  been  linearized  to  the 
first  order  of  the  various  unknown  variables.  The  updated  laqranglan 
description  of  the  magnetic  subsystem  and  the  approximations  Involved  In 
the  linearization  procedure  have  been  given  special  attention. 

Eddy  current  calculations  on  rigid  conducting  plates  have  been  pre¬ 
sented  In  Chapter  Three.  The  EE  Galerkln  formulations  and  the  treatment 
of  the  Integral  terms  In  the  eddy  current  equations  have  been  presented 
In  detail.  Comparisons  of  the  numerical  results  with  the  Infrared  and 
search  coll  experimental  results  have  been  made.  Analytical  solutions 
In  the  low  and  high  magnetic  Reynolds  number  limits  have  also  been  employ 
ed  to  verify  the  numerical  results.  Detailed  studies  Including  some 
dimensional  analysis  have  been  carried  out  for  the  steady  state  and 
transient  eddy  current  problems  on  Infinitely  long  plates.  For  finite 
plates,  only  the  steady  state  analysis  has  been  formulated  and  performed. 

Chapter  Four  treats  the  linear  vibrations  of  Infinitely  long  conduct 
Ing  plates  excited  by  single  pulsed  currents.  The  types  of  coupling 
studied  have  been  limited  to  the  effects  of  the  transverse  motions  of  the 
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plates.  The  staggered  transient  analysis  and  the  computational  procedure 
have  been  presented  In  detail.  The  calculated  moments  have  been  compared 
to  experimental  values  obtained  with  low  Inductance  strain  gage  readings. 
The  effect  of  pulse  duration  on  the  Induced  vibrations  and  energy  transfers 
have  been  demonstrated  through  the  three  pulsed  currents  studied.  The 
results  show  that  as  the  pulse  duration  decreases,  the  nonlocal  effect 
becomes  more  Important  and  the  magnetic  force  acts  more  like  an  Impact 
load.  A  quantitative  description  of  the  effect  of  pulse  duration  on  the 
Induced  vibration  has  not  been  presented  because  the  control  condition  of 
the  driving  currents  does  not  readily  permit  equivalence  of  the  Input 
energies  In  the  three  problems  studied.  A  rough  method  has  been  suggested 
to  estimate  the  effect  of  pulse  duration  on  the  contribution  of  nonlocal 
pushing  force  In  the  mechanical  energy  converted.  The  effect  of  the 
difference  In  the  time  characteristics  of  the  two  subsystems  on  the  Induced 
vibration  may  then  be  represented  quantitatively. 

The  Induced  nonlinear  vibrations  of  Infinitely  long  plates  have  been 
formulated  In  Chapter  Five  based  on  the  linearized  equations  for  the  up¬ 
dated  Incremental  analysis  procedure.  The  staggered  transient  analysis 
and  the  computational  scheme  have  been  presented.  An  Iterative  procedure 
has  been  proposed  for  the  motion-dependent  magnetic  load.  Some  preliminary 
results  for  the  nonlinear  static  problems  have  been  presented.  The 
convergence  difficulties  In  the  nonlinear  dynamic  problems  have  been 
discussed  and  some  suggestions  made  for  the  Improvement  of  the  Iteration 
scheme  used. 

6.2  Conclusions 

The  following  conclusions  may  be  drawn  from  the  work  presented: 


***** 
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1.  The  stream  function  method  presents  a  useful  tool  for  the  eddy  current 
calculations  on  thin  nonferrous  conducting  plates  In  the  low  and  Intermedi¬ 
ate  frequency  cases.  The  capability  of  this  method  for  curved  surfaces  Is 
limited  to  small  curvature  situations.  This  method  has  a  close  relation¬ 
ship  with  other  eddy  current  circuit  models.  The  advantages  and  dis¬ 
advantages  of  the  different  methods  need  further  comparative  studies. 

2.  From  the  experience  of  using  the  Fortran  programs  developed,  the  E0DY1 
and  EDDY1T  codes  appear  to  be  reasonably  efficient.  The  EDDY2  code  takes 
more  computation  time  because  of  the  complicated  forms  of  the  weighting 
functions  for  the  six-node  triangular  elements.  The  efficiency  of  the 
EDDY2  code  can  be  greatly  Improved  by  using  rectangular  elements  and/or 
linear  triangular  elements.  The  weighting  functions  for  these  elements 
can  be  easily  Integrated  or  obtained  from  the  welqhtlng  functions  for  the 
six-node  triangular  elements  presented  In  Appendix  B. 

3.  The  modelling  and  formulations  for  the  linear  coupled  problem  In 
Chapter  Four  have  been  verified  satisfactorily  through  comparisons  with 
low  Inductance  strain  gage  experiments.  The  velocity  and  displacement 
have  a  damping  effect  on  the  coupled  problem.  The  total  energy  transferred 
to  the  conducting  plate  Is  reduced  because  of  these  motional  effects,  but 
their  Influence  seems  to  be  small  for  the  problems  studied.  The  pulse 
duration  has  an  effect  on  both  the  magnetic  force  and  the  Induced  vibra¬ 
tions.  For  the  short  pulse  an  Impact  approximation  based  on  the  nonlocal 
pushing  force  will  be  reasonable. 

The  quasi -Newton  method  cannot  be  used  In  the  Iterative  procedure  together 
with  the  structural  elements  In  the  transient  nonlinear  problems.  The 
fictitious  In-plane  deformation  caused  by  the  simplifying  assumptions  In 
the  structural  theory  creates  large  In-plane  Inertia  forces  which  make 
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the  convergence  of  equilibrium  Iteration  difficult  when  the  quasi-Newton 
method  Is  used.  The  simplifying  kinematic  assumptions  In  the  structural 
elements  have  a  definite  Influence  on  their  performance  In  transient  non¬ 
linear  analysis.  The  full  Newton  or  a  BFGS  method  should  be  used  for 
these  problems  If  structural  elements  are  to  be  used.  The  study  of  the 
Induced  nonlinear  vibrations  of  conducting  plates  cannot  proceed  until 
these  numerical  problems  are  solved. 

6.3  Suggestions  for  Further  Research 

The  following  suggestions  are  made  for  future  work: 

1.  For  the  Immediate  continuation  of  the  present  work.  It  Is  suggested 
that  linear  elements  be  Implemented  In  the  program  FDDY2  to  achieve 
better  efficiency  of  the  computation;  that  the  effect  of  pulse  duration 
on  the  energy  transfer  and  Induced  vibration  of  Infinitely  lono  plates  be 
studied  In  more  depth;  and  that  the  full  Newton  and  BFGS  methods  be  tested 
for  the  transient  nonlinear  problems  to  enable  the  continuation  of  the 
study  for  the  Induced  nonlinear  vibrations  of  conducting  plates. 

2.  The  eddy  current  circuit  nature  of  the  stream  function  formulation 
should  be  explored  to  enable  a  modal  analysis  of  the  eddy  current  and  the 
coupled  problems.  The  magnetic  energy  stored  In  various  eddy  current  modes 
and  Its  conversion  Into  mechanical  energies  In  different  vibrational  modes 
of  the  plate  should  be  calculated  In  such  a  modal  theory  for  the  coupled 
problems.  The  modal  study  of  the  coupled  problem  described  basically  Is 
equivalent  to  the  magnetic  force-energy  method  discussed  In  Reference  25. 

3.  The  In-plane  force  and  the  magnetic  damping  moment  should  be  Included, 
and  their  effect  on  the  Induced  vibrations  and  stability  of  the  conductors 
studied.  Different  types  of  magnetic  fields,  such  as  transverse  or  Inclined 


uniform  fields,  should  also  be  Investigated. 

4.  Cyclically  pulsed  currents  and  their  effect  on  Induced  vibrations 
should  be  studied  for  possible  magneto- flutter  effects.  Heat  conduction 
and  thermoelastic  coupling  In  such  cases  may  be  Important  and  may  need  to 
be  Included  In  the  analysis. 


APPENDIX  A 


MATRICES  AND  NONLOCAL  INTEGRATIONS  FOR  ONE-DIMENSIONAL  PLATES 


This  appendix  presents  explicit  finite  element  matrices  including  the 
analytic  expressions  for  the  nonlocal  Integration  terms  described  In  sub- 
subsection  3.3.1.b  for  the  transient  one- dimensional  problems.  The 
matrices  and  expressions  for  the  steady  state  problems  described  In  sub¬ 
subsection  3. 3.1. a  may  be  obtained  from  these  by  setting  h  *  1,  ap  «  R, 
and  by  replacing  X  and  X'  by  the  nondlmenslonalized  quantites  x  and  £. 


1  -1 

-1  1 


(A.l ) 


[PE]  »  op  LE 


1/3  1/6 
1/6  1/3 


(A.2) 


*  {ln  r(X'-Xl>2  +  }  h2]  -  An  [(X'-XE)2  +  {  h2]}  + 
-  ^  (X*-x|)  {tan-1  [£  (X’-XE)]  -  tan"1  [jj-  (X'-XE)]> 


(A. 3a) 


W2  ”  "  {ln  f(X'-Xl)2  +  J  -  A"  [(X'-x|)2  +  |  h2]}  +  (A. 3b) 

♦  ^  (X'-XE)  (tan"1  [|  (X'-XE))  -  tan"1  (X*-x|)]> 

In  which  the  superscript  E  denotes  the  Eth  element,  LE  the  length  of  the 
element,  and  XE  and  XE  the  coordinates  of  the  two  end  nodes  of  element  E, 
as  shown  In  Figure  A.l. 
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The  nonlocal  Integration  Eq.  (3.41)  Is  performed  element-wise  and 
calculated  over  all  the  elements  for  each  element  E.  For  any  element  G, 
the  contributions  to  the  nonlocal  Integration  of  element  E  are.  Figure 
A.l, 

X2 

Qjk  '  U  N®  (X1)  WE  (X*)  dX\  j,  k-1,2  ( A. 4) 

G 

*1 

The  resulting  analytical  expressions  of  these  Integrations  are: 

"i i  ■  Jr  ^  ^  * 

♦  F  [-  U?-xE)2  +  (xG-xE)2  +  <XG-XE)2  -  (x|-x|)2]  + 

+  7  (XG-XE)2  -  2  (Xf-Xf)  (XG-XE»  tn  [  (xf-X^>2  +  + 

+  7  [‘T7'  (X?-^)2  ♦  2  (x|-x|)  <XG-XE)1  tn  f (X^-x|)2  +  J2]  + 

+  I  ['T7+  (XrXl)2  -  2  LE  (Xf-Xf)l  tn  [(XG-XE)2  ♦  fl  + 

+  7  [T7  "  ( X2” X^ ) 2 ]  tn  r(x|-x|)Z  +  J2]  + 

+  (|  (X?-XE)3  +  2  (LE+  LG)  (XG-XE)2  -  4  LE  (XG-XE)(XG-XE)  - 
-  j-  (x|-x|)]tan_1  (J)  (XG-XE)  + 

+  [-  |  <X?-XE)3  -  2  LG  (XG-XE)2  ♦  (xf-x|)l  tan’1  (|)  (XG-XE)  + 
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♦  [|  (X|-X^)3  -  2  (X|-X^)2  <x|-x|)  +  f  (x|-x|)]  tan-1  (£)  (x|-X^)  + 

♦  [§  <XG-XE)3  -  f  (x|-x|)]  tan'1  (£)  (x|-x|)  } 


1  OH  /  ol,  E.  6 
K  jT[S  t  -  2hl  L  4 


£  [(X?-Xf)2  -  (XG-XE)2  -  (x|-x^)2  4  ( x|- x| ) 2 ]  + 


7  [-  yjr  ♦  (x?-xE)2  -  2  le  (xg-xe)]  m  t(xG-xE)2 


+  7  [T7  •  (X^-X|>21  m  [(x^-x|)2  +  J2]  + 

+  7  [T?  +  <x2"xl)2  -  2  (xf-xf )(x|-x|)]  in  [ (x|-xf )2  4  J2]  + 

+  7  T7  '  (X2'X2)?  +  2  (X2G^><X?-XI>U"  t(x2-x2>2  +  J]  + 


+  [-  \  (xf-XE)3  4  2  lE  (XG-XE)2  4  £  (XG-XE)]  tan-1  (£)  (XG-XE)  4 
+  [|  <XG-XE)3  -  ^  (Xf-Xf)!  tan’1  (£)  (X^-x|)  4 
♦  (|  (x|-xf)3  4  2  (lE-lG)  (XG-XE)2  -  4  LE  (XG-XE)(XG-XE)  - 
-  J-  (xf-x|)l  tan'1  (£)  (XG-XE)  4 
4  [|  (XG-XE)3  -  2  (XG-XE)2  (XG-XE)  4  £  (Xf-X|)]. 

*  tan'1  (|)(XG-XE)  } 


k 
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»2?  ■  h  ^  ^  * 

♦  jf  [(xf-x^)2  -  (xf-x|)2  -  (x|-xf )2  +  (x|-x|)2]  + 

♦  \  t-  4  -  (xrxi>2  +  2  (xrx5)(x2“xi)i  tn  f(xrxi)2  +  r1  + 

♦  7  r^-+  (Xf-xf)2  -  2  (xf-x^){x|-x|)]  tn  [(xf-x|)2  +  Jp]  + 

+  \  -  (4-xf)2!  tn  [(XG-XE)2  +  Jp]  + 

+  J  r-  ^  +  (XG-XE)2  +  2  LE  (x|-x|)1  tn  «XG-xf)2  +  Tp]  + 

♦  [j  (X?-xf)3  -  2  (xf.xf)2  (xf.xf)  +  {£  (x^-x^ltan"1  (|)(X^-X^)  + 

+  [|  (Xf-X|)3  -  2  (LE-LG)  (XG-XE)2  +  4  LE  (X^-XE)(XG-XE)  - 

-  (XG-XE)]  tan'1  (£)  (XG-XE)  + 

+  [f  (XG-XE)3  -  tp  <xi~xi >1  (|)  (x2"xi)  + 

+  [-  \  (XG-XE)3  -  2  LE  (XG-XE)2  +  (Xf-Xf)l  tan'1  (§)  (x|-x|)  } 


0 


«!!  ■  h i ^  * 

+  J  [-  (XG-XE)2  +  (XG-XE>2  +  (Xg-Xf)*  -  (x|-x|)2]  + 

+  7  [t?  •  (xrxi)2]  ln  [(xrxi)2  +  ir]  * 

*  7  [-  ^  -  (xf-x|)2  +  2  LE  (x®-x|)]  in  [(xf-x|)2  +  5^]  + 

+  J  [-  ^  -  (XG-XE)2  +  2  (XG-XE)(XG-XE)]  tn  [(xf-xf)2  +  Jh  + 

+  7  [77  +  ^x2"x2)2  ■  2  (X1*X2)(X2‘X1)]  ln  KX|-X1)Z  +  r-3  + 

+  [f  (XG-XE)3  -  (XG-XE)]  tan"1  (|)  (XG-XE)  + 

+  [f  (XG-XE)3  -  2  (XG-XE)2  (XG-XE)  +  (X?-XE)]. 

•  tan"1  (^)(xf-x|) 

+  [-  |  (x|-X^)3  +  2  LG  (XG-XE)2  +  (XG-XE)]  tan"1  (£)  (XG-XE)  + 
+  t|  (X2G-X|)3  -  2  (LE+LG)  (XG-XE)2  +  4  LE  (XG-XE)(XG-XE)  - 
-  (xf-XE)]  tan"1  (|)  (x|-x|)  } 


This  appendix  presents  the  finite  element  matrices  and  the  analytic 
expressions  of  the  nonlocal  weighting  functions  described  In  subsection 
3.4.1  for  the  steady  state  two-dimensional  eddy  current  problems.  The 
nondl mens  1 onal 1  zed  finite  element  equations  are  those  presented  in  Eq. 
(3.53). 

By  using  the  natural  coordinates  L,  l_2  l-3  for  the  six-node  tri¬ 
angular  element  shown  In  Figure  B.l,  one  may  write  the  shape  functions  N, 
as: 


N,  »  1,(21.,-?) 

N4  “  4Ll4 

N2  "  L2(2L2‘1J 

N5  ”  4L2L3 

(B.l) 

N3  -  L3(2L3-1) 

N6  “  4L3L1 

In  which 


1 

TK 


al  bl  cl 
a2  b2  c2 
*3  b3  c3 


(B.2) 


ai  a  xi+i  yi+i  •  xi+2  *i+i 
b1  •  *1+1  -  *1*2 


1  -  1,2,3 


(B.3) 


C1  ’  x1+2  ’  x1+l 


— -*■  V  iii  4 
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Figure  B.l  Six-node  triangular  element  in 
global  coordinates. 


Element  E 


4 -  Q  - * 


Figure  B.2  Translated  and  rotated  coordinates 
used  in  nonlocal  integration. 
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A  *  \  (X2Y3  <  X3 Y,  +  X1V2  -  X^  -  X3Y2  -  X^)  (B.4) 

The  matrices  [K]  and  [P]  In  Eq.(3.53)  may  then  be  integrated  to  obtain 
6 


r(41  wRA 

tn  mw 


-1 

6 

syrnn, 

-1 

-1 

6 

0 

0 

-4 

32 

-4 

0 

0 

16 

32 

0 

-4 

0 

16 

16 

32 


(B.5) 


and 


0 

'  3  ^2^3  +  C?C3^ 

-  K  C h2b3 +  c2c3^ 

-  1  (b,b3  +  c^-j) 

0 

-  j  (bjbj  +  c^) 
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To  Integrate  the  weighting  functions  W;j  In  Eq.  (3.55),  the  element 
Is  referred  to  a  new  coordinate  system  (s,  t)  obtained  by  the  following 
relationship.  Figure  B.2, 


w. 

M 


cos  e  sin  e 
-sin  e  cos  e 


\  x-*'  ! 
i  y-yC)  \ 


(B.7) 


cos  9  -sin  9  ^  s  I  )  X]  I 

M  sin  9  cos  9  (  t  \  l  yE  ^ 

Under  this  transformation, 

( ( x- c ) 2  +  (y-n)2  +  1)3/2  *  [(s-s*)2  ♦  (t-t')2  +  l]3^2 
-  1 -u- ( 1-c) V ,  l2  -  u-cv,  L3  ■  V 
In  which 

u  ■  s/a,  v  ■  t/b,  c  *  s3/a 

The  six  weighting  functions  W*j  may  then  be  written  as 

*1  -  ab  II,  -  3  Iu  -  3<l-c)  Iy  +  2  Iuu  ♦  4< 1-c)  Iuy  ♦  2(1-c)2  Iyv) 


w!  ■  ab  [-  I  ♦  c  I  +  2  I  -  4c  I  +  2c2  I  ) 
2  1  u  v  uu  uv  vv 


(B.8) 


(B.9) 


(B. 10) 


WE3  -  ab  [-  Iy  ♦  2  Iyv] 


«4  »  4ab  IIU  -  c  Iv  -  Iuu  ♦  (2c*1)  Iuy  ♦  c(l-c)  Ivy] 
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K3  =  in 


/ a2(u‘-1)2-*-b2v,^l  -  a(u'-l) 
/  a2u'2+b^v,2+l 


-  au 


(B.20) 


K4  -  /a^(u,-c)2+b2{v'-l)2+1  -  / a^u'*+b*v'^+l 


K  -  / a*(l-u‘)2+b*v,Z+l  -  / aVc-u'^b^v'-l)2*! 

b 


(8.21 ) 

(B.22) 


The  nonlocal  matrix  [Q]  Is  calculated  numerically  by  using  these  six 
analytic  expressions  of  the  weighting  function  W*j  within  each  element. 

Pure  numerical  integrations  are  used  to  evaluate  the  nonlocal  Integrations 
for  (u',v')  outside  the  element. 


APPENDIX  C 

MATRICES  AND  VECTORS  FOR  NONLINEAR  DYNAMIC  ANALYSIS 


This  appendix  summarizes  the  matrices  and  vectors  for  the  finite 
element  UL  transient  nonlinear  analysis  described  In  Section  5.1.  Details 
of  the  derivations  are  given  by  Bathe  and  Bolourchi  (Ref.  87). 

The  Interpolation  functions  for  the  components  of  displacement 
Increment  In  Eg.  (5.1)  are 


(C.l) 


!-  £  6^ 

"*♦2  f 

-6^i  S*3 

h\ 

o  <t»4 

L*5  0 

*6  "^1 

In  which 

{u}  -  |U1  u2  u3  u4  u5  u6j 


(C.2) 


*1  *  f  “  (P  » 


*2  *  1  -  4  I  +  3  (£) 


<f>-j  *  2  f  -  3  (£)  ,  i>,  *  1  -  3  (£)  +  2  (£) 


( C.  3) 


*5  ■  f  -  2  (f)Z  ♦  (£)\  +6  -  3  (f)2  -  2  (f)3  ; 

and  In  which  L  Is  the  length  of  the  beam  element,  r,  s  are  the  beam  convect- 
ed  coordinates,  Figure  C.l.  The  shear  deformation  Is  not  Included  In  this 
approximation. 

The  linear  strain-displacement  transformation  matrix  [tBL]  relates 
the  linear  part  of  the  strain  components  to  the  nodal  degrees  of  freedom 
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Figure  C.l  Convected  coordinates  and  degrees  of 
freedom  of  beam  element. 


Figure  C.2  Current  and  original  local  coordinates 
of  beam  element. 


f  ^ 


?.]  7 


r  [tB|  ]  {ul,  (C.4) 

( ) 

and  Is 


[\] 


1  i9rs  al  aII  1  c!_  ,s  crs 

■r  4l-6[7  r  '6[7+1z[3  2r-^ 


(u) 


(C.5) 


0  0  0  0  0  0 

The  linear  stiffness  matrix  [^K^]  calculated  from  Eq.  (5.3)  uslnq 
[*8^]  given  above  and  Is  the  same  as  the  usual  linear  beam  stiffness  matrix: 


[\]  *  E 


A 

r 


121 

17 


61  41 

r 


Symm. 


A 

r 


0  - 


0 

0 

A 

r 

121 

61 

0 

121 

r 

"7 

~7 

61 

21 

A 

61 

7 

r- 

u 

41 

r 


(C.6) 


The  matrix  of  Cauchy  stresses  [^t]  In  Eq.  (5.4)  In  this  case  is 

Tn  0  0 

o  Tl1  0 


(C.7) 


0 


0 


0 


The  nonlinear  strain-displacement  transformation  matrix  [°NL1  Is 


[  bnl]  -  tbi j 33x6 


1  gS  i  o  rs  g  rs  1  /r  s  ,i  ^rs  «s  ^  rs 

r  6[?-12[5  r  _6^  12^  2r_6[? 


0  -6^?  +6\  1  -4f  0  6^?  -6^r  -2f 

l/  LJ  L  L  \r  LJ  L  l/ 

0  _1  +4f  -3[7  0  _6[?  +6jr  2f  '3f? 


(C.8) 


The  nonlinear  stiffness  matrix  [^K^]  In  Eq.  (5.4)  is  obtained  using  [*8^] 
given  above  and  Is  shown  in  Eq.  (C.9).  The  Newton-Cotes  quadrature  formula 
Is  used  to  Integrate  [*K^]  numerically. 


I 


The  vector  of  Cauchy  stresses  {^i)»1n  Eq.  (5.5)  Is 


{tT}  =  LT11  Tn  J 


(C. 10) 


In  which  has  been  set  to  zero.  The  Newton-Cotes  formula  Is  also  used 
to  Integrate  the  vector  of  nodal  point  forces  (V)  given  below. 

I  Tn{“f)  I 


Tn(6f?’12^) 


[tp}  =  ;v  Tn(4r  "6fl) 


(C.ll) 


tn(-6^  +i2^|) 


,H(2f-6l|) 


The  element  stiffness  equations  must  be  transformed  to  global  coordinates 
before  they  are  assembled  Into  a  set  of  global  equations.  The  transformation 


matrix  from  local  to  global  coordinates  Is 


[tR]  »  [**]  [°R]  (C,12) 

In  which  [°R]  Is  the  transformation  matrix  from  the  original  local  coordinates 
of  the  element  to  the  global  coordinates,  and  [tR]  Is  the  transformation 
matrix  from  the  current  to  the  original  local  coordinates  of  the  beam 
element.  [tR]  Is  given  by 


8 


(C. 13) 
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in  which 


tlR] 


cos^  s1n*e  0 

s1nte  coste  0 

0  0  1 


with 
s1ntB  * 


(C. 14) 


«  [ (°L  +  +  (o«2)2]1/2  (C.15) 

as  shown  in  Figure  C.2. 

The  element  matrices  and  vectors  In  global  coordinates  are 
[tR)T[tK][tR],[tR]T{tF>,  and  £tR]T{u>,  respectively.  In  which 

^K]  =  ♦  f1^). 

The  element  stress  Increments  are  calculated  by  using  the  following 
strain  increment 

e,,  *  I  \  u.  +  ~  (tL.t"Ati)  (C.16) 

11  j*2  Llj  3  °L 

The  total  element  stresses  are  updated  by  using 


t+At 


rll 


T11 


♦  E  e 


11 


/ 

J 


(C.17) 
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